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Introduction 


This  report  is  in  support  of  the  Integrated  High  Payoff  Rocket  Propulsion  Technology 
(IHPRPT)  Demonstration  Program,  which  provides  a  technology  roadmap  for  demonstrating 
significant  enhancements  in  mission  capabilities,  cost  efficiencies,  and  operability  for  various 
future  propulsion  needs.  The  scope  of  this  project  is  to  demonstrate  that  it  is  feasible  to 
numerically  construct  diagrams  of  flammability  limits  for  energetic  ionic  liquid  (IL) 
fuel/oxidizer  combinations,  which  will  significantly  aid  AFRL/RQRP  in  its  ability  to  identify 
hypergolic  propellants  and  predict  their  ignition  delay  times,  prior  to  carrying  out  any 
synthetic  or  empirical  screening  activities.  The  execution  of  the  proposed  modeling  and 
simulation  (M&S)  study  early  in  the  life-cycle  of  the  IL  propellant  development  program 
provides  a  mitigation  strategy  to  reduce  the  risk  of  failure  for  the  insertion  of  IL  fuel 
technology  into  the  small  satellite  market  as  AFRL/RQRP  continues  to  meet  IHPRPT  goals. 

Globally  in  the  last  few  decades,  pioneering  fuel  research  has  been  carried  out  in  order 
to  replace  existing  fuels  for  various  economic  and  ecologic  reasons.  This  development  has  not 
just  been  limited  to  auto  fuels  but  has  included  all  types  of  combustion  fuels,  including  space 
fuels.  One  of  the  most  utilized  bipropellant  combinations  for  satellite  propulsion  and 
interplanetary  missions,  monomethylhydrazine  (MMH)/nitrogen  tetroxide  (NTO),  falls  under 
this  scope.  The  mixture  is  very  energetic  and  presents  some  practical  advantages  in  space 
propulsion;  however,  its  components  at  the  same  time  possess  drawbacks.  MMH  and  NTO 
are  toxic  and  carcinogenic  compounds  and  MMH  vapors  are  flammable.  One  of  the 
alternatives  proposed  to  replace  MMH  are  ionic  liquid  fuels. 

Ionic  liquids  exhibit  physico-chemical  properties  (of  low  vapor  pressure,  good  thermal 

stability,  etc.)  that  make  them  suitable  for  replacing  MMH  while  also  maintaining  their 

hypergolicity  with  propellant  oxidizers  (such  as  nitric  acid  and  NTO).  Due  to  the  large 

18. 

number  of  ionic  liquids  (10  )  possible,  and  the  impossibility  for  scientists  to  test  them  all,  it 
is  necessary  to  determine  a  simple  physico-chemical  criterion  (like  the  adiabatic  flame 
temperature)  to  predict  mixture  flammability  limits  to  model  hypergolic  conditions,  as  well  as 
to  understand  the  chemical  risks  associated  during  the  handling  of  IL  propellants. 

The  part  of  this  study  is  divided  into  three  sections.  In  the  first  section,  a  bibliographic 
search  is  presented  and  a  determination  of  the  flame  temperature  criterion  is  also  discussed. 
In  the  second  section,  the  thennodynamic  data  for  these  ionic  liquids  and  the  methods  used  to 
obtain  the  data  is  discussed.  The  computed  results  obtained  for  ionic  liquid  +  oxidant  mixtures 
are  reported  in  the  third  section  of  the  report. 
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1  State  of  the  art 

In  most  cases,  liquids  or  gases  once  mixed  do  not  react  violently  at  room  temperature. 
However,  some  combinations  are  known  to  react  rapidly  and  under  appropriate  mixing  ratios 
can  lead  to  violent  explosions.  Monomethylhydrazine  and  nitrogen  tetroxide  (MMH  &  NTO) 
is  one  such  combination  and  is  the  preferred  choice  for  creating  propulsive  thrust  in  numerous 
spacecraft  applications.  These  low  temperature  (at  298  K  or  below)  self-  or  auto-ignitions  are 
often  referred  to  as  hypergolic  ignitions.  In  addition  to  the  MMH  &  NTO  combination, 
MMH,  hydrazine  and  unsymmetrical  dimethyl  hydrazine  have  also  been  used  with  nitric  acid. 
It  is  crucial  to  be  able  to  predict,  a  priori ,  the  auto-ignition  behavior  that  may  arise  from  the 
mixing  of  two  or  more  liquids,  either  for  the  development  of  new  chemical  rocket  propulsion 
technologies  or  for  the  purpose  of  mitigating  the  chemical  risks  associated  in  the 
transportation  and  management  of  rocket  fuels.  An  extensive  search  in  the  open  literature  has 
shown  that  there  are  no  computational  methods  currently  available  to  address  these  needs. 
The  objective  of  this  research  is  to  present  a  procedure  to  predict  mixture  flammability  limits 
to  model  hypergolic  conditions,  as  well  as  to  understand  the  chemical  risks  associated  when 
handling  hypergolic  propellants. 


1.1  Ionic  Liquids 

Since  the  beginning  of  the  early  2000s,  ionic  liquids  have  emerged  as  important 
classes  of  compounds  for  many  chemical  applications,  especially  as  industrial  solvents. 
Smiglak  and  co-workers  have  recently  reviewed  these  applications  and  determined  that 
>8,000  studies  have  been  published  on  ionic  liquids  [Smiglak  M.  et  ah,  2007]. 


1.1.1  Definition 

Salts  are  compounds  formed  by  the  association  of  ions.  Ionic  liquids  are  salts  with 
melting  points  equal  to  or  below  100°C.  Room  temperature  ionic  liquids  (RTILs)  are  liquid 
salts  at  ambient  temperatures  (20  -  25°C).  The  first  RTIL,  ethylammonium  nitrate,  was 
described  in  1914  by  Walden  and  has  a  melting  point  equal  to  12°C  [Pucheault  M,  2010]. 

Ionic  liquids  are  synthesized  by  the  reaction  of  a  base  with  a  strong  acid.  Ionic  liquids 
consist  of  negative  and  positive  ions.  It  is  generally  accepted  that  the  anion  (negative  ion)  is 
more  important  in  controlling  the  physical  properties  of  the  ionic  liquid  and  that  the  nature  of 
the  cation  (positive  ion)  affects  the  stability  and  reactivity  of  the  ionic  liquid.  Therefore,  it  is 
possible  to  adjust  IL  physico-chemical  properties  by  making  changes  in  their  component  ion 

i  o 

structures  to  suit  requirements  when  designing  energetic  materials.  Approximately  10 
different  ionic  liquids  are  possible.  The  number  of  known  ionic  liquids  has  expanded 
exponentially  over  the  past  two  decades  and  it  is  certain  that  new  classes  will  be  developed  in 
the  forthcoming  years.  Therefore,  it  is  not  practical  to  study  experimentally  all  of  these  ionic 
liquids.  Hence,  a  few  theoretical  tools  have  been  developed  for  some  of  the  properties  and 
new  ones  have  to  be  developed  to  predict  remaining  physico-chemical  properties  so  as  to 
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allow  for  an  efficient  approach  to  screen  and  design  energetic  ionic  liquids  [Catoire  L.  et  al., 
2008]  . 

Basic  scientific  data  concerning  ionic  liquids,  such  as  their  physico-chemical 
properties,  their  thermal  stability,  etc.,  are  reported  here.  The  ionic  liquids  considered  are  not 
all  energetic  ionic  liquids  (EILs),  but  trends  observed  with  non-energetic  ionic  liquids  may 
also  be  of  interest  to  EILs. 


1.1.2  Properties 

Smiglak  et  al.  have  reported  the  performance  requirements  for  EILs  (see  Table  1) 
[Smiglak  M.  et  al.,  2007]. 


Physical  properties 

Melting  Point 

Less  than  -  40°C 

Density 

>1.4  g.cm'3 

Viscosity 

As  low  as  possible 

Surface  tension 

<  1.10"3  N.cm'1 

Therma 

Stability 

TGA  (75°C,  isothermal) 

<1%  loss  of  material  over  24hr 

TGA  (rate  10°C/min) 

>  120°C  for  T decomp 

Thermodynamic  Properties 

Heat  of  Formation 

As  positive  as  possible 

Heat  of  Combustion 

>  6  kcal.g'1 

Toxicology 

LD5o 

>  0.5  g.kg'1 

AMES 

Negative 

Legend:  TGA:  Thermogravimetric  analysis;  LD50:  Median  Lethal  Dose;  AMES:  The  Ames  test  is  a 
widely  employed  method  that  uses  bacteria  to  test  whether  a  given  chemical  can  cause  cancer. 

Table  1:  Performance  requirements  for  EILs  [Smiglak  et  al.,  2007] 


Energetic  ionic  liquids  of  interest  are  salts  that  have  melting  points  that  fall  well  below 
room  temperature  (-40  °C).  Since  it  is  impossible  to  measure  the  melting  points  for  the  large 
number  of  ILs  and  EILs,  numerical  methods  have  to  be  derived  and  validated  (with  little  data 
available).  Several  attempts  have  been  made  [Catoire  et  al.,  2008]  but  all  of  them  suffer  from 
several  weaknesses.  The  most  unexpected  result  reported  by  Catoire  et  al.  concerns  the  fact 
that  for  some  of  the  ionic  liquids  studied,  the  melting  and  the  freezing  points  are  different. 
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This  is  due  to  the  fact  that  ionic  liquids  tend  to  super-cool  [Van  Valkenburg  et  ah,  2005]. 
Another  interesting  feature,  potentially  of  importance  for  propulsion  applications,  is  that  these 
authors  report  that  some  ionic  liquids  solidified  at  room  temperature  after  several  months. 

Another  interesting  feature  of  ionic  liquids  is  that  they  are  generally  denser  than 
organic  liquids.  For  propulsion  applications  this  is  of  great  interest.  Smiglak  et  al.  report  that  a 
density  above  1 .4  gem"  is  required  for  EILs.  One  other  important  fact  reported  is  that  ionic 
liquids  are  generally  not  pure  and  that  water  is  a  very  common  contaminant,  which  can 
decrease  the  density  of  the  ionic  liquid  [Van  Valkenburg  et  al.,  2005],  Fredlake  and  co¬ 
workers  have  measured  the  density  of  some  imidazolium-based  ionic  liquids  [Fredlake  et  al., 
2004].  A  decrease  in  the  density  has  been  observed  with  an  increase  in  temperature,  as 
expected  for  molecular  liquids.  Other  general  trends  for  these  imidazolium-based  ionic  liquids 
are:  the  density  falls  off  as  the  alkyl-substituted  chain  length  of  the  cation  increases  and  the 
density  rises  with  increasing  molecular  weight  of  the  anion.  Other  relevant  physio-chemical 
data  has  been  reported  by  Catoire  and  co-workers  [Catoire  et  al.,  2008,  Interim  report]. 


1.1.3  Interest  in  Ionic  Liquids 

Concerning  propulsion  applications,  EILs  have  several  advantages  compared  to 
existing  hypergolic  fuels  such  as  hydrazine  (N2H4)  and  its  methylated  derivatives  MMH 
(monomethylhydrazine,  CH3NHNH2)  and  UDMH  (unsymmetrical  dimethyl  hydrazine, 
(CH3)2NNH2).  These  ’traditional’  fuels  are  toxic  (short-term  effects)  and  carcinogenic  (long¬ 
term  effects)  at  very  short  exposures,  which  can  be  very  problematic  since  they  have 
relatively  high  vapor  pressures.  Moreover,  their  flashpoints  are  quite  low  and  this  can  be  a 
source  of  chemical  risk  (explosion)  during  handling  and  storage.  On  the  contrary,  ionic  liquids 
are  thermally  stable  and  have  a  negligible  vapor  pressure,  high  density  and  large  liquidous 
range.  Physico-chemical  properties  are  available  only  for  relatively  few  compounds,  even  for 
such  basic  properties  as  density  and  melting  point.  Performance  requirements  for  EILs  in 
terms  of  melting  point,  surface  tension,  density,  viscosity,  etc.,  are  given  in  Table  1,  [Smiglak 
et  al.,  2007]. 


1.1.4  Chemical  nature  of  vapor  above  ionic  liquids 

Different  experimental  studies  have  been  carried  out  in  the  last  few  years,  [Catoire  et 
al.,  2008  and  Catoire  et  al.,  2009],  to  determine  the  enthalpy  of  formation  of  ILs.  Enthalpy  of 
formation  is  not  only  required  in  the  liquid  phase  but  also  in  the  gas  phase  above  the  ionic 
liquid  since  ignition  occurs  in  the  gas  phase.  It  is  experimentally  far  from  being  obvious  to 
perform  constant-volume  combustion  calorimetry  with  the  gas  phase  above  an  ionic  liquid. 
First,  the  vapor  pressure  of  an  ionic  liquid  is  generally  negligible  at  room  temperature. 
Second,  once  the  ionic  liquid  is  heated,  the  vapor  pressure  above  the  liquid  increases  but  an 
analysis  has  to  be  performed  to  determine  whether  or  not  the  vapor  above  the  liquid  is 
representative  of  the  ionic  liquid.  In  other  words,  the  gas  phase  can  be  composed  of  gases 
resulting  from  the  thennal  decomposition  of  the  ionic  liquid.  For  computations,  this  is  not 
necessarily  a  problem  but  in  that  case  one  can  hardly  define  an  unambiguous  enthalpy  of 
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formation  of  the  vapor  above  the  liquid.  One  has  to  know  beforehand  the  chemical  nature  of 
the  gas  phase  above  the  ionic  liquid  [Catoire  et  ah,  2008  and  Catoire  et  ah,  2009]. 

The  chemical  nature  of  the  vapor  above  an  ionic  liquid  is  difficult  to  determine 
experimentally.  It  is,  however,  important  to  know  how  the  vapor  pressure  forms  above  an 
ionic  liquid.  Four  possibilities  exist.  At  equilibrium  above  the  liquid  phase,  ion  pairs  can  exist, 
or  separated  ions  can  exist,  or  neutral  pairs  can  exist,  or  finally  separated  neutrals  can  exist. 
Experimental  findings  and  theoretical  calculations  are  reported  in  the  literature  on  this  issue. 
It  appears  that  ionic  liquids  have  to  be  divided  into  two  classes:  pro  tic  ionic  liquids  and 
aprotic  ionic  liquids. 

The  data  reported  below  has  been  taken  from  our  previous  reports  [Catoire  et  ah,  2008 
and  Catoire  et  ah,  2009]. 


Protic  Ionic  Liquids 

According  to  Earle  and  co-workers,  for  protic  ionic  liquids,  equilibrium  exists  between 
ionic  and  non-ionic  species  [Earle  et  ah,  2006]: 

BH+  +  X  ({)  B(£)  +  HX(t)  <->  B(g)  +  HX(g) 

Where  (£)  represents  the  liquid  phase  and  (g)  represents  the  gas  phase.  Calculations  made  by 
Schmidt  and  co-workers  on  protic  triazolium-based  energetic  ionic  liquids  show  that  the 
neutral  pair  is  of  lower  energy  than  the  ionic  pair  [Schmidt  et  al,  2005].  Separation  of  the 
species  in  the  ion  pair  requires  a  significant  amount  of  energy;  >50  kJ  mol1.  One  of  the 
conclusions  of  this  paper  is  that  the  ionic  pair  once  formed  will  transfer  one  proton  to  form  a 
neutral  pair.  This  proton  transfer  is  therefore  fundamental  for  the  decomposition  of  the 
triazolium-based  species  considered.  However,  Baranyai  and  co-worker  did  not  observe  the 
formation  of  species  according  to  the  proton  transfer  process,  and  therefore,  their  mechanism 
is  questionable  [Baranyai  et  al.,  2004],  Zorn  et  al.  deduced  from  their  calculations  on  protic 
tetrazolium-based  ionic  liquids  that  when  a  cation  is  paired  with  an  oxygen  rich  anion,  an 
isolated  gas  phase  ion  pair  was  not  generally  found  to  be  stable  [Zorn  et  al,  2006].  In  fact,  a 
proton  transfers  barrierlessly  from  the  cation  to  the  anion  to  form  a  neutral  pair.  This  result  is 
consistent  with  the  experimental  findings  of  Leal  and  co-workers,  who  report  that  the  gas 
phase  above  a  protic  ionic  liquid  is  composed  of  separated  neutral  molecules  [Leal  et  al., 
2007]. 


Aprotic  Ionic  Liquids 

Calculations  made  by  Emel’yanenko  and  co-workers  on  the  gas  phase  equilibrium: 

ion  pair  <->  separated  ions 

for  l-butyl-3-methylimidazolium  dicyanamide  indicate  that  dissociation  of  the  ion  pair  is 
marginal  [Emel’yanenko  et  al.,  2007].  Experimental  and  theoretical  works  reported  by 
Strasser  and  co-workers  suggest  that  the  vapor  pressure  of  an  aprotic  ionic  liquid  (l-ethyl-3- 
methylimidazolium  bis(trifhroromethylsulfonyl)imide)  is  composed  of  isolated  ion  pairs 
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[Strasser.  et  al.,  2007].  The  presence  of  clusters  of  ion  pairs  seems  to  be  negligible.  This 
conclusion  is  supported  by  another  experimental  work  perfonned  by  Leal  et  al.,  who  find  that 
the  gas  phase  is  composed  of  isolated  ion  pairs  and  that  no  experimental  evidence  supports  the 
formation  of  ion  clusters  [Leal  et  al.,  2007].  For  aprotic  ionic  liquids,  the  only  possible 
mechanism  for  volatilization,  according  to  Earle  et  al.,  is  as  ion  pairs,  isolated  or  as  clusters 
[Earle  et  al.,  2006].  Armstrong  and  co-workers  perfonned  an  analysis  of  the  vapor  pressure  of 
1  -ethyl-3 -methylimidazolium  ethylsulfate  by  mass  spectrometry  (MS)  [Armstrong  et  al., 
2007].  They  deduced  from  their  experiments  that  the  vapor  pressure  of  these  ionic  liquids 
consists  of  ion  pairs,  although  they  do  not  observe  directly  these  ion  pairs  because  the  ionized 
ion  pair  (formed  in  the  ionization  chamber  of  the  mass  spectrometer  by  electron 
bombardment)  is  unstable  and,  therefore,  fragmentation  occurs  in  their  experiments.  One 
other  finding  by  Annstrong  and  co-workers  is  that  clusters  of  ion  pairs  do  not  seem  to  fonn  in 
the  gas  phase  [Armstrong  et  al.,  2007]. 


1.2  Hypergolic  mixtures  with  ionic  liquids 

The  aim  of  the  second  part  of  the  state  of  the  art  section  is  to  document  hypergolic  IL 
mixtures.  In  2011,  Schneider  and  co-workers  reported  hydrogen-rich  ionic  liquids  that  were 
hypergolic  with  hydrogen  peroxide  [Schneider  et  al.,  201 1].  Before  this  time,  no  ionic  liquids 
were  known  to  be  hypergolic  with  hydrogen  peroxide.  The  Schneider  et  al.  study  considered 
hydrogen-rich  ILs  with  metal  hydride  anions.  Two  different  ionic  liquids  were  tested  (see 
Figure  1). 

The  ionic  liquid,  listed  as  1,  is  the  trihexyltetradecylphosphonium  (THTDP)  cation 
with  borohydride  (BFLf)  anion.  The  second  ionic  liquid  listed  as  2  is  the  THTDP  cation  with 
aluminum  borohydride  ABH  (A1(BH4)4‘)  anion.  These  two  ionic  liquids  were  subjected  to 
drop  tests  to  detennine  their  reactivity  with  common  propulsion  oxidizers.  A  drop  test 
consists  of  dropping  the  ionic  liquid  into  the  oxidant  and  to  measure  the  time  between  the 
contact  of  the  two  liquids  and  the  time  when  ignition  is  observed.  The  (THTDP)(BH4~)  ignites 
only  3s  after  dropping  onto  H2O2,  while  the  ignition  delay  of  (THTDP)(A1(BH4)4")  was  quite 
short  (less  than  30  ms)  for  two  different  H2O2  solutions  (90  and  98%).  This  ionic  liquid  also 
reacted  with  N204  (NTO)  and  WFNA  (white  fuming  nitric  acid).  Schneider  et  al.  concluded 
that  a  RTIL  with  a  complex  ABH  anion  would  be  hypergolic  with  all  common  propellant 
oxidizers  including  H2O2  [Schneider  et  al.,  2011]. 
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Legend:  1:  Cation:  Trihexyltetradecylphosphonium  -  Anion:  Borohydride  ;  2:  Cation: 

Trihexyltetradecylphosphonium  -  Anion:  Tetra(tetrahydroborato)aluminide 

Figure  1:  Hydrogen-rich  ILs  with  metal  hydride  anions  [Schneider  el  al.,  2011] 


In  201 1,  Zhang  and  co-workers  studied  dicyanoborate-based  ionic  liquids  [Zhang  et  al. 
2011a].  Compared  to  borohydride-based  (BFLf)  and  cyanoborate-based  (BH3(CN)')  ionic 
liquids,  these  ionic  liquids  are  water-stable  and  can  be  synthesized  in  aqueous  media.  The  ten 
dicyanoborate-based  ionic  liquids  considered  are  shown  in  Figure  2,  consisting  of  five 
different  cations  (N,Ndimethylhydrazinium,  N-methylimidazolium,  pyridinium,  N- 
methylpyrrolidinium,  and  N-methyltriazolium)  with  two  different  substituents  (butyl  and  allyl 
groups)  [Zhang  et  al.  2011a].  All  of  these  BH2(CN)2"-based  ionic  liquids  were  found  to  be 
hypergolic  with  WFNA.  The  characteristics  of  these  ILs  are  given  in  the  Table  2. 
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Legend:  3:  1  -methyl-3 -butylimidazolium  ;  4:  l-butyl-N,N-dimethylhydrazinium  ;  5:  l-allyl-N,N- 
dimethylhydrazinium,  6:  1 -methyl-3 -ally limidazolium  ;  7:  1 -butylpyridinium  ;  8:  lallylpyridinium  ;  9:  1 -butyl  - 
N-methylpyrrolidinium  ;  10:  1-allyl-N-methylpyrrolidinium  ;  11:  l-methyl-4-butyl-l,2,4-triazolium  ;  12:  1- 
methyl-4-allyl- 1 ,2,4-triazolium. 

Figure  2:  Different  cations  tested  for  the  dicyanoborate-based  ionic  liquids  studied  by  Zhang 

et  al.  [Zhang  et  al.,  2011a] 


In  Table  2,  it  is  seen  that  these  ionic  liquids  have  a  long  liquidous  range,  between  - 
80°C  and  >200°C  for  many  of  them.  Compared  to  the  requirement  of  density  >1.4  g.cm"  for 
energetic  ionic  liquids  ([Smiglak  et  al.,  2007]),  the  densities  of  the  dicyanoborate-based 
materials  are  lower  and  vary  from  0.91  to  1.03  g.cm"  .  The  ignition  delay  times  of  these  ionic 
liquids  are  relatively  short  and  vary  over  a  narrow  range  from  4  to  32  ms.  When  the 
dicyanoborate  anion  of  the  ILs  is  replaced  by  the  dicyanamide  (N(CN)2_)  anion  and  by  the 
nitrocyanamide  (N(N02)(CN)‘)  anion,  the  ignition  delay  (measured  in  a  drop  test)  of  the 
mixture  increase  substantially.  For  example,  for  the  IL  with  the  butyl-dimethylhydrazinium 
cation  (4),  and  with  the  dicyanoborate  anion,  the  ignition  delay  is  equal  to  6  ms,  but  when  the 
anion  is  the  dicyanamide,  the  ignition  delay  is  equal  to  46,  and  228  ms  when  the  anion  is  the 
nitrocyanamide  (see  Table  2)  [Zhang  et  al.  2011a], 

With  good  physico-chemical  properties  (such  as  melting  point,  viscosity,  etc.)  and 
ignition  delay  time  as  short  as  4  ms,  dicyanoborate-based  ILs  appear  to  be  promising 
substitutes  for  the  replacement  of  hydrazine  and  its  derivatives  in  hypergolic  propellant 
systems  [Zhang  et  al.  2011a], 
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85 

207 

42 

43 

6  (NCA) 

91 

220 

1.11 

44 

46 

7<b| 

<  - 

80 

252 

0.96 

19.8 

18 

gib] 

<  - 

80 

203 

1.00 

13.5 

6 

g|b] 

<  - 

80 

303 

0.92 

22.3 

26 

9  (DCA) 

44 

10M 

<  - 

80 

259 

0.94 

16.6 

8 

n« 

<  - 

80 

220 

0.99 

29.9 

32 

12w 

<  - 

80 

217 

1.03 

21.0 

6 

[a]  Tm/Tt :  phase-transition  temperature,  Td:  decomposition  temperature 
(onset),  p:  Density  (25°C),  rj:  viscosity  (25  °C),  ID:  ignition  delay  time 
(WFNA),  DCA:  corresponding  dicyanamide-based  ionic  liquid,  NCA: 
corresponding  nitrocyanamide-based  ionic  liquid,  [b]  ~BH2(CN)2,  this 
work. 

Legend:  The  molecule  number  corresponds  to  the  one  given  in  the  Figure  2 

Table  2  :  Properties  of  the  dicyanoborate-based  ILs  and  the  corresponding  dicyanamide-  and 
nitrocyanamide-based  ionic  liquids  [Zhang  et  al.,  2011a] 

In  2009,  Gao  and  co-workers  studied  hypergolic  mixtures  with  2,2-dialkyltriazanium- 
based  ionic  liquids  [Gao  et  al.,  2009],  which  due  to  the  multiple  nitrogen-nitrogen  bond 
character  of  the  cation,  also  had  a  potential  for  having  a  high  heat  of  formation  (see  Figure  3). 
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Figure  3:  Synthesis  of  2,2-dimethyltriazanium  salts  [Gao  et  al.,  2009] 


The  physico-chemical  properties  of  the  ionic  liquids  with  a  2,2-dimethyltriazanium 
cation  are  given  in  the  Table  3.  The  results  show  that  the  ionic  liquids  #  2,  3  and  7  are 
hypergolic  with  NTO  (N2O4)  or  white  fuming  nitric  acid  (WFNA  =  100%  HNO3),  but  only 
the  ionic  salts  #  2  and  7  exhibit  very  short  ignition  delay  times  (less  than  10  ms).  These 
ignition  delay  times  were  measured  via  droplet  tests.  However,  the  ionic  salt  #  7  is  a  solid  at 
ambient  temperature  (Melting  Point  =  99.0  °C).  The  ionic  salts  #  1  and  7  have  densities  equal 
to  1.47  g.cm'  ,  and  meet  Smiglak  and  co-workers’  requirements  of  density  >1.4  g.cm' 
[Smiglak  et  al.,  2007].  ILs  #  2  and  3  have  densities  of,  respectively,  1.26  and  1.15  g.cm". 
Ionic  salts  #  4,  5,  and  6  do  not  exhibit  hypergolicity  with  N2O4  or  WFNA.  This  phenomenon 
was  not  explained  by  the  authors  (see  Table  3)  [Gao  et  al.,  2009]. 


Salts 

d'*l 

7  W 

•  m 

Td« 

Ignition  delayld| 
N204  WFNA 

1 

1.47 

— 

26 

nh 

2 

1.26 

-0.19 

145.7 

8 

16 

3 

1.15 

10.7° 

134.2 

sh 

22 

4 

1.48 

- 

134.2 

nh 

nh 

5 

1.35 

47.8 

142.5 

nh 

nh 

6 

1.20 

74.4 

153.3 

nh 

nh 

7 

1.47 

99.0 

145.6 

10 

4 
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Legend:  [a]  density  (g.cm"3)  ;  [b]  melting  point  (°C)  ;  [c]  Decomposition  temperature  (°C) ;  [d]  in  ms  ;  nh  =  not 
hypergolic  ;  sh  =  hypergolic  when  a  second  drop  of  fuel  was  dropped  into  N204.  The  molecule  numbers 

correspond  to  the  ones  given  in  the  Figure  3 

Table  3:  Properties  of  energetic  2,2-dimethyltriazanium  salts  [Gao  et  al.,  2009] 


In  2010,  Zhang  and  co-workers  reported  the  synthesis  of  N,N-dimethylhydrazinium- 
based  ionic  liquids  [Zhang  et  al.,  2010].  The  salts  were  obtained  by  quaternizing  N,N- 
dimethylhydrazine  with  alkyl  halides  followed  by  metathesis  reactions  with  silver 
dicyanamide  or  silver  nitrocyanamide  (see  Figure  4). 
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12:  R  =  CH2CCH 
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15:  R  =  CH3 
16:  R  =  C2H5 
17:  R  =  C4H9 
18:  R  =  H2CHCH2 
19:  R  =  CH2CCH 
20:  R  =  CH2CH2OH 
21:  R  =  CH,CN 


Figure  4:  Synthesis  of  N,N-dimethylhydrazinium  ionic  liquids  [Zhang  et  al,  2010] 


The  physico-chemical  properties  of  these  salts  are  given  in  Table  4.  The  N,N- 
dimethylhydrazinium  salts  are  liquids  at  ambient  temperature  except  for  salts  #  8,  9,  15  and 
16.  There  is  no  detectable  phase-transition  temperature  (Tm  and  Tg)  for  compound  #  11-14 
and  18-21  in  the  temperature  range  between  -25  and  60  °C.  The  decomposition  temperatures 
of  these  ILs  vary  from  144.8  (compound  #  14)  to  296.9  °C  (compound  #  16),  with  the 
nitrocyanamide  anion  ILs  exhibiting  greater  thermal  stability  than  the  dicyanamide  analogues 
[Zhang  et  al.,  2010].  Density  and  viscosity  of  the  N,N-dimethylhydrazinium  ionic  liquids 
vary  from  1.01  to  1.28  g.cm"  and  67.5  to  1310  mPa.s,  respectively.  The  nitrocyanamide  salts 
have  a  higher  density  and  viscosity  than  the  corresponding  dicyanamides  due  to  the  presence 
of  the  nitro  (NO2)  group  (see  Table  4). 

All  of  these  ionic  liquids  show  hypergolicity  with  WFNA.  Among  these,  salts  #9,  11, 
12  and  13  have  ignition  delay  times  below  or  equal  to  40  ms.  Presence  of  the  nitrocyanamide 
anion  instead  of  the  dicyanamide  anion  in  the  IL  results  in  an  increase  in  the  ignition  delay 
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time  for  a  given  cation.  For  nitrocyanamide  salts,  the  shortest  ignition  delay  (ID)  time  is 
observed  for  compound  #  15  with  an  ID  equal  to  126  ms  (see  Table  4)  [Zhang  et  ah,  2010]. 


Salts 

7m(78)l>l 

[°C] 

7/1 

f°C] 

[gem  3] 

,;tdl 

[mPas] 

ID'e> 

[ms] 

8 

60.5 

253.5 

1.10 

— 

58 

9 

30.9 

267.1 

1.06 

67.5 

22 

10 

20.4 

263.3 

1.01 

113.9 

46 

11 

- 

199.2 

1.05 

78.6 

24 

12 

- 

174.3 

1.13 

228.6 

30 

13 

- 

236.0 

1.15 

161.8 

40 

14 

- 

144.8 

1.17 

1057.0 

1286 

15 

35.2 

292.4 

1.24 

- 

126 

16 

25.4 

296.9 

1.17 

- 

198 

17 

9.0 

285.5 

1.11 

119.5 

228 

18 

- 

208.2 

1.16 

84.9 

130 

19 

- 

189.3 

1.21 

269.8 

134 

20 

- 

269.1 

1.26 

185.9 

247 

21 

- 

193.5 

1.28 

1310.0 

1642 

Legend:  [a]  Melting  point ;  [b]  Decomposition  temperature  ;  [c]  Density  ;  [d]  Viscosity  (25°)  ;  [e]  Ignition  delay 
time  with  WFNA.  The  molecule  number  corresponds  to  the  one  given  in  the  Figure  4. 

Table  4:  Properties  of  N,N-dimethylhydrazinium  ionic  liquids  [Zhang  et  al.,  2010] 


In  2012,  Wang  and  co-workers  investigated  hypergolicity  of  boronium-cation  based 
ionic  liquids  [Wang  et  al.,  2012].  With  additional  literature  data,  the  authors  concluded  that 
molecules  with  strongly  reduced  B-H  bonds  appear  to  enhance  hypergolic  ignition.  The 
synthesis  scheme  for  boronium-cation  based  ILs  used  by  Wang  et  al.  is  shown  in  Figure  5. 
Due  to  the  short  B-H  bond,  a  strong  interaction  between  B  and  H  exists,  which  leads  to  a 
weak  reducing  potential  of  hydrogen.  This  may  explain  the  excellent  stability  of  these  ionic 
liquids  in  water  [Wang  et  al.,  2012]. 
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Figure  5:  Synthesis  of  boronium-based  ionic  liquids  [Wang  et  al.,  2012] 


The  physico-chemical  properties  of  these  boronium-based  ionic  liquids  are  shown  in 
Table  5.  All  of  these  eight  compounds  have  melting  points  below  100  °C,  but  only  5  of  them 
are  liquids  at  room  temperature  (2b,  3b,  4a,  4b  and  5b).  The  decomposition  temperatures  of 
these  ionic  liquids  range  from  190  to  306  °C,  showing  that  they  are  thermally  stable  [Wang  et 
al.,  2012]. 

High  density  and  low  viscosity  are  two  important  properties  of  hypergolic  ionic 
liquids.  However,  high  (or  low)  density  is  often  concomitant  with  high  (or  low)  viscosity.  For 
example,  IL  2b  has  the  second  highest  density  (1.15  g.cnf  )  and  the  highest  viscosity  (844 

•5 

mPa.s),  on  the  contrary,  IL  5b  possesses  the  lowest  density  (1.05  g.cnf  )  and  viscosity  (35 
mPa.s).  Ignition  delays  of  these  ILs  with  WFNA  were  measured  via  droplet  tests.  A 
characteristic  light-green  flame  due  to  boron  combustion  was  observed  upon  ignition.  Ionic 
liquid  5b  has  the  shortest  ignition  delay  (14  ms)  followed  by  ILs  3b  and  4b  (both  18  ms)  and 
3a  (26  ms).  If  we  compare,  the  methyl-substituted  series-a  and  the  allyl-substituted  series-b, 
the  first  one  gives  longer  ignition  delay  times,  suggesting  the  positive  impact  of  the  allyl 
substituent  on  the  ignition  properties  of  hypergolic  mixtures  with  WFNA.  The  authors 
concluded  that  ionic  liquid  5b  is  the  most  promising  hypergolic  fuel  (short  ID,  low  viscosity 
and  liquid  over  a  wide  range  of  temperature)  (see  Table  5)  [Wang  et  al.,  2012]. 
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Salt 

T  M 
*  m 

[*q 

7/' 

ra 

p" 

[gem  3] 

ID  [d] 

[ms-1] 

[s] 

//  [0 
[mPas] 
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1.28 

nh 

183.1 
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2  b 

<-80 

198 

1.15 

64 

178.2 

844 

3  a 

30 
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1.12 

26 

160.4 

43 

3  b 

<-80 

190 

1.09 

18 

159.0 

69 

4a 

<-80 

227 

1.20 

34 

181.3 

60 

4h 

<-80 

236 

1.14 

18 

176.5 

103 

5a 

62 

306 

1.12 

h“ 

164.5 

- 

5  b 

<-80 

266 

1.05 

14 

162.4 

35 

Legend:  [a]  Melting  point ;  [b]  Decomposition  temperature  ;  [c]  Density  ;  [d]  Ignition  delay  time  with  WFNA  ; 
[e]  Specific  impulse  ;  [f]  Viscosity  (25°) ;  nh  =  not  hypergolic  ;  h  =  hypergolic.  The  molecule  numbers 

correspond  to  the  ones  given  in  the  Figure  5. 

Table  5:  Properties  of  boronium  cation-based  ionic  liquids  [Wang  et  al.,  2012] 


In  2010,  He  and  co-workers  studied  hypergolicity  of  nitrocyanamide  anion-based  ionic 
liquids  [He  et  al.,  2010].  Synthesis  of  imidazolium-based  salts  with  the  nitrocyanamide  anion 
is  shown  in  Figure  6.  These  ILs  were  synthesized  by  anion  exchange  between  imidazolium- 
based  halides  and  silver  nitrocyanamide  in  methanol.  After  filtration,  the  filtrates  were  dried 
to  give  light-yellow  liquids  (compounds  1  to  4)  or  a  light-brown  solid  (compound  5). 


Imidazolium-based  salts 


X  =CI,  Br  or  I 

1:  R  =  ethyl;  2:  R  =  n-butyl;  3:  R  =  allyl; 

4:  R  =  2-methoxylethyl;  5:  R  =  cyanomethyl 

Legend:  NCA"  =  nitrocyanamide  anion  N(CN)(N02)" 

Figure  6:  Synthesis  of  nitrocyanamide  anion-based  ionic  liquids  [Fie  et  al.,  2010] 


The  physico-chemical  properties  of  these  compounds  are  summarized  in  Table  6.  All 
of  these  have  melting  points  below  100  °C  but  only  ILs  1  to  4  are  liquids  at  room  temperature. 
Compound  5  possesses  a  melting  point  of  83  °C,  which  is  due  to  the  cyano-methyl  group 
interactions.  All  of  these  ionic  liquids  have  high  decomposition  temperatures  (Td  above  200 
°C)  and  low  viscosities  (q  below  60  cP),  which  are  the  preferred  attributes  for  any  new 
bipropellant  application.  Compounds  #  1,  2,  3  and  4  were  screened  for  hypergolic  activity  in 
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drop  tests.  All  were  found  to  be  hypergolic  with  WFNA  with  values  of  ignition  delay  times 
ranging  between  46  (compound  3)  to  81  ms  (compound  2)  (see  Table  6)  [He  et  ah,  2010]. 


Compound 

Tm/Tg|a|  (°C) 

Td|b|  (°C) 

p|c|  (g.cnT3) 

Vd|  (cP) 

ID|e|  (ms) 

1 

-1-13 

253 

1.18 

23 

78 

2 

-1-90 

256 

1.13 

57 

81 

3 

-/- 91 

220 

1.11 

44 

46 

4 

-/-83 

266 

1.21 

54 

65 

5 

83/- 

203 

1.47 

Legend:  [a]:  Melting  point/phase  transition  temperature  ;  [b] :  Decomposition  temperature  ;  [c] :  Density, 
25°C  ;  [d] :  Viscosity,  25°C  ;  [e] :  Ignition  delay  time  with  WFNA  (100%  HN03).  The  molecule  numbers 

correspond  to  the  ones  given  in  the  Figure  6. 

Table  6:  Physicochemical  properties  of  nitrocyanamide  anion-based  ionic  liquids  [He  et  al., 

2010] 


In  2008,  Schneider  and  co-workers  studied  the  hypergolicity  of  different  dicyanamide 
anion-based  ionic  liquids  [Schneider  et  al.,  2008].  The  authors  chose  these  dicyanamide- 
based  ionic  liquids  because  these  ionic  liquids  possess  some  of  the  lowest  viscosities  among 
known  ionic  liquids.  The  imidazolium  cation  was  considered  for  this  study  because  they 
generally  possess  greater  stability  than  their  triazolium  or  tetrazolium  analogues.  The 
synthesis  of  these  ILs  is  shown  in  Figure  7. 


\©y 


N  RrO 


AgN(CN)2 


Br  CHjOH,  -AgBr 


N  ©  N 
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CN 


CN 

R  =  -allyl(l),  -(3-butenyl)(2),  -propargyl(3),  -(2-butynyl)(4),  -(2-pentynyl)(5) 


Cation  6  :  l-methyl-4-amino-l,2,4-triazolium  (the  formation  pathway  is  identical  to  the  others) 

Figure  7:  Synthesis  of  dicyanamide  anion-based  ionic  liquids  [Schneider  et  al.,  2008] 


Table  7  shows  the  physico-chemical  properties  of  these  dicyanamide-anion  based 
systems.  All  of  these  compounds  are  ionic  liquids  (melting  points  below  100  °C),  but 
compounds  #  4  and  5  are  solids  at  room  temperature.  Only  ILs  1  and  2  show  good  thermal 
stability  (Td  above  200  °C)  and  low  viscosity  with  values  equal  to  42  and  27  cP,  respectively. 
The  ignition  delay  times  of  these  ILs  were  measured  via  droplet  tests  with  WFNA. 
Compounds  #  1,  3  and  6  were  found  to  be  hypergolic  with  the  WFNA.  Compounds  4  and  5 
were  not  tested  because  of  their  unsuitability  as  liquid  propellants  (see  Table  7)  [Schneider  et 
al.,  2008], 


Compound 

Tg|a|  (°C) 

Tm|b|  (°C) 

Td|c|  (°C) 

q  |d|  (cP) 

ID  |e|  (ms) 
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1 

-85 

207 

42 

43 

2 

-90 

210 

27 

Not 

measured 

3 

-61 

17 

144 

110 

15 

4 

49 

179 

5 

59 

184 

6 

-66 

143 

92 

31 

Legend:  [a]  Glass  transition  temperature  [b] :  Melting  point ;  [c] :  Decomposition  temperature  ;  [d]:  Viscosity, 
25°C  ;  [e]:  Ignition  delay  time  with  WFNA  (100%  HNO3).  The  molecule  numbers  correspond  to  the  ones  given 

in  the  Figure  7. 

Table  7:  Properties  of  dicyanamide -anion  based  ionic  liquids  [Schneider  el  al.,  2008] 


In  201 1,  Zhang  and  co-workers  reviewed  all  of  the  recent  work  on  hypergolic  mixtures 
[Zhang  et  al,  2011b],  Also  in  this  review,  they  extended  the  work  of  Joo  and  co-workers  on 
azide-containing  hypergolic  mixtures  [Joo  et  al.,  2010], 

The  azido  group  is  a  highly  energetic  moiety,  which  can  add  approximately  280 
kJ.mol1  to  the  energy  content  of  a  molecule.  2-azido-N,N-dimethylethylamine  (DMAZ)  is  a 
reduced-hazard  liquid  fuel  with  many  of  its  physical  properties  comparable  to  those  of 
monomethylhydrazine  (MMH)  fuel.  Literature  data  is  available  suggesting  DMAZ  is  less 
toxic  than  MMH.  Methylation  of  DMAZ  with  iodomethane,  with  a  slight  excess  of  the 
corresponding  silver  salts  (silver  nitrocyanamide  and  silver  dicyanamide)  results  in  ILs  23  to 
26.  When  hydrazoic  acid  (HN3)  was  reacted  with  N,Ndimethylisopropylamine  or 
triethylamine,  the  corresponding,  white-solid,  ammonium  azides  (27  &  28)  were  obtained. 
The  latter  two  compounds  were  shown  to  be  first  examples  of  hypergolic  salts  containing  an 
azide  anion  (see  Figure  8)  [Zhang  et  al,  2011b]. 


(a) 


R  =  CH3;  Y  =  N(N02)(CN)  (23),  N(CN)2  (24) 

R  =  CH2CH2N3;  Y  =  N(N02)(CN)  (25),  N(CN)2  (26) 


(27) 


(b) 


hn3 

- ► 


(28) 
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Figure  8:  Synthesis  of  hypergolic  ionic  liquids  (a)  with  mono-  or  bis-  (azidoethyl)ammonium 
cations  and  (b)  with  azide  anion  [Zhang  et  al.,  2011b] 


Table  8  summarizes  the  physico-chemical  properties  of  these  azide  cation  (and  anion) 
-based  salts.  Compounds  23,  24,  27  and  28  are  ionic  liquids  (Tm  below  100  °C)  but  only  one 
of  them  (#  24)  is  a  liquid  at  ambient  temperature.  Compounds  23  to  26  have  relatively  high 
decomposition  temperatures  (Td  above  200°C)  and  high  densities.  On  the  contrary,  ILs  27  and 

•3 

28  have  low  densities  with  values  around  1  g.cm'  .  Hypergolicity  of  these  ILs  was  evaluated 
via  droplet  testing  with  WFNA  (23  to  26)  or  N2O4  (27  and  28).  Compounds  27  and  28  are 
protic  liquids  and  were  found  to  be  hypergolic  only  with  N2O4  (see  Table  8)  [Joo  et  al.,  2010 
and  Zhang  et  al,  2011b]. 


Compd 

7  hi 

1  m 

[°C] 

T/l 

[°q 

dM 

[gem"3] 

ID[el 

[ms] 

A  H/n 
[kj  g"1] 

/  Is) 
sp 

[s] 

23 

28 

245 

1.24 

8 

1.9 

218 

24 

9 

235 

1.15 

20 

2.4 

202 

25 

- 

222 

1.32 

226 

3.2 

231 

26 

- 

222 

1.21 

16 

3.8 

221 

27 

75 

subl. 

0.99 

hgw 

(n2o4) 

3.3 

249 

28 

80 

subl. 

1.01 

hgM 

(n2o4) 

3.2 

245 

[a]  Melting  point,  DSC,  10°Cmin-1.  [b]  Decomposition  onset,  DSC, 
10°Cmin_1;  subl.  =  sublimes.  [c]  Density,  25°C.  [d]  Viscosity,  25 °C. 

[e]  Ignition  delay  (WFNA).  [f]  Enthalpy  of  formation,  [g]  Specific  impulse, 
calculated  isobarically  at  68  atm  (Cheetah  5.0).  [h]  Hypergolic. 

Legend:  The  molecule  numbers  correspond  to  the  ones  given  in  the  Figure  8. 

Table  8  :  Properties  of  azide-functionalized  ionic  liquids  [Zhang  et  al.,  2011b] 


Table  9  shows  all  of  the  systems  for  which  hypergolicity  has  been  observed 
experimentally.  It  is,  however,  worthy  to  note  that  no  composition  has  been  reported  either  in 
the  liquid  phase  or  in  the  gas  phase  for  these  cases. 


Hypergolic  mixtures 

Cation 

Anion 

Oxidant 

NH9 

1  +  2 

l-LC-N  — NH9 

3  !  2 

ch3 

N(CN)(N02)- 

WFNA  and  N204 

NCV 

WFNA 

N(CN)2- 

WFNA  and  N204 

cr 

N204 
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NH9 

1  +  2 

H3C-N-CH3 

ch3 

N(CN)(N02)' 

WFNA 

N(CN)2- 

WFNA 

Vn2 

h3c-n— ch3 

h3c-h2 

BH2(CN)2- 

WFNA 

N(CN)(N02)' 

WFNA 

N(CN)2- 

WFNA 

NH, 

1  +  2 

H3C-N-CH3 

hc^h> 

II 

ch2 

BH2(CN)2- 

WFNA 

N(CN)(N02)' 

WFNA 

N(CN)2- 

WFNA 

NH, 

1  + 

H3C-N-CH3 

If™2 

CH 

BH2(CN)2- 

WFNA 

N(CN)(N02)' 

WFNA 

N(CN)2- 

WFNA 

NH, 

H,C  — N  — CH, 

O  |  O 

h2c-ch2 

2I 

OH 

N(CN)(N02)' 

WFNA 

N(CN)2- 

WFNA 

H 

h5c2—  n-c2h5 
c2h5 

N3' 

N204 

H  CH3 

1  +  /  3 

H,C— N— CH 

3  1  \ 
ch3  ch3 

N3‘ 

n2o4 

CH, 

1  + 

H3C-N-C-C-N=N=N 

1  H,H, 
ch3  2  2 

N(CN)(N02)' 

WFNA 

N(CN)2- 

WFNA 

CH, 

1  + 

H,C-N  —  C  —  C  —  N  =  N  =  N 

1  H,  H, 

N(CN)(N02)' 

WFNA 

N(CN)2- 

WFNA 

I 

N 

I  I 
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H 

H  C  C 

1  1  ♦ 

H  C  ^  C  H 

3  H  1  2 

H,C  . 

2  ch2 

CH3 

BH2(CN)2‘ 

WFNA 

N(CN)(N02)- 

WFNA 

N(CN)2- 

WFNA 

H 

HC  C 

1  l  + 

H,C  CH? 

3  HI2 

HC^ 

"ch2 

BH2(CN)2- 

WFNA 

N(CN)(N02)‘ 

WFNA 

N(CN)2- 

WFNA 

1 

H 

H  C  C 

1  1  ♦ 

H  ,C  C  ^  CH, 

3  H  1  2 

^CH 

N(CN)2- 

WFNA 

H 

HC  C 

1  1  + 

>L 

H„C  CH„ 

3  HI2 

ch3 

N(CN)(N02)‘ 

WFNA 

H 

H  C  C 

1  1  ♦ 

H  ,C  CH, 

3  H  1  2 

HC  ^ 

2  0 

1 

ch3 

N(CN)(N02)- 

WFNA 

//  %  +  H2 

\  N-C  H 

N — 2  c— C 

H  N 

n2  CH3 

BH2(CN)2- 

WFNA 

_ 1 

//  ^  +  H2 
<7  N-C  h 

N - 2  c— C 

n2  CH3 

BH2(CN)2- 

WFNA 

Hj, 

H  2  C 

H,C  /  '"'CH, 

3  \  +'  ^C  3 

£j  "2 

BH2(CN)2- 

WFNA 

N(CN)2- 

WFNA 
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m.c.  .A,/"- 

77  H 

BH2(CN)2- 

WFNA 

H 

N  C 

1  1  + 

MM  ^ 

H  ,C  C  ^  CH, 

3  H  1  2 

H  C 

2  ch2 
ch3 

BH2(CN)2- 

WFNA 

H 

N  C 

1  1  + 

H,C  CH  „ 

3  HI2 

HC^ 

^ch2 

BH2(CN)2- 

WFNA 

Table  9:  Summary  of  all  the  systems  found  as  hyper golic  and  presented  in  this  report. 


1.3  Flame  Temperature  Criterion 

In  this  part,  the  flame  temperature  criterion  considered  for  the  construction  of 
flammability  limits  diagram  is  explained.  The  data  reported  here  are  taken  from  a  research 
proposal  by  Cato  ire  [Cato  ire,  2012], 

Previous  studies  on  non-hypergolic  systems  have  shown  that  the  adiabatic  flame 
temperature  at  the  LFL  and  UFL  remains  approximately  constant,  and  that  the  reported 
experimental  flammability  limits  correspond  to  adiabatic  flame  temperatures  in  the  1000-1500 
K  range  for  hydrocarbon/air  systems.  One  may  therefore  assume  that  flammability  diagrams 
for  hypergolic  mixtures  can  also  be  rationalized  with  such  a  simple  thermochemical  criterion. 
This  assumption  has  been  applied  to  the  monomethylhydrazine  (MMH)  -  nitrogen  tetroxide 
(NTO)  system  by  using  available  gas  phase  experimental  flammability  diagrams  at  room 
temperature  in  which  the  condition  of  rapid  adiabatic  mixing  was  assumed.  High  speed 
videography  was  used  to  record  the  auto-ignition  events.  An  example  of  igniting  and  non¬ 
igniting  zones  as  a  function  of  total  pressure  and  equivalence  ratio  at  a  given  mol%  of  He 
diluent  is  given  in  Figure  9.  Additional  data  is  available  in  the  research  proposal  [Catoire, 
2012]. 

As  a  proof-of-concept,  Catoire  and  co-workers  used  their  MMH/NTO  detailed 
chemical  kinetics  model  [Catoire  et  ah,  2004]  to  establish  a  numerical  criterion  to  predict 
flammability  diagrams  for  MMH/NTO  hypergolic  combinations.  Calculations  were  carried 
out  to  detennine  the  adiabatic  flame  temperature  at  constant  volume.  For  non-igniting 
mixtures,  the  calculated  flame  temperature  was  below  the  flame  temperature  calculated  for 
igniting  mixtures  (Figures  6  to  8  in  ref  Catoire,  2012).  For  a  stoichiometric  mixture  at  low 
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total  pressure,  in  which  ignition  was  observed,  the  temperature  profile  reported  in  Figure  10 
exhibits  a  flame  temperature  of  about  3000  K.  Moreover,  this  point  is  located  right  on  the 
ignition  limit  (pressure  of  6  kPa  and  equivalence  ratio  of  1.0).  This  observation  and  similar 
others  allow  for  the  establishment  of  a  simple  numerical  criterion  to  predict  ignition  events 
under  non-adiabatic  conditions,  which  are  more  representative  of  the  actual  experimental 
conditions  used  in  the  MMH/NTO/He  study;  i.e.,  a  rapid  one-step  or  two-step  temperature 
profile  leading  to  a  flame  temperature  at  least  equal  to  2900  K.  Near  the  limit,  in  the  non¬ 
ignition  zone  (just  below  the  lines  in  Figures  1-5  in  ref  Cato  ire,  2012),  flame  temperatures  of 
about  2700-2800  K  are  computed.  However,  it  is  to  be  noted  that  for  adiabatic  experimental 
conditions,  flames  are  likely  to  be  observed  in  this  zone  [Catoire,  2012]. 

The  temperature  criterion  detennined  here  is  dependent  on  the  experimental  setup 
used  (i.e.,  rapid  adiabatic  mixing),  and  it  is  assumed  that  the  threshold  for  the  theoretical 
flame  temperature  to  estimate  flammability  limits  is  much  higher  for  hypergolic  mixtures  than 
it  is  for  non-hypergolic  mixtures,  such  as  those  in  hydrocarbon-air  systems. 


Equivalence  ratio 


Figure  9:  Flammability  as  a  function  of  total  mixture  pressure  and  equivalence  ratio  obtained 
by  assuming  rapid  mixing  for  MMH/NTO/He  mixtures.  □  :  no  ignition;  w  :  ignition 
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Figure  10:  Typical  temperature  profile  computed  under  adiabatic  and  constant  volume 
assumptions  for  an  igniting  mixture  ( equivalence  ratio  1.0  and  low  total  pressure  of  6  kPa). 


Figure  1 1  shows  that  the  use  of  this  numerical  criterion  is  quite  consistent  with  the 
experimental  results.  This  agreement  between  experiments  and  computations  is  to  be 
expected  for  this  mixture  because  the  criterion  was  established  according  to  experimental 
results  of  Figure  9.  Furthermore,  the  good  agreement  also  observed  between  experiments  and 
computations  for  other  mixtures  (Figure  1 1  in  ref  Catoire,  2012)  shows  that  the  criterion  can 
be  extrapolated  to  other  MMH/NTO  mixtures  with  a  high  degree  of  confidence. 
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Equivalence  ratio 


Figure  11:  Computed  flammability  diagram  (mol%  He  =  53).  The  line  is  the 
experimental  limit  between  ignition  and  non-ignition  regions  given  in  Figure  9.  □  :  no 
ignition;  I :  ignition 


The  predictive  ability  of  this  thermodynamic  method  is  extended  successfully  to 
numerous  other  non-ionic  liquid/oxidant  mixtures  used  for  space  applications.  These  results 
have  been  published  in  Combustion  and  Flame  in  2017. 


To  summarize  this  part  of  the  report,  new  hypergolic  ionic  liquid/oxidant  mixtures  that 
have  recently  been  reported  in  the  literature  were  surveyed.  A  simple  thermodynamic 
criterion  for  the  adiabatic  flame  temperature  at  constant  volume  is  proposed  herein  to 
ascertain  if  a  given  IF/oxidant  mixture  will  be  hypergolic  or  not. 

In  the  present  work,  a  mixture  is  considered  hypergolic  if  its  adiabatic  flame 
temperature  at  constant  volume  is  equal  or  above  2700  K.  When  this  is  identified,  the 
flammability  limits  are  then  computed  as  described  below. 
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2  Thermochemistry  of  ionic  liquids 


Ionic  liquids  are  important  emerging  compounds  because  of  their  potential 
applications  as  energetic  propellants.  In  the  last  few  years,  efforts  have  been  made  to 
determine  their  thermodynamic  data,  but  data  are  still  scarce  even  for  the  simplest 
compounds. 

Thermodynamic  data,  such  as  standard  enthalpy  of  formation  at  298.15  K,  standard 
heat  capacities  (at  constant  pressure,  Cp°  or  at  constant  volume,  Cv°)  over  a  wide  range  of 
temperatures  and  standard  absolute  entropies  S°  over  a  wide  range  of  temperature  are  needed 
for  chemical  equilibrium  modeling  studies  [Catoire  et  ah,  2009]. 


2.1  Estimation  of  thermodynamic  data 


Ab  initio  calculations  have  been  shown  to  be  reliable  methods  for  estimating  the 
thermochemistry  of  many  compounds.  Three  methods  for  estimation  of  gas-phase  standard 
enthalpy  of  formation  at  298.15  K  have  been  proposed  by  Osmont  and  co-workers  [Osmont, 
2007  and  Osmont  et  ah,  2007]  by  using  quantum  chemistry  computations  at  the  B3LYP/6- 
31G(d,p)  level.  Two  of  these  methods  can  be  applied  to  ionic  liquids  .  One  of  these  methods 
is  devoted  to  general  organic  compounds  and  the  other  to  energetic  materials,  in  particular 
nitro  and  nitrate  compounds.  The  gas-phase  standard  enthalpy  of  formation  of  molecule  j  at 
298.15  K  can  be  determined  from  the  following  equation  [Catoire  et  ah,  2009  and  Osmont  et 
al„  2007]: 


627.51  x 


f 

Ej  +  ZPEj  +  thermal  corrections  +  ^ 


\ 

* 

CCi  C  • 


V 


J 


where  ai  is  the  number  of  atoms  i  in  molecule  j  and  C,*  the  atomic  correction  for  atom  i.  Ej 

and  ZPEj  denote,  respectively,  the  absolute  electronic  energy  and  zero-point  energy, 
calculated  using  the  GAUSSIAN  03  software  package.  The  units  are  Hartree  molecule'1  for 
Ej,  ZPEj  and  thermal  corrections,  and  Hartree  atom"1  for  Q*,  whereas  AfH°298.i5  K(g)  is  in  kcal 
mol'1.  Four  atomic  corrections  (for  C,  H,  O  and  N)  are  considered  here  for  the  method 
established  for  organic  compounds  and  seven  atomic  corrections  (for  C,  H,  O,  N,  Ca,  Oa,  and 
Na)  are  considered  for  the  method  established  for  energetic  compounds  with  one  more  atomic 
correction  for  hypervalent  nitrogen  (noted  Nhyp)  as  encountered  in  ionic  liquids.  The  values  of 
these  coefficients  are  given  for  the  energetic  compounds  and  the  organic  compounds, 
respectively  in  Table  10  and  Table  11. 
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Cp,  Cv  and  S  are  calculated  from  statistical  thermodynamic  equations  once  the 
vibration  frequencies  and  inertia  moments  have  been  determined  in  the  Gaussian  package 
[Catoire  et  ah,  2009]. 


Atom 

Q* 

H 

0.58188704 

C 

38.11447874 

N 

54.76547900 

O 

75.15682281 

Nd 

54.75925621 

Cd 

38.11542265 

Od 

75.15439445 

Nhyp 

54.75799267 

Table  10:  Atomic  corrections  used  in  the  energetic  compounds  method  [Osmont,  2007  and 

Osmont  et  al.,  2007] 


Atom 

G* 

H 

0.582238 

C 

38.114781 

N 

54.760090 

O 

75.150299 

Table  11:  Atomic  corrections  used  in  the  organic  compounds  method  [Osmont,  2007  and 

Osmont  et  al.,  2007] 


2.2  Validation  of  the  method 

The  enthalpy  of  formation  calculations  are  validated  by  comparison  with  literature 
values.  Table  12  compares  the  experimental  enthalpy  of  formation  values  with  those 
calculated  using  both  methods  for  different  imidazoles,  triazoles  and  tetrazoles.  These 
compounds  are  not  ionic,  but  they  are  structurally  quite  close  to  the  cations  present  in  the  ILs 
considered  here  [Catoire  et  al.,  2009]. 

Table  12  shows,  unambiguously,  that  the  method  established  for  energetic  compounds 
performs  better  than  the  method  established  for  general  organic  compounds.  The  former 
method  has  previously  been  established  for  nitro  and  nitrate  compounds  and  therefore  should 
provide  reasonably  accurate  results  for  the  ionic  liquids  [Catoire  et  al.,  2009]. 
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Compound 

Experimental 

AfH°298.15  K(g) 

(kcaLmof1) 

Calculated  AfH°298.i5  k 
(g)  with  the  method  for 
organic  compounds 
(kcaLmof1) 

Calculated  AfH°298.i5  k 
(g)  with  the  method  for 
energetic  compounds 
(kcaLmof1) 

1H- 1,2, 4-triazole 

46.1±0.2 

45.05 

40.7 

43.4 

2-ethyl- 1  H-imidazole 

16.3±0.31 

13.5 

15.6 

lH-tetrazole 

76.6±0.7 

71.7 

74.0 

1  H-imidazole 

33.29±0.45 

30.6±1.8 

29.1 

32.3 

2-methyl- 1  H-imidazole 

21.5±0.26 

18.5 

21.1 

1H- 1,2, 3 -triazole 

59 

57 

59.7 

lH-pyrazole 

42.9 

43.3 

39.2 

42.4 

Pyrrole 

25.9±0.1 

34.23 

26.4 

30.0 

Table  12:  Comparison  between  experimental  gas-phase  standard  enthalpy  of  formation  at 


298.15  K  and  calculated  gas-phase  standard  enthalpy  of  formation  at  298.15  Kfor 
imidazoles,  triazoles  and  tetrazoles  taken  from  [Ccitoire  el  al.,  2009]. 


However,  as  seen  previously  for  aprotic  ionic  liquids,  the  chemical  vapor  above  the  IL 
can  be  composed  of  ions.  It  is  therefore  important  to  also  validate  these  methods  for  cations 
and  anions.  Table  13  compares  the  experimental  enthalpy  of  formation  values  with  those 
calculated  with  both  methods  for  different  cations  and  anions.  The  experimental  enthalpies  of 
formation  are  taken  from  Mathieu  and  co-workers  [Mathieu  et  al.,  2010]  and  from  Halim  and 
co-workers  [Halim  et  al.,  1986]. 

The  results  show  that  for  the  cations,  the  method  for  organic  compounds  gives  better 
agreement.  This  result  is  excepted  because  the  cations  considered  are  not  energetic  ions. 
However,  for  the  anions,  both  methods  give  values  in  considerable  disagreement  with 
experimental  values.  Errors  of  up  to  145  %  (for  OH')  are  seen,  which  lead  us  to  conclude  that 
the  use  of  the  B3LYP/6-31G(d,p)  method  for  the  detennination  of  the  enthalpy  of  fonnation 
of  anions  is  not  reliable  (see  Table  13). 
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Compound 

Experimental 

AfH°298.i5K(g) 

(kcal.moT1)  [Mathieu 
et  ah,  2010  and  Halim 
etal.,  1986] 

Calculated 

AfH°298.i5K(g)  with  the 
method  for  organic 
compounds 
(kcal.moT1) 

Calculated 

AfH°298.l5  K  (g)  with 

the  method  for 
energetic  compounds 
(kcahmof1) 

NH4+(s) 

155 

153.9 

146.3 

CH(OH)2+(s) 

96 

96.1 

103.4 

CH3CO+(s) 

156 

155.3 

157.4 

N02'(s) 

-45.2 

-30.0 

-25.4 

N03-(s) 

-74.7 

-61.4 

-60.4 

CN'(s) 

17.7 

36.7 

36.6 

H3C-CH2-(s) 

35.1 

53.5 

52.1 

H3C-N‘-CH3(s) 

26.1 

35.2 

36.9 

OH'(s) 

-33.2 

11.0 

14.9 

C6H50'(s) 

(phenoxy) 

-40.5 

-28.5 

-23.1 

Table  13:  Comparison  between  experimental  gas-phase  standard  enthalpy  of  formation  at 
298.15  K  and  calculated  gas-phase  standard  enthalpy  of  formation  at  298.15  Kfor  selected 
cations  and  anions.  Calculation  method:  B3LYP/6-31G(d,p)  corrected  on  an  atomic  basis,  (s) 

stands  for  the  singlet  state. 


Another  ab  initio  method  (G2)  to  detennine  the  enthalpy  of  formation  for  anions  was 
tested.  This  method  is  more  precise  but  computationally  can  be  time  consuming.  The  enthalpy 
of  formation  is  given  by: 


A  fH°  298.15  K  (g)(IL)  = 


( 

( 

627.51  x 

G2  Enthalpy  (IL)  -  V  ai  x  G2Enthalpy(i ) 

V 

V  i 

JJ 

+  X«;XA  fH  °298.15K  (g)(i) 
i 


Here,  a;  is  the  number  of  atoms  i  in  the  IL  and  G2Enthalpy  denotes  the  enthalpy  at  298  K  of 
the  molecule,  calculated  using  the  GAUSSIAN  03  software  package.  The  units  are 
Hartree.molecule'1  for  G2Enthalpy,  whereas  AfH°298.i5K(g)  is  in  kcalmol  '. 


Table  14  compares  the  experimental  enthalpy  of  formation  values  and  the  calculated 
ones  for  different  anions.  The  results  show  good  agreement  between  the  calculated  and  the 
experimental  values. 
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Compounds 

G2Enthalpy 

AfH°298.15K(g)  Of 

atom  (kcal.mof 

*) 

Calculated 
AfH°298.15  K(g) 

(kcal.mol1) 

Experimental 
AfH°298.i5K(g) 
(kcal.mof1) 
[Mathieu  et  al., 
2010  and  Halim  et 
al.,  1986] 

C 

-37.785278 

171.29 

H 

-0.497639 

52.1 

0 

-74.980024 

59.43 

N 

-54.515599 

112.97 

B 

-24.599675 

133.84 

CN'(s) 

-92.725481 

17.8 

17.7 

N02'(s) 

-204.919145 

-46.4 

-45.2 

N03'(s) 

-280.040762 

-75.8 

-74.7 

OH'(s) 

-75.709467 

-33.9 

-33.2 

H3C-CH2(s) 

-78.95849 

38.4 

35.1 

H3C-N-CH3(s) 

-134.248824 

29.6 

26.1 

Table  14:  Comparison  between  experimental  gas-phase  standard  enthalpy  of  formation  at 
298.15  K  and  calculated  gas-phase  standard  enthalpy  of  formation  at  298.15  Kfor  different 

anions.  Calculation  method:  G2. 


2.3  Results 

In  this  part,  the  enthalpy  of  formation  and  the  thermodynamic  data  of  the  IL  studied 
are  presented.  The  enthalpy  of  formation  of  the  cations  are  presented  in  Table  15,  which  have 
been  determined  via  B3LYP/6-31G(d,p)  calculations  and  corrected  by  using  the  Osmont 
method  for  energetic  compounds.  The  anion  values  were  determined  via  G2  calculations. 

Due  to  excessive  computational  time  encountered  for  the  aprotic  ionic  liquids,  only 
thermodynamic  data  of  the  separated  ions  have  been  calculated. 

It  must  be  noted  that  we  were  not  able  to  determine  the  enthalpy  of  fonnation  for 
N(N02)(CN)"  due  to  convergence  problems  within  the  G2  method  as  implemented  in  the 
Gaussian  03  software  package. 
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Table  15:  Calculations  of  the  gas-phase  standard  enthalpy  of  formation  at  298.15  Kfor 
different  ions  of  interest  for  ionic  liquids  studied  in  this  report.  Calculation  methods: 
corrected  B3LYP/6-3 1  G(d,p)  for  the  cations  and  molecules;  and  G2  method  for  the  anions. 

Once  the  standard  enthalpy  of  formation  at  298.15  K  as  well  as  the  Cp,  Cv  and  S  had 
been  calculated,  the  14  coefficients  in  the  NASA  format  needed  for  thermodynamic 
calculations  were  determined  according  to  the  CHEMKIN-NASA  format. 

These  NASA  polynomials  have  the  form: 

Cp/R  =  al  +  a2  T  +  a3  T2  +  a4  T3  +  a5  T4 

H/RT  =  al  +  a2  T/2  +  a3  T2  /3  +  a4  T3/4  +  a5  T4/5  +  a6/T 

S/R  =  al  ln(T)  +  a2  T  +  a3  T2  /2  +  a4  T3  /3  +  a5  T4  /4  +  a7 
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Where  al,  a2,  a3,  a4,  a5,  a6,  and  a7  are  the  numerical  coefficients  supplied  in  the  NASA 
thermodynamic  files.  The  first  7  numbers  starting  on  the  second  line  of  each  species  entry 
(five  of  the  second  line  and  the  first  two  of  the  third  line)  are  the  seven  coefficients  (al 
through  al,  respectively)  for  the  high-temperature  range  (above  1500  K,  the  upper  boundary 
is  specified  on  the  first  line  of  the  species  entry).  The  remaining  seven  numbers  are  the 
coefficients  (al  through  al ,  respectively)  for  the  low-temperature  range  (below  1500  K,  the 
lower  boundary  is  specified  on  the  first  line  of  the  species  entry). 


The  thermodynamic  data  of  the  species  studied  here  are  given  below.  The  species  are 
named  MOL1,  MOL2,  MOL3,  etc.  The  correspondence  between  the  molecules  is  given  in 
Table  15.  Thennodynamic  data  for  N(N02)(CN)"  (Molecule  n°2)  have  not  been  determined. 


MOL1  H  10C  2N  3  G  298.150  5000.000  1500.00  1 

0. 1222135 1E+02  0.27003260E-01-0.99627348E-05  0.17079181E-08-0.1 1 138762E-12  2 
0.85934737E+05-0.421 10438E+02-0.56400683E+00  0.57280381E-01-0.36422690E-04  3 
0.1 1771380E-07-0.15058781E-1 1  0.90199776E+05  0.26186977E+02  4 

MOL3  H  11C  3N  2  G  298.150  5000.000  1500.000  1 

0.12506042E+02  0.29490288E-01-0.10890801E-04  0.18664465E-08-0.12162023E-12  2 
0.6954 1 1 82E+05-0.44885430E+02-0. 1 0450059E+0 1  0.58765690E-01-0.33305478E-04  3 
0.87793 139E-08-0.76556296E-12  0.74378328E+05  0.28557790E+02  4 

MOL4  C  2N  3  G  298.150  5000.000  1500.000  1 

8.39498290E+00  4.53143670E-03-1.81293070E-06  3.29552030E-10-2.24501790E-14  2 
1 .33068850E+04-1 .66833200E+0 1  3.69004260E+00  1.70717050E-02-1.43469600E-05  3 
5.89745910E-09-9.49972540E-13  1.47190690E+04  7.92532660E+00  4 

MOL5  H  17C  6N  2  G  298.150  5000.000  1500.000  1 

0.2076984 1E+02  0.46325712E-01-0.17189680E-04  0.29552207E-08-0.19297896E-12  2 
0.55083923E+05-0.88213717E+02-0.29713969E+01  0.97424363E-01-0.5607701  IE-04  3 
0.1481 1226E-07-0.12644377E-1 1  0.63580068E+05  0.40529523E+02  4 


MOL6  H  13C  5N  2  G  298.150  5000.000  1500.000  1 
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0.17281353E+02  0.36143 137E-01-0. 134373 19E-04  0.23154042E-08-0.15154042E-12  2 
0. 77 194076E+05-0. 6898044 1E+02-0. 2095 1885E+01  0.81214599E-01-0.51909467E-04  3 
0.16481 108E-07-0.20236530E-1 1  0.83749462E+05  0.3483 1424E+02  4 

MOL7  H  11C  5N  2  G  298.150  5000.000  1500.000  1 

0. 16535 182E+02  0.31575764E-01-0.1 1778552E-04  0.20360063E-08-0.13362984E-12  2 
0.10 148 196E+06-0.63303688E+02-0.16266850E+00  0.72 173440E-0 1-0. 48446302E-04  3 
0.16586150E-07-0.22675481E-1 1  0.10693343E+06  0.25497909E+02  4 

MOL8  H  13C  4N  20  1  G  298.150  5000.000  1500.000  1 

0.17374545E+02  0.35902592E-01-0.13305052E-04  0.22871019E-08-0.14941239E-12  2 
0.46989499E+05-0.70084578E+02-0.19549970E+01  0.80494649E-01-0.50943720E-04  3 
0.15925004E-07-0.19073303E-1 1  0.53570618E+05  0.33614265E+02  4 

MOL9  H  16C  6N  1  G  298.150  5000.000  1500.000  1 

1.82134750E+01  4.32403920E-02  -1.60232010E-05  2.75054700E-09-1.79354230E-13  2 

3.96208820E+04-7.48831380E+01-3.15401830E+00  8.67997730E-02-4.61619800E-05  3 
1.01808480E-08-4.236191 10E-13  4.75409470E+04  4.18991630E+01  4 

MOL  10  H  14C  5N  1  G  298.150  5000.000  1500.000  1 

1.56143130E+01  3.74427360E-02-1. 384085 10E-05  2.37154620E-09-1.54421660E-13  2 
4.51906300E+04-6.14899080E+01-1.81484190E+00  7.221 15700E-02-3.69007720E-05  3 
7.41 64381 0E-09- 1 .279 1 0030E- 1 3  5.17366300E+04  3.40534480E+01  4 

MOL11  H  13C  5N  1  G  298.150  5000.000  1500.000  1 

1.59482810E+01  3.47910660E-02-1.29295800E-05  2.2231 1060E-09-1.45094350E-13  2 
- 1 .84222690E+04-6.30692 1 1 0E+0 1  -9.59869520E-0 1  6.71 093780E-02-3 .247780 1 0E-05  3 
5.23561280E-09  2.98754580E-13-1.19131880E+04  3.01472910E+01  4 


MOL  12  H  15C  6N  1  G  298.150  5000.000  1500.000  1 
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1.85598390E+01  4.05535670E-02-1.50935760E-05  2.59837650E-09-L69756920E-13  2 
-2.25364320E+04-7.60810400E+01-2.26125710E+00  8.18745290E-02-4.22125790E-05  3 

8.33928440E-09-7.457791 10E-14-L46923830E+04  3.81366730E+01  4 

MOL  13  H  13C  5N  4  G  298.150  5000.000  1500.000  1 

2.1 1309450E+01  3.83671830E-02-1.43254990E-05  2.47370000E-09-1.62055080E-13  2 
8.77922240E+04-8.53840660E+01-5.90900300E-01  8.54137440E-02-5.04936990E-05  3 
1 .37 136290E-08- 1 .22957220E-12  9.55325930E+04  3.22981 140E+01  4 

MOL  14  H  14C  6N  7  G  298.150  5000.000  1500.000  1 

0.29266684E+02  0.44981381E-01-0.16930399E-04  0.2941581 1E-08-0.19363340E-12  2 
0. 12 167229E+06-0. 1252 1235E+03-0.259 1 5707E+00  0. 1 123 1959E+00-0. 7287 1240E-04  3 
0.22738680E-07-0.26489004E-1 1  0.13181224E+06  0.33478473E+02  4 

MOL  15  H  15C  8N  2  G  298.150  5000.000  1500.000  1 

0.23633619E+02  0.44316180E-01-0.165881 15E-04  0.28701773E-08-0.18833317E-12  2 
0.5313761 1E+05-0. 101 19307E+03-0.50748197E+01  0.10959076E+00-0.70581446E-04  3 
0.21853323E-07-0.25165423E-1 1  0.63019284E+05  0.53179294E+02  4 

MOL  16  H  11C  7N  2  G  298.150  5000.000  1500.000  1 

2.01407130E+01  3.41313630E-02-1.28330030E-05  2.22960970E-09-1.46829910E-13  2 
7.56107620E+04-8.22553380E+01-4.22537250E+00  9.34198090E-02-6.64336670E-05  3 
2.35242240E-08-3.27461520E-12  8.35604630E+04  4.73094550E+01  4 

MOL  17  H  14C  9N  1  G  298.150  5000.000  1500.000  1 

0.22887047E+02  0.42541565E-01-0.16001824E-04  0.27798010E-08-0.18299558E-12  2 
0.57807766E+05-0.9861 1 143E+02-0.70273065E+01  0.1 1350452E+00-0.78156131E-04  3 
0.2648898 1E-07-0.34820330E-1 1  0.67773046E+05  0.61 140793E+02  4 


MOL  18  H  10C  8N  1  G  298.150  5000.000  1500.000  1 
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0.19385432E+02  0.32357553E-01-0.12244457E-04  0.2 1 384603E-08-0. 14 142 1 1  IE-12  2 
0.80444642E+05-0.8032 1 39 1E+02-0.609669 1 9E+0 1  0.969362 1 8E-0 1  -0.73449456E-04  3 
0.27841274E-07-0.41753222E-1 1  0.884688 17E+05  0.542 12366E+02  4 

MOL  19  H  20C  9N  1  G  298.150  5000.000  1500.000  1 

2.51021820E+01  5.63194830E-02-2.10361390E-05  3.63349540E-09-2.38086630E-13  2 
3.65946780E+04-1.14534790E+02-8.25590260E+00  1.29996640E-01-7.94355540E-05  3 
2.27986850E-08-2.30060050E-12  4.83208490E+04  6.56533820E+01  4 

MOL20  H  16C  8N  1  G  298.150  5000.000  1500.000  1 

0.21607868E+02  0.46141348E-01-0.17284782E-04  0.29936958E-08-0.19663707E-12  2 
0.58673424E+05-0.95582664E+02-0.74028215E+01  0.1 1382253E+00-0.75285304E-04  3 
0.24469191E-07-0.30588004E-1 1  0.68465705E+05  0.597702 14E+02  4 

MOL21  H  14C  7N  3  G  298.150  5000.000  1500.000  1 

0.23383677E+02  0.41886809E-01-0.15683501E-04  0.27136129E-08-0.17803078E-12  2 
0.63734 166E+05-0.99096867E+02-0.44766222E+01  0.10463770E+00-0.66891 146E-04  3 
0.20342236E-07-0.22610795E-1 1  0.73390870E+05  0.50938069E+02  4 

MOL22  H  10C  6N  3  G  298.150  5000.000  1500.000  1 

0.1988181 1E+02  0.31710093E-01-0.1 1931202E-04  0.20734749E-08-0.13655142E-12  2 
0.86239 1 1 9E+05-0.80433040E+02-0.364 1 749 1E+0 1  0.88471 133E-01-0.62723789E-04  3 
0.21995312E-07-0.30147499E-1 1  0.93967638E+05  0.44830692E+02  4 

MOL23  H  2B  1C  2N  2  G  298.150  5000.000  1500.000  1 

1.0461 1980E+01  8.23 1 96630E-03-3 .244532 1 0E-06  5.83262300E-10-3.93984190E-14  2 
-1.19524750E+04-2.73473330E+01  3.47564280E+00  2.58450550E-02-1.98425640E-05  3 
7.50902870E-09- 1 . 1 1 850340E- 12-9.7426 1 500E+03  9.56697690E+00  4 


MOL24  H  9C  7N  2  G  298.150  5000.000  1500.000  1 
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1.93946860E+01  2.95388380E-02-1.1 1578610E-05  1.94645600E-09-1.28625530E-13  2 
LO 1429660E+05-7. 592706 10E+01 -2. 14987260E+00  8.39291620E-02-6.24863520E-05  3 
2.33983040E-08-3.477131 10E-12  L08237480E+05  3.78967370E+01  4 

MOL25  H  11C  6N  2  G  298.150  5000.000  1500.00  1 

0.18124559E+02  0.33085667E-01-0.12386017E-04  0.21439698E-08-0.14075077E-12  2 
0. 621 38460E+05-0. 7255 1 170E+02-0. 3702838 1E+01  0.83822661E-01-0.55653614E-04  3 
0.18054281E-07-0.22391805E-1 1  0.69526986E+05  0.44404675E+02  4 

MOL26  H  13C  7N  20  1  G  298.150  5000.000  1500.00  1 

0.23062 140E+02  0.39492745E-01-0.14793921E-04  0.25609108E-08-0.16808653E-12  2 
0.43554033E+05-0.97103934E+02-0.39133912E+01  0.10105465E+00-0.65982974E-04  3 

0.20701446E-07-0.24231499E-1 1  0.52813638E+05  0.47864867E+02  4 


Because  of  excessive  computational  time,  the  enthalpy  of  formation  and  therefore,  the 
thermodynamic  data  of  ions  pairs  or  neutral  pairs  were  not  calculated  in  this  work. 


43 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


3  Equilibrium  Calculations 

The  third  part  of  this  report  presents  results  obtained  from  the  equilibrium 
calculations.  At  this  time,  only  mixtures  which  were  found  to  be  hypergolic  experimentally 
are  considered  here. 


3.1  Initial  conditions 

The  determination  of  the  adiabatic  Flame  Temperature  (TF)  and  of  the  Adiabatic 
Isochoric  Complete  Combustion  Pressure  (Paicc)  was  carried  out  via  adiabatic  constant 
volume  calculations  with  the  "Equilibrium  Calculations"  suite  of  the  COSILAB  software 
package  [COSILAB,  2010], 


The  species  taken  into  account  for  the  equilibrium  calculations  are  listed  below: 

-  For  H,  C,  N  and  O: 

H,  C,  N,  O,  H+,  OH-,  H2,  CO,  N204,  C02,  H20,  N2,  03,  OH,  NO,  N02,  N20,  H202, 

ho2,  ch4,  hno2,  ch,  hcn,  ch2,  ch3,  CN,  CN-,  n3-,  n3h,  hno,  nh3,  ch3o,  ch2oh, 
ch3oh,  c2h2,  c2h4,  c3h6,  c4h8,  hocn,  hcho,  ch3cho,  hno3,  c2h5oh,  c3h8,  c2h6, 

C4Hio,  hnco,  NOT,  N03- 

-  For  chlorine  (Cl): 

Cl,  CF,  HCL,  C102,  CIO,  HOC1,  Cl2 

-  For  boron  (B): 

B,  B203,  BH2(CN)2\  BO,  BH,  bh2,  bh3,  bh4,  bh5,  b2,  b2h,  b2h2,  B2H6,  b2o,  b2o2 


For  all  of  these  species,  the  thermodynamic  data  was  taken  from  the  database  of 
Burcat  [Goos  et  ah,  2013], 


All  calculations  were  carried  out  with  same  starting  conditions  (initial  temperature  of 
298.15  K  and  initial  pressure  of  1  Atm). 

For  the  mixtures  for  which  charged  species  are  involved,  the  initial  conditions  are: 
X(cation)  =  X(anion) 

X(cation)  +  X(anion)  +  X(oxidant)  =1 

where  X  represents  the  mole  fraction  of  the  concerned  species.  X  varies  from  0  to  0.5. 
For  mixture  where  neutral  species  are  considered: 
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X(neutral  1)  =  X(neutral  2) 

X(neutral  1)  +  X(neutral  2)+  X(oxidant)  =1 

All  of  the  mixtures  considered  are  hypothetical  and  it  is  tenable  that  hypergolic  ignition  is  also 
possible  due  to  species  present  in  the  gas  phase  but  fonned  by  reactions  taking  place  in  the 
liquids.  This  possibility  is  examined  in  parts  two  and  three  of  this  report.  Furthermore,  the 
temperature  of  both  the  liquid  and  the  gas  is  certainly  higher  than  298  K  but  neither  the 
temperature  nor  the  compositions  are  known  when  experimental  hypergolic  ignition  is 
observed. 

As  seen  in  the  first  part  of  this  report,  ionic  liquids  can  be  divided  into  two  classes, 
namely  pro  tic  and  aprotic  ILs.  The  results  obtained  for  each  of  these  two  classes  are  presented 
separately  below. 


3.2  Protic  Ionic  Liquid 


Protic  Ionic  Liquids  are  ionic  liquids  for  which  the  cation  can  transfer  a  proton  (H+)  to 
the  anion  to  form  two  neutral  compounds.  As  seen  in  the  first  part  of  this  report,  the  chemical 
vapor  above  the  ionic  liquid  can  be  composed  of  separated  ions,  separated  neutrals,  ion  pairs, 
and  pair  of  neutrals  or  a  mixture  of  all  of  these.  The  chemical  vapor  above  the  IL/oxidant 
liquid  system,  once  mixed,  can  also  be  composed  of  a  combination  of  all  of  these  as  well  as 
oxidant  species  (and  its  fragments).  The  nature  of  the  vapor  remains  to  be  determined 
experimentally,  both  qualitatively  and  quantitatively.  In  this  study,  three  different  cases  have 
been  considered  for  equilibrium  calculations: 

-  In  the  vapor,  the  cation  and  the  anion  both  separated  are  present. 


-  In  the  vapor,  two  neutral  molecules  both  separated  are  present. 

-  In  the  vapor,  the  cation,  the  anion  and  the  neutral  molecules  all  of  them  separated  are 
present. 


Two  different  protic  ionic  liquids  were  found  to  be  hypergolic  with  N2O4  in  the 
literature.  These  systems  are: 


CH, 

I  2  3 

CH, 


(NC6Hi6+)  +  N3"  hypergolic  with  N2O4 


(NCsHi4+)  +  N3‘  hypergolic  with  N2O4 
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The  first  mixture  considered  is  NCeHi6++  N3"  with  N2O4.  The  results  obtained  for  the 
mixtures  with  only  the  cation  and  the  anion  separated  in  the  vapor  are  given  in  Figure  12  and 
Table  16. 

As  described  in  the  second  part  of  this  report,  a  temperature  criterion  of  2700  K  was 
proposed  for  determining  the  flammability  limits  of  potential  hypergolic  systems.  In  this  case, 
the  lower  flammability  limit  (LFL)  and  the  upper  flammability  limit  (UFL)  are  determined, 
respectively,  for  a  mole  fraction  of  the  cation,  X(cation),  between  0.039  and  0.048  for  the 
LFL  and  0.310  and  0.315  for  the  UFL. 

It  can  be  seen,  that  at  the  beginning  with  an  increase  of  the  mole  fraction  of  the  cation, 
the  flame  temperature  (TF)  and  the  adiabatic  isochoric  complete  combustion  pressure  (Paicc) 
rise  until  they  reached  their  maximum  for  a  mole  fraction  of  cation  around  0.140.  At  this  point 
Tp  is  around  4200  K  and  the  Paicc  is  around  57  Atm.  Then  when  the  mole  fraction  continues 
to  increase,  TF  and  the  Paicc  decrease  over  the  remaining  composition  domain. 
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Figure  12:  Flame  temperature  and  Paicc  far  different  +  N3  /  N2O4  mixtures  at 

298.15  K  and  1  Atm. 
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Mixture 

Cation 

CHS 

H  —  N — CH3 — CH, 

1  2  3 

CHj 

nc6h16+ 

Anion 

N3' 

Oxidant 

N204 

Stoichiometric 

Equation 

1  [NC6H16+  N3']  +  5  N204  ->  6  C02  +  8  H20  +  7  N2 

Initial 

Conditions 

P=1  Atm  ;  T  =  298.15  K 

X(NC6H16+)  =  X(N3") 

X(NC6H16+)  +  X(N3-)  +  X(N204)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0.039 

(TF  =  2661  K  ;  PA1CC 
=  28.2  Atm) 

<X(NC6H16+)<  < 
X(N3‘)  < 

0.048 

(Tf  =  2969  K  ;  PAICC 
=  32.0  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.310 

(Tf  =  2720  K  ;  PAICC 
=  45.7  Atm) 

<X(NC6H16+)<  < 
X(Nf)  < 

0.315 

(Tf  =  2650  K  ;  PA1CC 
=  44.5  Atm 

Paicc  Max 
Composition 

57.9  Atm 

(X(Ca>  0.145) 

Tf  Max 
Composition 

4266  K 

(X(Ca+)=  0.130) 

Table  16:  Equilibrium  Calculations  results  for  the  NCf,Hif  +  N3'  /  N2O4  system.  X( Ca+)  is 

the  mole  fraction  of  cation. 


In  Figure  13  and  in  Table  17  are  the  results  if  we  consider  the  chemical  vapor  above 
the  ionic  liquid  composed  of  neutral  compounds.  The  mixture  can  be  represented  as: 


(NC6Hi5)  +  HN3  with  N204 


A  diminution  in  the  value  of  the  UFL  can  be  observed  compared  to  the  value  obtained 
with  the  mixtures  with  charged  species.  On  the  contrary,  no  change  is  noted  in  the  value  for 
LFL.  The  UFL  is  now  encompassed  between  0.259  and  0.265,  leading  to  a  reduction  of  the 
flammability  domain  compared  to  the  one  observed  for  the  charged  species.  The  maximum 
flame  temperature  and  Paicc  are  obtained  for  mixtures  with  a  mole  fraction  X(Ca+)  around 
0.130. 
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Molar  Fraction  of  NC6H15  (or  N3H) 


Figure  13:  Flame  temperature  and  Paicc  far  different  NC (fall 5  +  HN3/  N2O4  mixtures  at 

298.15  K  and  1  Atm. 


Mixture 

Neutral 

Hj  H, 

3  |  3 

„CH, 

h3c"  s 

nc.6H15 

Neutral 

hn3 

Oxidant 

N204 

Stoichiometric 

Equation 

1  [NC6H15  N3H]  +  5  N2O4  -»  6  C02  +  8  H20  +  7  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(NC6H15)  =  X(HN3) 

X(NC6H15)  +  X(HN3)  +  X(N204)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0.039 

(Tf  =  2534  K  ;  PA1CC 
=  26.8  Atm) 

<X(NC6H15)<  < 

X(HN3)  < 

0.048 

(Tf  =  2842  K  ;  Pajcc 
=  30.6  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.259 

(Tf  =  2722  K  ;  PA1CC 
=  43.7  Atm) 

<X(NC6H15)<  < 

X(HN3)  < 

0.265 

(Tf  =  2658  K  ;  PA1CC 
=  43.0  Atm) 

Paicc  Max 
Composition 

54.0  Atm 

(X(Ne)=  0.138) 

TF  Max 
Composition 

4055  K 

(X(Ne)=  0.130) 

Table  17:  Equilibrium  Calculations  results  for  the  NC(HI5  +  HN3/  N2O4  system.  X(Ne)  is  the 

mole  fraction  of  neutrals. 
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The  chemical  vapor  above  the  ionic  liquid  can  be  composed  of  a  mixture  of  charged 
and  neutrals  species.  The  mixture  is  then  supposed  to  be: 


-0-0 

h2  h, 

H-N-CH3-CH3 

H.C'  'N'  'CH. 

1 

u  ^CH2 

CH3 

(NC6Hi6+)  +  N3  + 

(NC6Hi5)  +  HN3  and  this 


mixture  may  be  hypergolic  with  N2O4 


In  our  study,  we  have  considered  that  the  mole  ratio  between  the  neutral  and  the 
charged  species  is  equal  to  one.  The  initial  conditions  are  the  following: 

X(NC6H16+)  =  X(N3‘)  =  X(NC6H15)  =  X(HN3) 

X(NC6H16+)  +  X(N3-)  +  X(NC6H15)  +  X(HN3)  +  X(N204)  =  1 


The  results  obtained  for  this  mixture  are  presented  in  Figure  14  and  in  Table  18.  The 
lower  flammability  limit  is  obtained  for  a  mole  fraction  of  NC6Hi6+  between  0.019  and  0.024, 
and  the  upper  flammability  limit  is  between  0.143  and  0.145.  The  maximum  flame 
temperature  of  the  system  is  equal  to  4164  K  and  is  obtained  for  a  mole  fraction  equal  to 
0.065. 
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Figure  14:  Flame  temperature  and  P mcc  for  different  NCfJiis  +  HNj  +  NC6Hi6+  +  Nf  / N2O4 

mixtures  at  298.15  K  and  1  Atm. 


Mixture 


Neutral 


h3c'c' 

h3c' 


N' 

I 

CH 


C 


2 


ch3 


Cation 


Neutral 

HN3 


?!> 

H  —  N  —  CH, —  CH, 
I  2  3 

CH, 


Anion 

n3‘ 


Oxidant 

N204 


NCfiH 


6^15 


NCfiH 


6n16 


Stoichiometric 

Equation 


1  [NC6Hi6+  N3-]  +1  [NC6H15  HN3]  +  10  N204  ->  12  C02  +  16  H20  +  14  N2 


Initial 

Conditions 


P=1  Atm  ;T  =  298.15  K 
X(NCfiH16+)  =  X(N3-)  =  X(NC6H15)  =  X(HN3) 


X(NC6H16+)  +  X(N3-)  +  X(NC6H15)  +  X(HN3)  +  X(N204)  =  1 


Lower 

0.019 

<X(NC6H15)< 

0.024 

Flammability 
Limit  (LFL) 

(Tf  =  2599  K  ;  PAICC 
=  27.5  Atm) 

<  X(HN3)  or  X(NC6H16+)  or 
X(N3-)  < 

(Tf  =  2907  K  ;  PAICC 

=  31.3  Atm) 

Upper 

0.143 

<X(NC6H15)< 

0.145 

Flammability 
Limit  (UFL) 

(Tf  =  2751  K;Paicc 
=  46.1  Atm) 

<  X(HN3)  or  X(NC6H16+)  or 
X(Nf)  < 

(Tf  =  2671  K  ;  PAICC  = 

44.9  Atm) 

Paicc  Max 
Composition 

55.9  Atm 

(X(Ne)=  0.069) 

TF  Max  Composition 

4164  K 

(X(Ne)=  0.065) 

Table  18:  Equilibrium  calculations  results  for  the  NCeflis  +  FIN 3  +  NC6Hi6+  +  Nf  /N2O4 

system. 


Like  the  above  system,  for  the  second  reported  hypergol  case,  three  different  test 
compositions  of  the  chemical  vapor  above  the  ionic  liquid  once  mixed  with  the  oxidizer 
(charged  species,  neutral  species,  and  a  mixture  of  both  charged  and  neutral  species)  were 
similarly  considered. 


The  results  for  the  charged  species  are  given  in  Table  19.  The  system  considered  is: 


(NC5Hi4+)  +  N3‘  with  N204 


The  lower  flammability  limit  is  obtained  at  a  mole  fraction  of  cation  (or  anion) 
between  0.039  and  0.048,  and  the  upper  flammability  limit  at  a  mole  fraction  of  cation  (or 
anion)  between  0.342  and  0.346.  For  a  mixture  with  neutral  species  NC5FL3  and  HN3  (results 
given  in  Table  20),  the  flammability  domain  is  smaller.  Indeed,  the  LFL  is  obtained  for  a  mole 
fraction  of  the  neutral  species  between  0.048  and  0.057  and  the  UFL  is  between  0.286  and 
0.291.  For  the  mixtures  with  the  charged  species,  the  maximum  flame  temperature  is  equal  to 
4310  K  for  a  mole  fraction  of  cation  equal  to  0.153  (see  Table  19).  For  the  mixture  with 
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neutral  species,  this  maximum  flame  temperature  drops  to  4064  K  for  almost  the  same  mole 
fraction.  This  behavior  is  also  observed  for  the  value  of  the  maximum  Paicc- 


Mixture 

Cation 

H-N  —  CH 

CH3  CH3 

nc5h14+ 

Anion 

n3- 

Oxidant 

n2o4 

Stoichiometric 

Equation 

1  [NC5H14+  N3-]  +  17/4  N204  ->  5  C02  +  7  H20  +  25/4  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(NC5H14+)  =  X(N3-) 

X(NC5H14+)  +  X(N3-)  +  X(N204)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0.039 

(Tf  =  2482  K  ;  PAICC  = 
26.0  Atm) 

<X(NC5H14+)<  < 
X(N3-)  < 

0.048 

(Tf  =  2797  K  ;  PAICC 
=  29.6  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.342 

(Tf  =  2749  K  ;  PAICC  = 
43.0  Atm) 

<X(NC5H14+)<  < 
X(N3‘)  < 

0.346 

(Tf  =  2690  K  ;  PAICC 
=  42.0  Atm 

Paicc  Max 

Composition 

57.0  Atm 

(X(Ca+)=  0.160) 

TF  Max 
Composition 

4310  K 

(X(Ca+)=  0.153) 

Table  19:  Equilibrium  calculations  results  for  the  NC5H14'  +  Nf /N2O4  system. 


Mixture 

Neutral 

ch3 

6h3 

nc5h13 

Neutral 

N3H 

Oxidant 

N204 

Stoichiometric 

Equation 

1  [NC5H13  N3H]  +  17/4  N204  ->  5  C02  +  7  H20  +  25/4  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(NC5H13)  =  X(N3H) 

X(NC5H13)  +  X(N3H)  +  X(N204)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0.048 

(Tf  =  2646  K  ;  PAICC  = 
28.0  Atm) 

<X(NC5H13)<  < 

X(N3H)  < 

0.057 

(Tf  =  2897  K  ;  PAICC 
=  31.0  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.286 

(Tf  =  2700  K  ;  PA1CC  = 
41.2  Atm) 

<X(NC5H13)<  < 

X(N3H)  < 

0.291 

(Tf  =  2642  K  ;  PAICC 
=  40.6  Atm) 
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Paicc  Max 

52.5  Atm 

Tf  Max 

4064  K 

Composition 

(X(Ne)=  0.153) 

Composition 

(X(Ne)=  0.145) 

Table  20:  Equilibrium  calculations  results  for  the  NC5H13  +  HN3/N2O4  system. 


The  results  obtained  for  a  mixture  of  charged  and  neutral  species  are  given  in  Table 
2 1 .  The  initial  conditions  are  given  below: 


9Hs  ch3 

cn3 

1  +  / 

H-N  —  CH 

I  \ 

CH3  CH3 

6h3 

(NC5H14+)  +  N3-  + 

(NC5H13)  +  HN3  hypergolic 


with  N2O4 


-  X(NC5H14+)  =  X(N3‘)  =  X(NC5H13)  =  X(HN3) 

-  X(NC5H14+)  +  X(NT)  +  X(NC5Hi3)  +  X(HN3)  +  X(N204)  =  1 


For  this  mixture,  the  lower  flammability  limit  is  obtained  for  a  mole  fraction  between 
0.019  and  0.024,  and  between  0.158  and  0.160  for  the  upper  flammability  limit. 

The  maximum  flame  temperature  is  obtained  for  a  mole  fraction  of  0.073  and  is  equal 
to  4191  K.  This  value  is  between  the  one  obtained  for  the  mixture  with  charged  species 
(4310  K)  and  the  one  obtained  for  the  mixture  with  neutral  species  (4064  K). 


Mixture 

Neutral 

Neutral 

hn3 

Cation 

Anion 

Nf 

Oxidant 

N204 

CH3 

h’c-n^ch3 

6h3 

W  7H3 

H-N  —  CH 

CH3  CH3 

nc5h13 

nc5h14+ 

Stoichiometric 

Equation 

1  [NC5H14+N3-]+l 

->  10  C02 

[NC5H13  HN3]  +  17/2  N204 

+  14  H20  +  25/2  N2 

Initial 

Conditions 

P=1  Atm;  T  =  298.15  K 

X(NC5H14+)  =  X(N3-)  =  X(NC5H13)  =  X(HN3) 

X(NC5H14+)  +  X(Nf)  +  X(NC5H13)  +  X(HN3)  +  X(N204)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0.019 

(Tf  =  2411  K  ;  PAICC  = 
25.2  Atm) 

<X(NC5H13)< 

<  X(HN3)  or  X(NC5H14+)  or 
X(N3-)  < 

0.024 

(Tf  =  2723  K  ; 

Paicc  =  28.8  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.158 

(Tf  =  2740  K  ;  PA1CC  = 
43.1  Atm) 

<X(NC5H13)< 

<  X(HN3)  or  X(NC5H14+)  or 
X(N3-)  < 

0.160 

(Tf  =  2660  K  ; 

Paicc  =  41.8  Atm) 

Paicc  Max 
Composition 

54.7  Atm 

(X(Ca+)=  0.076) 

Tf  Max  Composition 

4191  K 

(X(Ca+)=  0.073) 
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Table  21:  Equilibrium  calculations  results  for  the  NC5H14*  +  Nf  +  NC5H13  +  HN3/N2O4 

system. 


3.3  Aprotic  Ionic  Liquid 


For  these  ionic  liquids,  only  mixtures  in  which  the  cation  and  the  anion  are  separated 
in  the  vapor  are  considered. 


For  these  ionic  liquids,  the  results  will  be  presented  for  each  series  of  cations. 


3.3.1  N’N-dimethyltriazanium  cations 


The  first  ionic  liquid  tested  is  the  one  with  a  N’N-dimethyltriazanium  cation 


CH-j 

H_N— N— CH, 
2  !  a 

NH2 


(N3C2H10A  The  results  obtained  for  the  system  N3C2Hio+  /  N(CN)2"  which  is 


hypergolic  with  N2O4  are  presented  in  Figure  15  and  in  Table  22. 


The  results  show  that  the  flame  temperature  is  above  2700  K  for  the  entire  domain. 
According  to  the  criterion  proposed,  this  mixture  is  hypergolic  for  any  initial  composition 
(pure  IL  or  pure  oxidizer  excluded).  The  maximum  flame  temperature  is  equal  to  4450  K  and 
is  obtained  for  a  mole  fraction  equal  to  0.187.  The  maximum  value  for  the  Paicc  is  equal  to 
58.3  Atm  and  is  obtained  for  a  mole  fraction  equal  to  0.206. 


60 
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20 


Figure  15:  Flame  temperature  and  Paicc  for  different  N3C2H  uf  +  N(  CN)f  /  N2O4  mixtures  at 

298.15  K  and  1  Atm. 


Mixture 

Cation 

Anion 

Oxidant 
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ch3 

H_N— N—  CH, 

2  1  3 

NH2 

N3C2Hio+ 

N(CN)2- 

n2o4 

Stoichiometric 

Equation 

1  [N3C2H10+  N(CN)2‘]  +  13/4  N2O4 

->  4  C02  +  5  H20  +  25/4  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(N3C2H10+)  =  X(N(CN)2-) 

X(N3C2H,o+)  +  X(N(CN)2-)  +  X(N204)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<  X(N3C2Hio  )  <  < 
X(N(CN)2-)  < 

9.90E-03 

(Tf  =  3066  K  ;  PAICC 
=  25.7  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.497 

(Tf  =  2903  K  ;  PAICC  = 
48.9  Atm) 

<X(N3C2H10+)<  < 
X(N(CN)2-)  < 

0.500 

Paicc  Max 

Composition 

58.3  Atm 

(X(Ca+)=  0.206) 

Tf  Max 
Composition 

4450  K 

(X(Ca+)=  0.187) 

Table  22:  Equilibrium  calculations  results  for  the  NjC^H /o+  +  N(  CN)f  /  N2O4  system. 


Results  obtained  for  the  system  N3C2Hi0+  /  N(CN)2‘  which  is  hypergolic  with  HN03 
are  given  in  Table  23.  Like  the  previous  case,  the  flame  temperature  is  above  2700  K  for  the 
entire  domain  tested.  For  this  system  the  maximum  Paicc  is  obtained  for  a  mixture  with  a 
mole  fraction  of  cation  (or  anion)  equal  to  0.497.  The  maximum  flame  temperature  is  equal  to 
3902  K  and  is  obtained  when  the  mole  fraction  of  the  cation  is  0.1 15.  This  value  is  lower  than 
the  one  obtained  for  the  system  with  N2O4. 


Mixture 

Cation 

ch3 

H,N  — N— CH-i 

2  !  3 

NH2 

N3C2Hio+ 

Anion 

N(CN)2- 

Oxidant 

hno3 

Stoichiometric 

Equation 

1  [N3C2H10+  N(CN)2‘]  +  26/5  HN03 

->  4  C02  +  38/5  H20  +  28/5  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(N3C2H10+)  =  X(N(CN)2-) 

X(N3C2H10+)  +  X(N(CN)2')  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(N3C2H10+)<  < 
X(N(CN)2-)  < 

9.90E-03 

(Tf  =  2724  K  ;  PAICC 
=  20.3  Atm) 

Upper 

Flammability 

0.497 

<  X(N3C2Hio  )  <  < 

0.500 

54 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


Limit  (UFL) 

(TF  =  2895  K  ;  PAICC  = 
48.7  Atm) 

X(N(CNV)  < 

Paicc  Max 
Composition 

48.7  Atm 

(X(Ca+)=  0.497) 

Tf  Max 
Composition 

3902  K 

(X(Ca+)=  0.1 15) 

Table  23:  Equilibrium  calculations  results  for  the  NjC2Hio+  +  N(CN)f  /  HNO3  system. 


The  anion  N(CN)2_  was  replaced  by  NO3".  For  the  ionic  liquid  N3C2Hio+  /  NO3"  which 
is  hypergolic  with  HNO3,  the  results  are  presented  in  Table  24.  On  the  contrary  to  the 
previous  cases  just  presented  above,  this  mixture  is  not  hypergolic  for  the  entire  domain.  The 
lower  flammability  limit  is  obtained  for  a  mole  fraction  of  cation  between  9.90E-03  and 
0.020.  However,  the  upper  flammability  limit  is  between  0.497  and  0.500  for  the  cation  mole 
fraction.  The  maximum  flame  temperature  is  equal  to  3904  K  for  a  mole  fraction  of  cation  (or 
anion)  equal  to  0.160. 

The  results  obtained  for  the  N3C2H10  NO3"  ionic  liquid  with  N2O4  are  presented  in 
Table  25.  For  the  entire  domain,  the  adiabatic  flame  temperature  at  constant  volume  is  above 
2700  K.  As  seen  before  for  the  mixture  with  the  N(CN)2"  anion,  higher  flame  temperatures  are 
obtained  with  N2O4  compared  to  HNO3.  For  this  mixture,  the  maximum  flame  temperature  is 
equal  to  4344  K  which  is  400  K  higher  than  the  one  obtained  for  the  mixture  with  HNO3. 


Mixture 

Cation 

ch3 

H-N-N  —  CH, 

2  1  3 

nh2 

N3C2Hio+ 

Anion 

N03‘ 

Oxidant 

HNO3 

Stoichiometric 

Equation 

1  [N3C2H10+  NOf]  +  12/5  HN03 

->  2  C02  +  3 1/5  H20  +  16/5  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(N3C2H10+)  =  X(NO3-) 

X(N3C2H10+)  +  X(N03‘)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

9.90E-03 

(Tf  =  2675  K  ;  PAICC  = 
19.8  Atm) 

<X(N3C2H10+)<  < 
X(N03-)  < 

0.020 

(Tf  =  2769  K  ;  PAICC 
=  21.2  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.497 

(Tf  =  3385  K  ;  PAICC  = 

54. 1  Atm) 

<  X(N3C2Hio  )  <  < 
X(NCV)  < 

0.500 

Paicc  Max 
Composition 

54.1  Atm 

(X(Ca+)=  0.497) 

Tf  Max 
Composition 

3904  K 

(X(Ca+)=  0.160) 

Table  24:  Equilibrium  calculations  results  for  the  NjC^fTio*  +  NOf  /  HNO3  system. 
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Mixture 

Cation 

ch3 

H-N-N-  CH, 

2  1  3 

nh2 

N3C2Hio+ 

Anion 

NOf 

Oxidant 

N204 

Stoichiometric 

Equation 

1  [N3C2H10+  NOf]  +  3/2  N2O4  ->  2  C02  +  5  H20  +  7/2  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(N3C2H10+)  =  X(NO3-) 

X(N3C2H10+)  +  X(N03-)  +  X(N204)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<  X(N3C2H10+)  <  < 
X(NOf)  < 

9.90E-03 

(Tf  =  3036  K  ;  PAICC 
=  25.2  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.497 

(Tf  =  3390  K  ;  PAICC  =  54.2 
Atm) 

<  X(N3C2Hio  )  <  < 
X(NOf)  < 

0.500 

Paicc  Max 
Composition 

58.2  Atm 

(X(Ca+)=  0.270) 

TF  Max 
Composition 

4344  K 

(X(Ca+)=  0.242) 

Table  25:  Equilibrium  calculations  results  for  the  N3C2H uf  +  NOf  /  N2O4  system. 


The  N'N-dimethyltriazanium  cation  was  also  tested  with  the  chloride  anion  (Cl  ).  This 
ionic  liquid  is  reported  hypergolic  with  N2O4.  The  results  obtained  for  this  mixture  are  given 
in  Table  26. 
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For  this  mixture,  the  flame  temperature  is  above  2700  K  for  the  entire  domain.  The 
maximum  Paicc  and  flame  temperature  are  equal,  respectively,  to  55.7  Atm  and  4509  K  for 
mole  fractions  of  cation  equal  to  0.224  and  0.206. 

In  summary,  for  this  cation,  hypergolicity  is  predicted  for  almost  the  entire  domain  for 
all  the  systems  tested. 


Mixture 

Cation 

CHa 

H_N  — -N— CH, 

2  !  3 

nh2 

N3C2Hio+ 

Anion 

cr 

Oxidant 

n2o4 

Stoichiometric 

Equation 

1  [N3C2H10+  Cl  ]  +  1 1/4  N204 

2  C02  +  1  C102  +  5  H20  +  17/4  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(N3C2H10+)  =  X(Cl') 

X(N3C2H10+)  +  X(Cl')  +  X(N204)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(N3C2H10+)<  < 
X(Cl')  < 

0.020 

(Tf  =  3107  K  ;  PAICC 
=  26.6  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.497 

(TF  =  3249  K  ;  PAICC  = 
47.1  Atm) 

<X(N3C2H10+)<  < 
X(Cl')  < 

0.500 

Paicc  Max 

Composition 

55.7  Atm 

(X(Ca)=  0.224) 

TF  Max 
Composition 

4509  K 

(X(Ca+)=  0.206) 

Table  26:  Equilibrium  calculations  results  for  the  NjC2Hio+  +  Cl~ /N2O4  system. 


3.3.2  N’N-dimethylhydrazinium  cations 

In  this  part,  the  N'N-dimethylhydrazinium  cation  with  different  alkyl  substituents  was 

tested. 


The  first  ionic  liquid  tested,  shown  below,  is  reported  hypergolic  with  HNO3. 


57 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


CH3 

H2N-N-CH3 

CH3 


(N2C3Hn+)  +  N(CN)2- 


For  this  case,  the  results  are  given  in  Figure  16  and  in  Table  27.  Flypergolicity  is 
predicted  until  a  mole  fraction  of  cation  (or  anion)  between  0.371  and  0.375  (the  upper 
flammability  limit).  The  maximum  flame  temperature  is  found  for  a  mole  fraction  of  cation 
equal  to  0.099  and  is  equal  to  3795  K. 


Figure  16:  Flame  temperature  and  Paicc  for  different  N2C3H.11*  +  N(CN)2  / HNO3  mixtures 

at  298.15  K  and  1  Atm. 


Mixture 

Cation 

CH3 

H2N-N-CH3 

CH3 

N2C3Hh+ 

Anion 

N(CN)2- 

Oxidant 

hno3 

Stoichiometric 

Equation 

1  [N2C3Hn+  N(CN)2"]  +  93/15  HN03 
->  5  C02  +  43/5  H20  +  84/15  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(N2C3H„+)  =  X(N(CN)2-) 

X(N2C3H„+)  +  X(N(CN)2-)  +  X(HN03)  =  1 
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Lower 

Flammability 
Limit  (LFL) 

0 

<X(N2C3H„+)<  < 
X(N(CN)2")  < 

9.90E-03 

(Tf  =  2734  K  ;  PAICC  = 
20.4  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.371 

(Tf  -2701  K  ;  PAICC 
=  40.9  Atm) 

<X(N2C3H„+)<  < 
X(N(CN)2-)  < 

0.375 

(Tf  =  2686  K  ;  PAICC  = 
40.9  Atm) 

Paicc  Max 

41.0  Atm 

Tf  Max 

3795  K 

Composition 

(X(Ca+)=  0.320) 

Composition 

(X(Ca+)=  0.099) 

Table  27:  Equilibrium  calculations  results  for  the  fyCsHu*  +  N(CN)2  /  HNO3  system. 


In  Table  28  are  the  results  for  the  ionic  liquid  N2C3Hn+  /  BFbfCNT"  which  was  found 
hypergolic  with  HNO3.  Compared  to  the  above  case,  the  dicyanamide  anion  was  replaced  by 
the  dicyanoborate  anion.  For  this  system,  the  upper  flammability  limit  was  found  for  a  mole 
fraction  of  cation  between  0.291  and  0.296  (with  Paicc  around  38  Atm).  The  presence  of  the 
dicyanoborate  anion  instead  of  the  dicyanamide  anion  leads  to  a  decrease  of  the  flammability 
domain  for  this  IL. 


Mixture 

Cation 

CH3 

H2N-N-CH3 

CH3 

N2C3Hn+ 

Anion 

BH2(CN)2- 

Oxidant 

HN03 

Stoichiometric 

Equation 

1  [N2C3Hu+  BH2(CN)2]  +  36/5  HNO3 

->  5  C02  +  1/2  B2C>3  +101/10  H20  +  56/10  N2 

Initial 

Conditions 

P=1  Atm;  T  =  298.15  K 

X(N2C3H„+)  =  X(BH2(CN)2-) 

X(N2C3H„+)  +  X(BH2(CN)2)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(N2C3H„+)<  < 
X(BH2(CN)2-)  < 

9.90E-03 

(Tf  =  2746  K  ;  PAICC  = 
20.6  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.291 

(Tf  =  2716K;PAICc 
=  37.9  Atm) 

<X(N2C3H„+)<  < 
X(BH2(CN)2-)  < 

0.296 

(Tf  =  2694  K  ;  PAICC  = 
37.8  Atm) 

Paicc  Max 
Composition 

38.1  Atm 

(X(Ca+)=  0.237) 

TF  Max 
Composition 

3716  K 

(X(Ca+)=  0.083) 

Table  28:  Equilibrium  calculations  results  for  the  N2C;i/T  if  +  BH2(  CN)2~  /  N2O4  system. 
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In  Table  29  and  30  are  the  results  for  the  1-butyl-N'N-dimethylhydrazinium  cation 
based  ionic  liquids.  The  first  system  tested  is  with  the  dicyanamide  anion  (Table  29)  and  the 
second  one  with  the  dicyanoborate  anion  (Table  30).  Both  of  these  ionic  liquids  are 
hypergolic  with  WFNA. 

For  these  two  systems,  hypergolicity  (TF  above  2700  K)  is  predicted.  For  the  mixture 
with  the  dicyanamide  anion,  the  upper  flammability  limit  is  found  for  a  mole  fraction  of 
cation  between  0.194  and  0.200.  As  seen  before,  the  presence  of  the  dicyanoborate  anion 
induces  smaller  flammability  domain.  The  upper  flammability  limit  is  found  for  mole  fraction 
of  cation  between  0.174  and  0.180. 

For  these  two  cases,  the  maximum  flame  temperature  is  around  3600  K  for  small 
values  of  mole  fraction  of  cation. 


Mixture 

Cation 

<??> 

HaN  -N-CHj 
„CH, 

h2c-  3 

_ch2 
h3c  3 

n2c6h17+ 

Anion 

N(CN)2- 

Oxidant 

hno3 

Stoichiometric 

Equation 

1  [N2C6H17+N(CN)21  +  49/5  HN03 
->  8  C02  +  67/5  H20  +  37/5  N2 

Initial 

Conditions 

P=1  Atm;  T  =  298.15  K 

X(N2C6H17+)  =  X(N(CN)2-) 

X(N2C6H17+)  +  X(N(CN)2')  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(N2C6H17+)<  < 
X(N(CN)2')  < 

9.90E-03 

(Tf  =  2790  K  ;  PAICC  = 
21.4  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.194 

(Tf  =  2723  K  ;  PAICC 
=  39.6  Atm) 

<X(N2C6H17+)<  < 
X(N(CN)2')  < 

0.200 

(Tf  =  2683  K  ;  PAICC  = 
39.6  Atm) 

Paicc  Max 
Composition 

39.7  Atm 

(X(Ca+)=  0.174) 

Tf  Max 
Composition 

3646  K 

(X(Ca+)=  0.065) 

Table  29:  Equilibrium  calculations  results  for  the  N2C(,Hn+  +  N( CN)f  /  HNO3  system. 
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Mixture 

Cation 

HaN  -N-CHj 

,CHj 

H2c-  3 

.ch2 
h3c  2 

N2C.6H17+ 

Anion 

BH2(CN)2- 

Oxidant 

hno3 

Stoichiometric 

Equation 

1  [N2C6H17+  BH2(CN)2']  +  54/5  HN03 

->  8  C02  +  1/2  B203  +  149/10  H20  +  37/5  N2 

Initial 

Conditions 

P=1  Atm;  T  =  298.15  K 

X(N2C6H17+)  =  X(BH2(CN)2-) 

X(N2C6Hn+)  +  X(BH2(CN)2)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(N2C6H17+)<  < 
X(BH2(CN)2-)  < 

9.90E-03 

(Tf  =  2803  K  ;  PA1CC  = 
21.6  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.174 

(Tf  =  2710  K  ;  PAicc 
=  38.1  Atm) 

<X(N2C6H17+)<  < 
X(BH2(CN)2-)  < 

0.180 

(Tf  =  2665  K  ;  PAICC  = 
38.0  Atm) 

Paicc  Max 
Composition 

38.2  Atm 

(X(Ca+)=  0.145) 

TF  Max 
Composition 

3605  K 

(X(Ca+)=  0.057) 

Table  30:  Equilibrium  calculations  results  for  the  N2C(,H]tV  +  BHjf  CN)f  /  HNO3  system. 


1-allyl-N'N-dimethylhydrazinium  cation  based  ionic  liquids  are  reported  hypergolic 
with  WFNA  with  two  different  anions:  the  nitrocyanamide  anion  (see  Table  31)  and  the 
dicyanoborate  anion  (see  Table  32). 

These  two  mixtures  are  predicted  hypergolic.  For  the  dicyanamide  anion  based  ILs, 
the  upper  flammability  limit  is  predicted  to  be  equal  to  a  mole  fraction  of  cation  between 
0.291  and  0.296.  As  seen  for  the  other  N'N-dimethylhydrazinium  cation  based  ILs,  the 
replacement  of  the  N(CN)2_  by  BFLICNL"  leads  to  a  decrease  of  the  flammability  domain 
(UFL  :  0.242  <  X(cation)  <  0.248  ;  see  Table  32). 
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Mixture 

Cation 

ch3 

h,n-n-ch, 

3  1  3 
^CH, 

HC  3 

II 

CH, 

N2C-5Hi3+ 

Anion 

N(CN)2- 

Oxidant 

hno3 

Stoichiometric 

Equation 

1  [N2C5H13+  N(CN)2~]  +  41/5  HNO3 

->  7  C02  +  53/5  H20  +  33/5  N2 

Initial 

Conditions 

P=1  Atm;  T  =  298.15  K 

X(N2C5H13+)  =  X(N(CN)2-) 

X(N2C5H13+)  +  X(N(CN)2‘)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(N2C5H13+)<  < 
X(N(CN)2')  < 

9.90E-03 

(TF  =  2772  K  ;  PAICC  = 
21.0  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.291 

(Tf  =  2719K;Paicc 
=  42.2  Atm) 

<  X(N2C5H13+)  < 
X(N(CN)2-) 

0.296 

(TF  =  2697  K  ;  PAICC  = 
42.2  Atm) 

P\icc  Max 
Composition 

42.2  Atm 

(X(Ca+)=  0.286) 

Tf  Max 
Composition 

3767  K 

(X(Ca+)=  0.083) 

Table  31:  Equilibrium  calculations  results  for  the  /V2C5//7/'  +  N(CN)f  /  HNO3  system. 
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Mixture 

Cation 

CHa 

FLN-N— CH_ 

3  1  3 
^CH, 

HC  2 

II 

CK, 

n2c5h13+ 

Anion 

BH2(CN)2- 

Oxidant 

hno3 

Stoichiometric 

Equation 

1  [N2C5H13+  BH2(CN)2‘]  +  46/5  HN03 

7  C02  +  1/2  B203  +  121/10  H20  +  33/5  N2 

Initial 

Conditions 

P=1  Atm;  T  =  298.15  K 

X(N2C5H13+)  =  X(BH2  (CNV) 

X(N2C5H13+)  +  X(BH2  (CNV)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(N2C5H13+)<  < 
X(BH2  (CNV)  < 

9.90E-03 

(Tf  =  2784  K  ;  PAICC  = 
21.2  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.242 

(Tf  =  2716K;PAICc 
=  39.6  Atm) 

<X(N2C5H13+)<  < 
X(BH2  (CNV)  < 

0.248 

(Tf  =  2687  K  ;  PAICC  = 
39.6  Atm) 

Paicc  Max 
Composition 

39.6  Atm 

(X(Ca+)=  0.225) 

TF  Max 
Composition 

3700  K 

(X(Ca+)=  0.074) 

Table  32:  Equilibrium  calculations  results  for  the  N2CsHj3+  +  BH2( CN)f  /  HNO3  system. 


l-(2-propargyl)-N'N-dimethylhydrazinium  cation  based  ionic  liquids  were  found 
hypergolic  with  WFNA  with  two  different  anions:  the  nitrocyanamide  anion  (see  Table  33) 
and  the  dicyanoborate  anion  (see  Table  34).  For  these  two  cases,  the  systems  are  predicted 
hypergolic.  For  the  dicyanamide  anion  based  IL,  the  upper  flammability  limit  was  found  for  a 
mole  fraction  of  cation  between  0.415  and  0.419. 

For  the  N2CsFIii+  /  BFI2(CN)2'  ionic  liquid,  the  upper  limit  was  found  for  a  mole 
fraction  of  cation  between  0.315  and  0.320.  For  this  cation,  replacement  of  the  N(CN)2"  anion 
by  the  BH2(CN)2‘  anion  again  induces  a  decrease  of  the  flammability  domain. 
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Mixture 

Cation 

CH, 

1 * 

HaN  — N— CH, 

2  1  3 
CHa 

Ilf 

CH 

N2CsHii+ 

Anion 

N(CN)2- 

Oxidant 

HN03 

Stoichiometric 

Equation 

1  [N2C5H11+N(CN)21  +  39/5  HNO3 
->  7  C02  +  47/5  H20  +  32/5  N2 

Initial 

Conditions 

P=1  Atm;  T  =  298.15  K 

X(N2C5H„+)  =  X(N(CN)2-) 

X(N2C5Hh+)  +  X(N(CN)2‘)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(N2C5H„+)<  < 
X(N(CN)2-)  < 

9.90E-03 

(Tf  =  2775  K  ;  PAICC  = 
20.9  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.415 

(Tf  =  2715K;Paicc 
=  46.7  Atm) 

<X(N2C5H„+)<  < 
X(N(CN)2‘)  < 

0.419 

(Tf  =  2692  K  ;  PAICC  = 
46.4  Atm) 

Paicc  Max 
Composition 

48.4  Atm 

(X(Ca+)=  0.398) 

Tf  Max 
Composition 

3892  K 

(X(Ca+)=  0.091) 

Table  33:  Equilibrium  calculations  results  for  the  N2C5H1S  +  N(CN)f  /  HNO3  system. 
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Mixture 

Cation 

CH, 

1  * 

H,N-N— CH, 

3  1  3 
CHa 

Ilf 

CH 

N2CsHn+ 

Anion 

BH2(CN)2- 

Oxidant 

HN03 

Stoichiometric 

Equation 

1  [N2C5H„+  BH2(CN)2']  +  44/5  HN03 

->  7  C02  +  1/2  B203  +  109/10  H20  +  32/5  N2 

Initial 

Conditions 

P=1  Atm;  T  =  298.15  K 

X(N2C5H„+)  =  X(BH2(CN)2-) 

X(N2C5H„+)  +  X(BH2  (CN)2-)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(N2C5H„+)<  < 
X(BH2(CN)2-)  < 

9.90E-03 

(Tf  =  2787  K  ;  PAICC  = 
21.1  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.315 

(Tf  =  2720  K  ;  PAICC 
=  42.7  Atm) 

<X(N2C5H„+)<< 
X(BH2(CN)2-)  < 

0.320 

(Tf  =  2688  K  ;  PAICC  = 
42.5  Atm) 

Paicc  Max 
Composition 

42.8  Atm 

(X(Ca+)=  0.296) 

TF  Max 
Composition 

3792  K 

(X(Ca+)=  0.083) 

Table  34:  Equilibrium  calculations  results  for  the  N2C2H/2+  +  BH2( CN)f  /  HNO3  system. 


The  last  N'N-dimethylhydrazinium  cation  under  consideration  in  this  study  is  the  l-(2- 
hydroxyethyl)-N'N-dimethylhydrazinium  cation.  Together  with  a  dicyanamide  anion,  this 
ionic  liquid  was  reported  hypergolic  with  nitric  acid  (HNO3). 

The  system  is  predicted  to  be  hypergolic.  The  upper  flammability  limit  was  found  for 
a  mole  fraction  of  cation  between  0.286  and  0.291  (see  Table  35). 


65 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


Mixture 

Cation 

CH, 

1 . 

HaN-N— CH, 

2  1  3 
^ch2 
h3c" 

OH 

N2C4OH13+ 

Anion 

N(CN)2- 

Oxidant 

hno3 

Stoichiometric 

Equation 

1  [N2C40H13+  N(CN)2']  +  7  HN03  ->  6  C02  +10  H20  +  6  N2 

Initial 

Conditions 

P=1  Atm;  T  =  298.15  K 

X(N2C40H13+)  =  X(N(CN)2-) 

X(N2C4OHi3+)  +  X(N(CN)+)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<  X(N2C4OH13+)  <  < 
X(N(CN)2-)  < 

9.90E-03 

(Tf  =  2746  K  ;  PAICC  = 
20.7  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.286 

(Tf  =  2705  K  ;  PAICC 
=  41.4  Atm) 

<  X(N2C4OH13+)  <  < 
X(N(CN)2‘)  < 

0.291 

(Tf  =  2683  K  ;  PAICC  = 
41.4  Atm) 

Paicc  Max 
Composition 

41.4  Atm 

(X(Ca>  0.281) 

Tf  Max 
Composition 

3695  K 

(X(Ca+)=  0.091) 

Table  35:  Equilibrium  calculations  results  for  the  N2C.4OH  n+  +  N(CN)f  /  HNO3  system. 


For  the  N'N-dimethylhydrazinium  cations,  different  alkyl  substituents  were  considered 
(methyl,  butyl,  allyl,  propargyl,  and  hydroxy  ethyl).  For  the  ionic  liquid  with  a  dicyanamide 
anion,  the  impact  on  the  flammability  domains  of  these  alkyl  substituents  in  the  cation  is  as 
follows:  butyl  <  hydroxy  ethyl  <  allyl  <  methyl  <  propargyl. 

3.3.3  Cations  with  an  azide  function 


Different  ionic  liquids  with  a  cation  containing  an  azide  function  were  studied  and 
reported  hypergolic  with  HNO3. 


The  first  system  tested  is  given  below: 


CH3 

1  + 

H.C-N  — 


CH, 


C— C 

Hj  H 


N=N=N 


(N4C5H13+)  +  N(CN)2‘  with  HN03 


The  results  obtained  for  this  are  given  in  Figure  17  and  in  Table  37.  Again, 
hypergolicity  is  predicted.  When  the  mole  fraction  of  the  cation  increases,  the  flame 
temperature  rises  until  it  reaches  a  maximum  (3771  K)  at  a  mole  fraction  of  0.083.  Then,  the 
flame  temperature  decreases  and  the  upper  flammability  is  predicted  for  a  mole  fraction  of 
cation  between  0.310  and  0.315. 
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A  second  similar  ionic  liquid  is  also  reported  hypergolic  with  HNO3.  The  ionic  liquid 


is: 


Cation: 


(N7C6H14+) ;  Anion:  N(CN)2' 


The  trends  predicted  for  the  previous  system  are  also  predicted  here.  The  results  for 
this  system  are  given  in  Table  37. 


50 


40 


o 

o 

< 

Q. 


30 


20 


Figure  17:  Flame  temperature  and  PAiccfor  different  N4.C5H.13*  +  N(CN)2  / HNO3  mixtures 

at  298.15  K  and  1  Atm. 
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Mixture 

Cation 

|Ha 

H.C-N— C-C—  N=N=N 

V’ H’ 

N4C-5Hi3+ 

Anion 

N(CN)2- 

Oxidant 

hno3 

Stoichiometric 

Equation 

1  [N4C5H13+  N(CN)2-]  +  41/5  HNO3 

->  7  C02  +  53/5  H20  +  38/5  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(N4C5H13+)  =  X(N(CNV) 

X(N4C5H13+)  +  X(N(CN)2‘)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<  X(N4C5H13+)  <  < 
X(N(CN)2-)  < 

9.90E-03 

(Tf  =  2777  K  ;  PAICC  = 
21.1  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.310 

(Tf  =  2710  K  ;  PAICC 
=  46.3  Atm) 

<  X(N4C5H13+)  <  < 
X(N(CN)2-)  < 

0.315 

(Tf  =  2692  K  ;  PAICC  = 
46.3  Atm) 

Paicc  Max 

Composition 

46.6  Atm 

(X(Ca)=  0.379) 

TF  Max 
Composition 

3771  K 

(X(Ca+)=  0.083) 

Table  36:  Equilibrium  calculations  results  for  the  N4C.5H n+  +  N(CN)f  /  HNO3  system. 


Mixture 

Cation 

Anion 

Oxidant 
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®HJ 

H.C — N  * — C — C N=N  =  N 

1  H,  H, 

<?H, 

N 

11 

N 

11 

N 

n7c6h14+ 

N(CN)2- 

HN03 

Stoichiometric 

Equation 

1  [N7C6H14+  N(CN)2"]  +  46/5  HN03 

->  8  C02  +  58/5  H20  +  48/5  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(N7C6H14+)  =  X(N(CN)2-) 

X(N7C6IIi4:)  +  X(N(CN)2-)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(N7C6H14+)<  < 
X(N(CN)2-)  < 

9.90E-03 

(TF  =  2810K;PAICC  = 
21.6  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.363 

(Tf  =  2707  K  ;  PAICC  = 
58.7  Atm) 

<X(N7C6H14+)<  < 
X(N(CN)2‘)  < 

0.367 

(Tf  =  2696  K  ;  PAICC  = 
58.8  Atm) 

Paicc  Max 
Composition 

63.1  Atm 

(X(Ca+)=  0.497) 

TF  Max 
Composition 

3807  K 

(X(Ca+)=  0.074) 

Table  37:  Equilibrium  calculations  results  for  the  NjCftH  if"  +  N(CN)f  /  HNO3  system. 


3.3.4  Imidazolium  cations 

As  previously  discussed  in  the  first  part  of  the  report,  different  alkyl  substituted  1- 
methylimidazolium  cation-based  ionic  liquids  have  been  reported  hypergolic  with  white 
fuming  nitric  acid. 

The  first  system  tested  was  the  3-butyl-l-methylimidazolium  cation-based  ionic  liquid. 
Two  different  anions  were  considered:  the  dicyanamide  anion  (see  Figure  18  and  Table  38) 
and  the  dicyanoborate  anion  (see  Table  39). 

For  both  of  these  systems,  the  flame  temperature  is  above  2700  K  for  a  given  domain 
of  compositions.  For  the  ionic  liquid  with  the  dicyanamide  anion,  the  upper  flammability 
limit  is  found  for  a  mole  fraction  of  cation  between  0. 174  and  0. 180.  For  the  ionic  liquid  with 
the  dicyanoborate  anion,  the  UFL  is  found  for  a  mole  fraction  of  cation  between  0.153  and 
0.160.  As  was  seen  before  for  the  N'N-dimethylhydrazinium  cations,  replacement  of  the 
N(CN)7"  anion  by  the  BFflCNfi"  anion  leads  to  a  smaller  flammability  domain. 
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Figure  18:  Flame  temperature  and  P,\iccfor  different  N2CkFI +  N( CN){  /  HNO3  mixtures 

at  298.15  K  and  1  Atm. 


Mixture 

Cation 

HC - £ 

h3g-N-c-N-ch2 

H2  ^CHZ 

6h3 

n2c8h15+ 

Anion 

N(CN)2- 

Oxidant 

HN03 

Stoichiometric 

Equation 

1  [N2C8H15+  N(CN)2-]  +  46/5  HN03 

-*•  8  C02  +  58/5  H20  +  48/5  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(N2C8H15+)  =  X(N(CN)2-) 

X(N2CsH15+)  +  X(N(CN)2‘)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<  X(N2C8H15+)  <  < 
X(N(CN)2-)  < 

9.90E-03 

(Tf  =  2807  K  ;  PAICC  = 
21.5  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.174 

(Tf  =  2718  K  ;  PAICC 
=  37.8  Atm) 

<  X(N2C8H15+)  <  < 
X(N(CN)2-)  < 

0.180 

(Tf  =  2668  K  ;  PAICC  = 
37.7  Atm) 

Paicc  Max 

Composition 

38.6  Atm 

(X(Ca+)=  0.74) 

TF  Max 
Composition 

3632  K 

(X(Ca+)=  0.065) 

Table  38:  Equilibrium  calculations  results  for  the  A7Cs'/7/5+  +  N(CN)f  /  HNO3  system. 
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Mixture 

Cation 

HC - 8 

H^-ch2 

6h3 

N2C8Hi5+ 

Anion 

BH2(CN)2- 

Oxidant 

hno3 

Stoichiometric 

Equation 

1  [N2C8H15+  BH2(CN)2‘]  +  12  hno3 

->  10  C02  +  1/2  B203  +  145/10  H20  +  8  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(N2C8H15+)  =  X(BH2(CN)2-) 

X(N2C8H15+)  +  X(BH2(CN)2‘)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<  X(N2C8H15+)  <  < 
X(BH2(CN)2-)  < 

9.90E-03 

(Tf  =  2820  K  ;  PAICC  = 
21.7  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.153 

(TF  =  2745  K  ;  PAICC 
=  36.7  Atm) 

<X(N2C8H15+)<  < 
X(BH2(CN)2-)  , 

0.160 

(Tf  =  2689  K  ;  PAICC  = 
36.6  Atm) 

Paicc  Max 

Composition 

37.6  Atm 

(X(Ca+)=  0.065) 

Tf  Max 
Composition 

3600  K 

(X(Ca+)=  0.057) 

Table  39:  Equilibrium  calculations  results  for  the  N2CsEli5¥  +  BH2( CN)f  /  HNO3  system. 


Other  alkyl  substituents  were  tested  for  1-methylimidazolium  cation-based  ionic 
liquids.  3 -allyl- 1-methylimidazolium  cation-based  ionic  liquids  were  tested  with  two  different 
anions:  the  dicyanamide  anion  (see  Table  40)  and  the  dicyanoborate  anion  (see  Table  41). 
Both  of  these  ionic  liquids  were  found  hypergolic  experimentally  with  nitric  acid. 

The  3 -propargyl- 1-methylimidazolium  /  nitrocyanamide  ionic  liquid  was  also  tested  in 
this  study  (see  Table  42)  with  nitric  acid. 

The  lower  flammability  limits  of  these  three  ionic  liquid  are  predicted  for  a  mole 
fraction  of  cation  (or  anion)  below  0.0099.  For  the  3 -allyl- 1-methylimidazolium  cation-based 
ionic  liquids,  replacement  of  the  N(CN)2"  anion  by  the  BH2(CN)2"  anion  induces  a  decrease  of 
the  flammability  domain  (see  Table  40  and  Table  41). 

For  the  3 -propargyl- 1-methylimidazolium  /  nitrocyanamide  ionic  liquid,  the  upper 
flammability  limit  is  found  for  a  mole  fraction  of  cation  between  0.342  and  0.346. 


Mixture 

Cation 

Anion 

Oxidant 
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H 

HC=C 

1  1  * 

H,C  C'  CH, 

3  HI2 

HC> 
^CH3 

N2C7Hn+ 

N(CN)2- 

HNO3 

Stoichiometric 

Equation 

1  [N2C7Hii+  N(CN)2"]  +  47/5  HN03 

->  9  C02  +  51/5  H20  +  36/5  N2 

Initial 

Conditions 

P=1  Atm  ;  T  =  298.15  K 

X(N2C7H„+)  =  X(N(CN)2-) 

X(N2C7H„+)  +  X(N(CN)2‘)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(N2C7H„+)<  < 
X(N(CN)2-)  < 

9.90E-03 

(Tf  =  2789  K  ;  PAICC  = 
21.1  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.248 

(Tf  =  2719  K  ;  PAICC 
=  39.3  Atm) 

<X(N2C7Hn  )  <  < 
X(N(CN)2")  < 

0.254 

(Tf  =  2688  K  ;  PA1CC  = 
39.3  Atm) 

Paicc  Max 
Composition 

39.7  Atm 

(X(Ca+)=  0.187) 

Tf  Max 
Composition 

3743  K 

(X(Ca+)=  0.083) 

Table  40:  Equilibrium  calculations  results  for  the  NyCyHn*  +  N(CN)f  /  HNO3  system. 


Mixture 


Cation 


Anion  Oxidant 

BH2(CN)2-  hno3 
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H 

HC=C 

1  1  * 

H,C  C"  CH, 

3  H  1  2 

HC> 

^CHj 

N2C7Hh+ 

Stoichiometric 

Equation 

1  [N2C7H„+  BH2(CN)2-]  +  52/5  HN03 

->  9  C02  +  1/2  B203  +  117/10  H20  +  36/5  N2 

Initial 

Conditions 

P  =1  Atm  ;T  =  298.15  K 

X(N2C7H„+)  =  X(BH2(CN)2) 

X(N2C7Hh+)  +  X(BH2(CN)2-)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(N2C7H„+)<  < 
X(BH2(CN)2‘)  < 

9.90E-03 

(Tf  =  2801  K  ;  PAICC  = 
21.3  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.213 

(Tf  =  2716  K  ;  PAICC 
=  37.5  Atm) 

<X(N2C7H„+)<  < 
X(BH2(CN)2-)  < 

0.219 

(Tf  =  2679  K  ;  PAICC  = 
37.4  Atm) 

Paicc  Max 
Composition 

38.1  Atm 

(X(Ca+)=  0.083) 

TF  Max 
Composition 

3676  K 

(X(Ca+)=  0.074) 

Table  41:  Equilibrium  calculations  results  for  the  N2C7H11'  +  BHyi CN)!  /  HNO3  system. 


Anion 

Oxidant 

Mixture 

Cation 

N(CN)2- 

hno3 
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H 

HC  C 

1  1  ♦ 

H,C  CH. 

3  H  I  2 

■CH 

N2C7H9 

Stoichiometric 

Equation 

1  [N2C7H9+  N(CN)2-]  +  9  HN03  ->  9  CO2  +  9  H20  +  7  N2 

Initial 

Conditions 

P=1  Atm  ;  T  =  298.15  K 

X(N2C7H9+)  =  X(N(CN)2-) 

X(N2C7H9+)  +  X(N(CN)2-)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(N2C7H9+)<  < 
X(N(CN)2-)  < 

9.90E-03 

(Tf  =  2792  K  ;  PA1CC  = 
21.0  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.342 

(Tf  =  2722  K  ;  PAICC 
=  42.7  Atm) 

<X(N2C7H9+)<  < 
X(N(CN)2-)  < 

0.346 

(Tf  =  2682  K  ;  PA1CC  = 
42.1  Atm) 

Pa icc  Max 
Composition 

44.0  Atm 

(X(Ca+)=  0.329) 

Tf  Max 
Composition 

3873  K 

(X(Ca+)=  0.091) 

Table  42:  Equilibrium  calculations  results  for  the  N2C7H9'  +  N(CN)f  /  HNO3  system. 


For  the  1-methylimidazolium  cation-based  ionic  liquids,  different  alkyl  substituents 
were  tested  (butyl,  allyl,  and  propargyl).  For  the  ionic  liquid  with  a  dicyanamide  anion,  the 
impact  on  the  flammability  domain  of  these  alkyl  substituents  in  the  cation  is  as  follows:  butyl 
<  allyl  <  propargyl. 


3.3.5  Pyridinium  cations 

Two  different  pyridinium  cation-based  ionic  liquids  were  studied  in  this  work.  The 
first  one  is  composed  of  the  1-butylpyridinium  cation  and  the  dicyanoborate  anion.  This  ionic 
liquid  has  been  found  experimentally  hypergolic  with  HNO3.  The  results  obtained  from  the 
equilibrium  calculations  are  given  in  Figure  19  and  in  Table  43. 

The  lower  flammability  limit  of  this  system  is  found  for  a  mole  fraction  of  cation  (or 
anion)  below  0.0099.  The  upper  flammability  limit  is  obtained  for  a  mole  fraction  between 
0.153  and  0.160. 

The  maximum  flame  temperature  is  equal  to  3600  K  and  is  obtained  for  a  small 
amount  of  ionic  liquid  (mole  fraction  =  0.057).  The  same  trend  is  predicted  for  the  1- 
allylpyridinium  cation-based  ionic  liquid  (TF  Max  =  3700  K  for  a  mole  fraction  =  0.065)  (see 
Table  44). 

The  ionic  liquid  (NC'sHuf  +  BH2(CN)2")  was  also  determined  hypergolic 
experimentally  with  HNO3.  For  this  system,  the  lower  flammability  limit  was  also  found  for  a 
mole  fraction  of  cation  (or  anion)  below  0.0099.  However,  we  observed  an  increase  of  the 
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upper  flammability  limit.  It  is  predicted  at  a  mole  fraction  of  cation  between  0.213  and  0.219 
(see  Table  44). 


As  seen  previously  for  imidazolium  and  N’N-dimethylhydrazinium  cations  the 
presence  of  the  allyl  substituent  instead  of  the  butyl  substituent  in  the  cation  leads  to  an 
increase  of  the  hypergolicity  domain. 
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Figure  19:  Flame  temperature  and  Paicc  for  different  NC9H iff  +  BH2(CN)2  /  HNO3  mixtures 

at  298.15  K  and  1  Atm. 


Anion 

Oxidant 

Mixture 

Cation 

bh2(cnv 

HN03 
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c— c 

CHa 

nc9h14+ 

Stoichiometric 

Equation 

1  [NC9H14+  BH2(CN)2  ]  +  63/5  HN03 

->  1 1  C02  +  1/2  B203  +  143/10  H20  +  39/5  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(NC9H14+)  =  X(BH2(CN)2-) 

X(NC9H14+)  +  X(BH2(CN)2')  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(NC9H14+)<  < 
X(BH2(CN)2-)  < 

9.90E-03 

(Tf  =  2830  K  ;  PAICC  = 
21.8  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.153 

(Tf  =  2751  K  ;  PAICC 
=  36.2  Atm) 

<X(NC9H14+)<  < 
X(BH2(CN)2-)  < 

0.160 

(Tf  =  2694  K  ;  PAICC  = 
36.0  Atm) 

Paicc  Max 

Composition 

37.5  Atm 

(X(Ca+)=  0.065) 

Tf  Max 
Composition 

3622  K 

(X(Ca+)=  0.057) 

Table  43:  Equilibrium  calculations  results  for  the  NCgHi4+  +  BET( CN)f  /  HNO3  system. 


Mixture 


Cation 


Anion  Oxidant 

BH2(CN)2-  hno3 
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X-=/  C=CH, 

H  2 

nc8h10+ 

Stoichiometric 

Equation 

1  [NC8H10+  BH2(CN)2-]  +  1 1  HNO3 

10  C02  +  1/2  B203  +  115/10  H20  +  7  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(NC8Hk/)  =  X(BH2(CN)2-) 

X(NC8H10+)  +  X(BH2(CN)2-)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(NC8H10+)<  < 
X(BH2(CN)2-)  < 

9.90E-03 

(Tf  =  2812K;PAICC  = 
21.3  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.213 

(Tf  =  2718  K  ;  PAICC 
=  36.7  Atm) 

<X(NC8H10+)<  < 
X(BH2(CN)2-)  < 

0.219 

(Tf  =  2676  K  ;  PAICC  = 
36.6  Atm) 

Paicc  Max 
Composition 

37.8  Atm 

(X(Ca+)=  0.083) 

TF  Max 
Composition 

3700  K 

(X(Ca+)=  0.065) 

Table  44:  Equilibrium  calculations  results  for  the  NCsH]0+  +  BH2( CN)f  /  HNO3  system. 


3.3.6  N-methylpyrrololidinium  cations 

Different  N-methylpyrrololidinium  cation-based  ionic  liquids  are  reported  hypergolic 
experimentally  with  nitric  acid. 

The  first  ones  tested  here  are  1-butyl-N-methylpyrrolidinium  cation-based  ionic 
liquids  with  different  anions  (the  dicyanamide  anion  N(CN)2"  (see  Figure  20  and  Table  45) 
and  the  dicyanoborate  anion  (see  Table  46)).  Equilibrium  calculations  were  carried  out  for 
these  two  ionic  liquids  with  HNO3  as  the  oxidant. 

For  both  systems,  for  a  given  domain  of  compositions,  the  flame  temperature  is  above 
2700  K  and  therefore  the  mixture  is  predicted  hypergolic.  The  lower  flammability  limit  is 
found  for  a  mole  fraction  of  cation  below  0.0099.  For  the  NCgFEo4^  +  N(CN)2"  ionic  liquid, 
the  upper  flammability  is  predicted  for  a  mole  fraction  of  cation  between  0.138  and  0.145.  For 
the  NCylEf/  +  BH2(CN)2-  ionic  liquid,  the  upper  flammability  limit  was  found  for  a  mole 
fraction  between  0.123  and  0.130.  As  has  been  seen  before  for  N-methylpyrrololidinium 
cation-based  ionic  liquids,  replacement  of  the  N(CN)2‘  anion  by  the  BH2(CN)2‘  anion  leads  to 
a  smaller  flammability  domain. 

For  both  of  these  systems,  the  maximum  flame  temperature  was  similar  and  around 
3550  K,  and  was  observed  for  small  mole  fraction  of  cation  (around  0.05  -  0.06)  (see  Table  45 
and  Table  46). 
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Equilibrium  calculations  for  1-allyl-N-methylpyrrolidinium  cation  dicyanoborate 
anion  ionic  liquid  (NCxHA  +  BH^GSTh')  were  also  carried  out.  Like  with  the  previous 
systems,  the  lower  flammability  limit  was  found  for  mole  fraction  of  cation  below  0.0099  (see 
Table  47).  For  this  ionic  liquid,  the  upper  flammability  limit  was  found  for  a  mole  fraction 
between  0.160  and  0.167.  As  seen  previously  for  the  other  cations,  the  presence  of  the  allyl 
substituent  instead  of  the  butyl  substituent  in  the  cation  again  leads  to  an  increase  of  the 
hypergolicity  domain. 
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Figure  20:  Flame  temperature  and  P mcc  for  different  NC9FI20  +  N( CAfA  /  HNO3  mixtures  at 

298.15  K  and  1  Atm. 
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H.C  C^c/C^CH3 

nc9h20+ 

N(CN)2- 

HN03 

Stoichiometric 

Equation 

1  [NC9H20+  N(CN)2-]  +  64/5  HNO3 

->  1 1  C02  +  82/5  H20  +  42/5  N2 

Initial 

Conditions 

P=1  Atm;  T  =  298.15  K 

X(NC9H20+)  =  X(N(CN)2-) 

X(NC9H2o+)  +  X(BH2(CN)2-)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(NC9H2o+)<  < 

X(N(CN)2-)  < 

9.90E-03 

(Tf  =  2832  K  ;  PAICC  = 
22.1  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.138 

(Tf  =  2710  K  ;  PA1CC  = 
37.4  Atm) 

<X(NC9H20+)<  < 
X(N(CN)2‘)  < 

0.145 

(Tf  =  2646  K  ;  PAICC  = 
37.2  Atm) 

Paicc  Max 
Composition 

38.2  Atm 

(X(Ca+)=  0.065) 

Tf  Max 
Composition 

3571  K 

(X(Ca+)=  0.057) 

Table  45:  Equilibrium  calculations  results  for  the  NCgHitf  +  N(CN)f  /  HNO3  system. 


Anion 

Oxidant 

Mixture 

Cation 

BH2(CN)2- 

HN03 
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H,C  C^c/C^CH3 

nc9h20+ 

Stoichiometric 

Equation 

1  [NC9H20+  BH2(CN)2]  +  69/5  HN03 

->  1 1  C02  +  1/2  B203  +  179/10  H20  +  42/5  N2 

Initial 

Conditions 

P=1  Atm;  T  =  298.15  K 

X(NC9H20+)  =  X(BH2(CN)2') 

X(NC9H20+)  +  X(BH2(CN)2')  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(NC9H20+)<  < 
X(BH2(CN)2-)  < 

9.90E-03 

(Tf  =  2845  K  ;  PAICC  = 
22.3  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.123 

(Tf  =  2748  K  ;  PA1CC  = 
36.5  Atm) 

<X(NC9H20+)<  < 
X(BH2(CN)2-)  < 

0.130 

(Tf  =  2677  K  ;  PAICC  = 
36.4  Atm) 

Paicc  Max 
Composition 

37.5  Atm 

(X(Ca+)=  0.057) 

Tf  Max 
Composition 

3559  K 

(X(Ca+)=  0.048) 

Table  46:  Equilibrium  calculations  results  for  the  NC9H20 +  +  BEli(CN)2  /  HNO 3  system. 


Anion 

Oxidant 

Mixture 

Cation 

BH2(CN)2- 

hno3 
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k k" 

O  H 

nc8h16+ 

Stoichiometric 

Equation 

1  [NC8H16+  BH2(CN)2-]  +  61/5  HN03 

10  C02  +  1/2  B203  +  151/10  H20  +  38/5  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(NC8H16+)  =  X(BH2(CN)2-) 

X(NC8H16+)  +  X(BH2(CN)2-)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(NC8H16+)<  < 
X(BH2(CN)2-)  < 

9.90E-03 

(Tf  =  2826  K  ;  PAICC  = 
21.8  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.160 

(Tf  =  2730  K  ;  PAICC 
=  37.2  Atm) 

<X(NC8H16+)<  < 
X(BH2(CN)2-)  < 

0.167 

(Tf  =  2677  K  ;  PAICC  = 
37.1  Atm) 

Paicc  Max 
Composition 

37.8  Atm 

(X(Ca+)=  0.065) 

TF  Max 
Composition 

3626  K 

(X(Ca+)=  0.057) 

Table  47:  Equilibrium  calculations  results  for  the  NCgHiff  +  BET( CN)f  /  HNO3  system. 


3.3.7  N-methyltriazolium  cations 


Two  different  N-methyltriazolium  cation-based  ionic  liquids  have  been  studied  in  this 
work.  The  first  one  is  the  ionic  liquid  composed  of  the  4-butyl- 1 -methyl- 1, 2, 4-triazolium 
cation  and  the  dicyanoborate  anion.  This  ionic  liquid  (N3C7Hi4+  +  BH2(CN)2  )  is  reported 
hypergolic  with  HNO3.  The  results  obtained  from  the  equilibrium  calculations  are  given  in 
Figure  21  and  in  Table  48. 

The  lower  flammability  limit  of  this  system  is  found  for  a  mole  fraction  of  cation  (or 
anion)  below  0.0099.  The  upper  flammability  limit  is  obtained  for  a  mole  fraction  between 
0.180  and  0.187. 

The  maximum  flame  temperature  is  equal  to  3600  K  and  is  obtained  for  a  small 
amount  of  ionic  liquid  (mole  fraction  =  0.065).  The  same  can  be  noted  for  the  4-allyl-l- 
methyl- 1,2, 4-triazolium  cation-based  ionic  liquid  (Tp  Max  =  3700  K  for  a  mole  fraction  = 
0.074)  (see  Table  44). 

This  ionic  liquid  (N3C6Hio+  +  BH2(CN)2')  was  also  determined  hypergolic 
experimentally  with  HNO3.  For  this  system,  the  lower  flammability  limit  was  also  found  for  a 
mole  fraction  of  cation  (or  anion)  below  0.0099.  However,  we  predict  an  increase  of  the  upper 
flammability  limit.  This  time,  it  was  determined  at  a  mole  fraction  of  cation  between  0.254 
and  0.259  (see  Table  49). 
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As  seen  previously  for  the  other  cations,  presence  of  the  allyl  substituent  instead  of 
the  butyl  substituent  in  the  cation  leads  to  an  increase  of  the  hypergolicity  domain. 
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Figure  21:  Flame  temperature  and  Paicc  for  different  N3C7H iff  +  BFFS CNff  /  HNO3 

mixtures  at  298.15  K  and  1  Atm. 
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H 

1  1  «. 

H,C  ch 

3  H  1  z 

H_C^ 
z  CH, 

1 

ch3 

N3C7Hi4+ 

Stoichiometric 

Equation 

1  [N3C7Hi4+  BH2(CN)2]  +  1 1  HN03 

-►  9  C02  +  1/2  B203  +  135/10  H20  +  8  N2 

Initial 

Conditions 

P=1  Atm  ;T  =  298.15  K 

X(N3C7H14+)  =  X(BH2(CN)2) 

X(N3C7H14+)  +  X(BH2(CN)2‘)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<  X(N3C7H14+)  <  < 
X(BH2(CN)2-)  < 

9.90E-03 

(Tf  =  2808  K  ;  PAICC  = 
21.6  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.180 

(Tf  =  2706  K  ;  PAICC 
=  38.0  Atm) 

<  X(N3C7H14+)  <  < 
X(BH2(CN)2-)  < 

0.187 

(Tf  =  2661  K  ;  PAICC  = 
37.9  Atm) 

Paicc  Max 

Composition 

38.2  Atm 

(X(Ca+)=  0.074) 

TF  Max 
Composition 

3620  K 

(X(Ca+)=  0.065) 

Table  48:  Equilibrium  calculations  results  for  the  N3C7H14  +  HEhf  CN)f  /  HNO3  system. 


Mixture 

Cation 

Anion 

Oxidant 
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N  S 

HaC^N''C^:'N''CH2 

3  H  1  2 

hc*ch2 

n3c6h10+ 

BH2(CN)2- 

hno3 

Stoichiometric 

Equation 

1  [N3C6H10+  BH2(CN)2]  +  47/5  HN03 

->  8  C02  +  1/2  B203  +  107/10  H20  +  36/5  N2 

Initial 

Conditions 

P=1  Atm  ;  T  =  298.15  K 

X(N3C6H10+)  =  X(BH2(CN)2-) 

X(N3C6H10+)  +  X(BH2(CN)2-)  +  X(HN03)  =  1 

Lower 

Flammability 
Limit  (LFL) 

0 

<X(N3C6H10+)<  < 
X(BH2(CN)2-)  < 

9.90E-03 

(Tf  =  2789  K  ;  PAICC  = 
21.1  Atm) 

Upper 

Flammability 
Limit  (UFL) 

0.254 

(Tf  =  2719  K  ;  PAICC 
=  39.5  Atm) 

<X(N3C6H10+)<  < 
X(BH2(CN)2-)  < 

0.259 

(Tf  =  2692  K  ;  PAICC  = 
39.4  Atm) 

Paicc  Max 
Composition 

39.5  Atm 

(X(Ca)=  0.225) 

TF  Max 
Composition 

3722  K 

(X(Ca+)=  0.074) 

Table  49:  Equilibrium  calculations  results  for  the  N2Cf,H  uf  +  BH2(CN)2  /  HNO3  system. 
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Conclusions  and  future  work 


This  work  presents  a  numerical  method  devoted  to  the  prediction  of  hypergolic 
mixture  flammability  limits.  A  flame  temperature  criterion  can  be  used  to  detennine 
flammability  limits  of  MMH/NTO  hypergolic  mixtures  and  other  non-ionic  fuel/oxidant 
systems.  In  this  part  were  also  presented  the  ionic  liquids  which  were  found  hypergolic 
experimentally  with  the  oxidants  N2O4  and  HNO3.  The  determination  of  the  thermodynamic 
data  of  cations,  anions  and  neutrals  in  the  gas  phase  was  presented. 

The  results  obtained  from  the  equilibrium  calculations  for  both  protic  and  aprotic  ionic 
liquids  are  presented  in  the  third  section  of  this  part.  For  all  of  the  systems  considered,  we 
were  able  to  determine  flammability  limits.  All  of  the  systems  reported  experimentally  to  be 
hypergolic  are  also  predicted  hypergolic  for  a  given  domain  of  compositions  using  our 
numerical  approach  of  adiabatic  Flame  Temperature  (Tp).  Some  of  these  systems  are 
predicted  hypergolic  for  almost  the  entire  domain  of  compositions.  However,  complementary 
calculations  are  needed  to  assess  these  results. 

Refinements  to  this  approach  are  also  needed  since  the  gas-phase  compositions 
considered  here  for  the  systems  are  somewhat  simplistic.  It  is  tenable  then  to  examine  more 
complex  gas-phase  compositions  in  which  species  formed  in  the  liquid-phase  reactions  are 
also  considered.  Such  species  may  be  isocyanic  acid  (HNCO),  HCN,  NH3,  etc.,  some  of 
which  have  been  experimentally  observed  during  the  pre-ignition  phase. 

New  insight  concerning  the  mixtures  already  studied  in  this  report  may  be  gained  via 
the  detennination  of  thermodynamic  data  for  ion  pairs  for  aprotic  ionic  liquids,  and  neutral 
pairs  for  protic  ionic  liquids.  Additional  ionic  liquid  systems  can  be  also  examined  by 
considering  other  alkyl  substituents  (e.g.  methoxymethyl)  in  the  already  studied  cations 
(imidazolium,  triazolium,  N'N-dimethylhydrazinium,  etc.). 

Finally,  this  thermodynamic  approach  of  predicting  flammability  envelopes  could  be 
fully  developed  into  a  time -resolved  chemical  kinetics  model,  which  implicitly  couples  the 
dynamics  of  the  initial  liquid-phase  chemistry  to  the  subsequent  gas-phase  chemistry  to 
predict  ignition  delay  times  under  realistic  experimental  conditions.  The  last  points  are 
examined  in  parts  two  and  three  of  the  present  report. 
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List  of  Symbols,  Abbreviations  and 

Acronyms 


An  or  An':  Anion 
Ca  or  Ca+:  Cation 
IL(s):  Ionic  Liquid(s) 

LFL:  lower  flammability  limit 
Ne:  Neutral  species 
(s):  singlet  spectral  term 
UFL:  upper  flammability  limit 

WFNA:  White  Fuming  Nitric  Acid  (Pure  nitric  acid:  HNO3) 
X:  Mole  fraction 
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1.  Introduction 


In  the  frame  for  searching  new  hypergolic  (auto-igniting)  propellant  systems,  it  is 
necessary  to  discuss  their  reactions  in  both  the  liquid  and  gas  phases.  In  practical  hypergolic 
applications,  propellant  fuels  and  oxidants  are  liquids  and  upon  contact,  reactions  occur  in  the 
liquid  phase  and  also  in  the  gas  phase  above  the  liquids.  One  point  which  seems  clear  is  that 
ignition  occurs  in  the  gas  phase  even  if  the  reactants  come  into  contact  as  liquids.  The  issue 
of  interest  is  related  to  the  understanding  of  the  qualitative  and  quantitative  compositions  of 
both  phases  when  ignition  occurs.  Several  hypotheses  can  be  proposed.  For  instance,  both 
liquids  vaporize  upon  contact  and  therefore  the  igniting  gas  phase  is  composed  of  the  fuel  and 
oxidant  as  is.  This  is  tenable  for  non-ionic  propellants  because  the  vapor  pressures  of  both 
liquids  are  already  high  as  is  the  case  for  the  gas  phase  ignition  of  monomethyl  hydrazine 
(MMH)  and  nitrogen  tetroxide  (NTO)  [Catoire,  2006].  For  ionic  liquids  (ILs)/oxidant 
combinations,  this  is  also  tenable  but  less  obvious  because  ILs  have  negligible  vapor  pressure 
at  room  temperature.  Nevertheless,  experimental  observations  have  shown  that  isocyanic  acid 
(HNCO)  is  present  in  the  gas  phase  for  the  1  -butyl-3  -methylimidazolium  dicyanamide/white 
fuming  nitric  acid  (BMIM-dca/WFNA)  system  when  ignition  occurs;  presumably  HNCO  is 
formed  during  liquid  phase  reactions  (see  Figure  1),  which  means  that  reactions  in  the  liquid 
phase  together  with  reactions  in  the  gas  phase  are  playing  a  role  during  the  ignition  process. 
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Figure  1.  Interpretation  of  hypergolic  ignition  of  BMIM-dca  with  WFNA.  The  two 
species  fonned  from  the  liquid  phase  intermediate  further  decompose  to  give  gaseous  HNCO, 
which  is  present  when  hypergolic  ignition  occurs  [Cato ire,  2012]. 


For  some  non-ionic  liquid/oxidant  hypergolic  systems  one  elementary  reaction  in  the 
liquid  phase  and  one  elementary  reaction  in  the  gas  phase  have  been  shown  to  be  sufficient  to 
explain  ignition  [Zhang,  2015].  However,  as  underlined  by  the  authors,  their  approach  is  far 
from  being  detailed  and  hence  more  work  is  needed  to  conclude  unambiguously  about  this 
issue.  These  non-ionic  liquid/oxidant  combinations  are  complex  systems  because  chemical 
reaction  dynamics  both  in  the  liquid  phase  and  the  gas  phase  and  their  inter-coupling  need  to 
be  examined.  Furthermore,  in  drop-test  experiments  in  which  a  liquid  droplet  of  fuel  falls  into 
a  volume  of  liquid  oxidant,  the  chemical  dynamics  at  the  liquid  interfaces  also  needs  to  be 
examined.  One  other  issue  of  interest  concerns  the  existence,  or  not,  of  ions  in  the  gas  phase; 
either  fonned  in  the  gas  phase  or  originating  from  the  ionic  liquid  phase  and  their  potential 
role  on  ignition.  All  of  these  features  are  developed  in  this  report  to  allow  for  a  deeper 
understanding  of  the  role  they  play  in  controlling  the  composition  of  the  mixture  formed 
before  ignition  occurs.  This  approach  may  also  allow  for  a  discussion  on  the  ignition  delay, 
which  is  known  to  be  an  important  parameter  in  practical  hypergolic  devices.  Briefly,  if  the 
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ignition  delay  is  too  long,  the  IL/oxidant  combination  is  of  no  interest  in  real  bipropellant 
applications  even  if  its  theoretical  specific  impulse  is  promising. 

2.  Chemical  reaction  dynamics  for  liquid  hydrazines/liquid  NTO  (or  liquid 
HN03)  hypergolic  ignition 


Daimon  et  al.  [Daimon,  1984],  Tanaka  et  al.  [Tanaka,  1985],  and  Daimon  et  al. 
[Daimon,  1991]  observed  a  N2H4  droplet  falling  into  a  NTO  pool  (see  Figure  2)  and  a  NTO 
droplet  falling  into  a  N2H4  pool  (see  Figure  3).  In  both  cases,  they  found  that  a  layer  of  NTO 
was  vaporized  upon  contact.  Their  interpretations  of  the  visualizations  given  in  Figures  2  and 
3  are  reported  in  Figures  4  and  5,  respectively.  In  the  first  case,  ignition  is  supposed  to  be  due 
to  gas-gas  reactions  at  the  interface  of  the  two  liquids  whereas  in  the  second  case,  NTO  vapor 
penetrates  into  the  liquid  hydrazine  and,  although  not  evoked  explicitly  by  the  authors, 
however,  a  gas-liquid  interface  is  formed. 


95 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


(TIME) 


C.— 

(  0 


(  4 1 


(  82 


(124 


(165 


(289  jis) 


(495  ys) 


(412  us) 

D.L. — - 
(435  jjs) 

EX. - 

(454  jjs) 


(247  jjs) 


ENLARGED  PHOTO 


(330  jjs ) 


(371  jjs) 
(392  (is) 


I _ i _ j 

0  5  10mm 


0  10  mm 


C,  :  CONTACT  OF  DROPLET  AND  POOL  SURFACE 

D. L. :  APPEARANCE  OF  DISPERSED  LAYER 
EX.  :  OCCURRENCE  OF  EXPLOSION 


Figure  3.  High  speed  motion  pictures  (48500  frames  per  second)  of  the  contact  of  a 
NTO  droplet  falling  into  a  hydrazine  pool. 
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Figure  5.  Interpretation  of  Tanaka  et  al.  [Tanaka,  1985]  to  explain  hypergolic  ignition 


when  a  droplet  of  NTO  is  falling  into  a  pool  of  liquid  hydrazine. 


For  N2H4  and  its  methylated  derivatives  MMH  and  UDMH  (unsymmetrical  dimethyl 
hydrazine)  and  for  two  oxidizers  NTO  and  FNA  (fuming  nitric  acid)  it  is  reported  that  the 
explosion  phenomena  observed  can  be  classified  into  three  categories  [Daimon,  1991]: 

-  In  the  case  of  hydrazine/NTO,  sudden  gasification  of  a  superheated  liquid  layer 

formed  at  the  boundary  of  the  two  liquids  occurs  spontaneously  and  a  detonation-like 
reaction  proceeds  in  the  reactive  mixture  produced. 

In  the  cases  of  MMH/NTO  and  UDMH/NTO,  sudden  gasification  is  caused  by  the 
shock  from  a  local  ignition  source  and  a  turbulent-combustion  reaction  proceeds  in  the 
reactive  mixture  produced. 

In  the  cases  of  hydrazine  type  fuels/FNA,  sudden  gasification  occurs  spontaneously  as 
in  the  case  of  hydrazine/NTO,  but  it  is  not  augmented  by  a  chemical  reaction  and  in 
these  cases  the  observed  explosion  is  weak. 

In  all  these  cases,  the  vapor  layer,  which  is  formed  between  the  reactive  fuel  droplet 
and  the  liquid  pool,  plays  an  important  role  for  the  occurrence  of  the  explosion.  The 
composition  of  this  layer  (pure  oxidant,  pure  fuel  or  mixture  of  both)  is  not  clearly  discussed 
and  all  these  studies  do  not  identify  unambiguously  the  chemical  cause  of  ignition  or 
explosion.  The  sections  below  are  devoted  to  chemical  dynamics  and  kinetics  in  the  gas 
phase,  in  the  liquid  phase,  at  the  gas-liquid  interface,  and  at  the  liquid-liquid  interface  to 
obtain  a  deeper  understanding  of  the  hypergolic  ignition  phenomenon. 
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3.  Chemical  reaction  dynamics  and  kinetics  for  gas/gas  hypergolic  ignition 


Wasko  [Wasko,  1963],  Chuan  and  Wiber  [Chuan,  1967],  and  Miyajima  and  Sakamoto 
[Miyajima,  1973]  observed  gas  phase  N2H4/NTO  hypergolic  ignitions,  and  Catoire  et  al. 
[Catoire,  2004  &  2006]  studied  gas  phase  MMH/NTO  hypergolic  ignitions  (see  Figure  6)  and 
provided  a  detailed  MMH/NTO  chemical  kinetics  model  to  interpret  the  observed  room- 
temperature  and  low-pressure  self-ignition. 
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Figure  6.  Hypergolic  ignition  observed  when  a  jet  of  gaseous  NTO  enters  into  a 
volume  occupied  by  gaseous  MMH  [Cato ire,  2006], 


A  two-stage  ignition  occurs.  Initially,  a  flame  appears  followed  by  its  extinction. 
After  some  delay,  a  second  flame  appears.  This  two-step  ignition  is  very  interesting  as  it 
suggests  that  something  happens  as  soon  as  the  reactants  come  into  contact  showing  that  the 
interfaces  probably  play  an  important  role  both  in  the  production  of  reactive  species  and  in  the 
production  of  heat.  Very  recently,  Daimon  et  al.  [Daimon,  2014]  provide  a  detailed  chemical 
kinetics  model  for  pure  gas  phase  N2H4/NTO  hypergolic  ignition.  Zhang  et  al.  [Zhang,  2015] 
studied  theoretically  the  reaction  of  azido-N,N-dimethylethanamine  (DMAZ)  and  pure  nitric 
acid.  Two  important  low-temperature  reactions  pathways  were  identified  for  DMAZ/HNO3 
ignition;  one  in  the  liquid  phase  and  one  in  the  gas  phase.  In  the  liquid  phase,  it  was  proton 
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transfer  from  the  nitric  acid  to  the  amine  nitrogen  of  DMAZ.  In  the  gas  phase,  it  was  the 
reaction  between  DMAZ  and  NO2.  Similarly,  Catoire  et  al.  [Catoire,  2012]  have  discussed 
the  role  of  isocyanic  acid  formed  in  the  liquid  phase  on  the  kinetics  of  gas  phase  ignition. 

To  conclude,  it  appears  that  processes  at  interfaces  may  play  a  major  role  in  sustaining 
hypergolic  ignition.  Both  interface  types,  gas-gas  and  gas-liquid,  need  to  be  considered.  For 
non-ionic  liquid/liquid  oxidant  systems,  reactions  in  the  gas  phase  are  probably  of  major 
importance,  whereas  for  IL/liquid  oxidant  systems,  the  issue  of  which  phases  are  important 
remains  to  be  worked  out  and  is  further  examined  herein. 

4.  Chemical  reaction  dynamics  in  the  gas  phase 


The  types  of  reactions  in  the  gas  phase  are  generally  well  known  and  are  not 
exhaustively  presented  here.  Generally,  only  neutral  species  (molecules,  radicals  and  atoms) 
are  involved.  Even  if  ions  exist  in  the  gas  phase,  they  are  much  less  numerous  than  neutral 
species  and  their  kinetics  role  is  generally  negligible  and  therefore  neglected  in  non-ionic 
liquid/oxidant  systems.  For  NTO  and  WFNA  as  oxidants,  NCE  molecules  react  with  gas 
phase  molecules  of  the  neat  fuel  or  with  product  species  either  fonned  from  liquid  phase 
reactions  or  from  the  initial  reaction  between  the  fuel  and  NO2.  For  ILs/oxidant  systems,  the 
existence  of  ions  however  must  be  discussed.  If  one  assumes  that  the  vapor  phase  from  ILs 
plays  a  role,  four  possibilities  exist.  Above  the  liquid  phase,  ion  pairs  can  exist,  or  separated 
ions  can  exist,  or  neutral  pairs  can  exist,  or  finally  separated  neutrals  can  exist.  Experimental 
findings  and  theoretical  calculations  are  reported  in  the  literature  on  this  issue.  It  appears  that 
ionic  liquids  have  to  be  divided  into  two  classes;  protic  ionic  liquids  and  aprotic  ionic  liquids. 
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Protic  Ionic  Liquids 


According  to  Earle  and  co-workers,  for  protic  ionic  liquids,  an  equilibrium  exists 
between  ionic  and  non-ionic  species  [Earle,  2006]: 

BH+(t)  +  X  (£)  <-*■  B(()  +  HX({)  <-*■  B(g)  +  HX(g) 

where  (£)  represents  the  liquid  phase  and  (g)  represents  the  gas  phase.  Calculations  on 
protic  triazolium-based  energetic  ionic  liquids  show  that  the  neutral  pair  is  of  lower  energy 
than  the  ionic  pair  [Schmidt,  2005],  Separation  of  the  species  in  the  ion  pair  requires  a 
significant  amount  of  energy  and  therefore  the  ionic  pair  once  formed  will  transfer  one  proton 
to  form  a  neutral  pair.  This  proton  transfer  is  therefore  fundamental  for  the  decomposition  of 
the  triazolium-based  ILs  considered.  However,  Baranyai  et  al.  did  not  observe  the  fonnation 
of  species  according  to  the  proton  transfer  process,  and  therefore,  they  questioned  this 
mechanism  [Baranyai,  2004],  On  the  other  hand,  Zorn  et  al.  deduced  from  their  calculations 
on  protic  tetrazolium-based  ionic  liquids  that  when  a  cation  is  paired  with  an  oxygen  rich 
anion,  an  isolated  gas  phase  ion  pair  was  not  generally  found  to  be  stable  [Zorn,  2006].  In 
fact,  a  proton  transfers  barrierlessly  from  the  cation  to  the  anion  to  form  a  neutral  pair.  This 
result  is  also  consistent  with  the  experimental  findings  of  Leal  and  co-workers,  who  reported 
that  the  gas  phase  above  a  protic  ionic  liquid  is  composed  of  separated  neutral  molecules 
[Leal,  2007], 


Aprotic  Ionic  Liquids 

Calculations  carried  out  by  EmeTyanenko  et  al.  on  the  gas  phase  equilibrium: 

ion  pair  <-»•  separated  ions 

for  1  -butyl-3  -methylimidazolium  dicyanamide  indicate  that  dissociation  of  the  ion  pair 
is  marginal  [EmeTyanenko,  2007].  Experimental  and  theoretical  works  reported  by  Strasser 
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et  al.  suggest  that  the  vapor  pressure  of  an  aprotic  ionic  liquid  (l-ethyl-3-methylimidazolium 
bis(trifluoromethylsulfonyl)imide)  is  composed  of  isolated  ion  pairs  [Strasser,  2007].  The 
formation  of  clusters  of  ion  pairs  was  reported  to  be  negligible.  This  conclusion  is  supported 
by  another  experimental  work  perfonned  by  Leal  et  al.,  who  found  that  the  gas  phase  is 
composed  of  only  isolated  ion  pairs  and  that  their  experimental  observations  did  not  support 
the  formation  of  ion  clusters  [Leal,  2007].  For  aprotic  ionic  liquids,  the  only  possible 
mechanism  for  volatilization,  according  to  Earle  et  al.,  is  as  ion  pairs,  isolated  or  as  clusters 
[Earle,  2006].  Armstrong  et  al.  performed  an  analysis  of  the  vapor  pressure  of  l-ethyl-3- 
methylimidazolium  ethylsulfate  by  mass  spectrometry  [Armstrong,  2007].  They  deduced 
from  their  experiments  that  the  vapor  pressure  of  these  ionic  liquids  was  due  to  ion  pairs, 
although  they  did  not  observe  directly  these  ion  pairs  because  the  ionized  ion  pair  (fonned  in 
the  ionization  chamber  of  the  mass  spectrometer  by  electron  bombardment)  was  unstable  and, 
therefore,  fragmentation  occurred  in  their  experiments.  One  other  finding  by  Armstrong  et  al. 
is  that  clusters  of  ion  pairs  do  not  seem  to  form  in  the  gas  phase  [Armstrong,  2007]. 

5.  Gas  diffusion  in  bulk  liquids 


Condensed  phase  concentration  of  a  gaseous  species  is  limited  by  its  solubility  in  the 
liquid  of  interest.  At  equilibrium,  the  Henry’s  law,  kh  =  c,/pL  links  partial  pressure  pj  of  a 
species  i  to  its  liquid  phase  concentration  Ci.  Henry’s  constants  kh  are  available  for  many 
gases  (H2,  He,  Ne,  etc.)  in  many  solvents  (alkanes,  alcohols,  etc.).  However,  Henry’s  law 
describes  the  steady  state  but  for  real  reacting  systems  it  is  doubtful  that  this  equilibrium  is 
reached  due  to  diffusion  limitations.  The  diffusion  inside  the  bulk  is  linked  to  the  interfacial 
area  which  can  be  “natural”  (planar  interface  between  gas  and  liquid)  or  “forced”  (as  it  is 
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probably  the  case  when  a  droplet  free  falls  into  a  pool).  In  both  cases,  the  diffusion  process 
can  be  expressed  using  a  film  model  such  as  the  Whitman  model  presented  in  Figure  7  for  a 
reacting  system.  Furthermore,  a  competition  between  diffusion  into  the  bulk  and 
consumption  in  the  interface  may  lead  to  smaller  concentrations  that  are  finally  reached  in  the 
bulk.  Thus  it  could  well  be  a  fact  that  the  reaction  phenomenon  in  the  bulk  is  somehow 
negligible  before  ignition,  i.e.,  everything  happens  in  the  interface  and  this  finally  leads  to 
ignition. 


Interface 

I 


Figure  7.  Whitman  model  for  a  reacting  system  consuming  the  species  A  in  the  liquid 
phase.  Pa  is  the  partial  pressure  of  species  A.  g,  l  and  i  stand  for  gas,  liquid  and  interface, 
respectively. 

6.  Chemical  reaction  dynamics  in  solution  (in  liquids) 


For  the  problem  at  hand,  it  is  essential  to  have  an  understanding  of  the  nature  of 
chemical  reactions  in  solution.  Much  less  is  known  about  the  reaction  dynamics  in  liquids  at 
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the  microscopic  level  than  about  the  dynamics  of  gas  phase  reactions.  This  statement 
[Schroeder,  1987]  was  valid  30  years  ago  and  is  most  certainly  valid  today  [Elies,  2006],  and 
also  because  in  recent  decades,  ab  initio  computations  of  gas  phase  thermochemistry  and  rate 
coefficients  have  become  routine  activities.  Some  concepts  derived  from  gas  phase  studies 
are  also  valid  for  reactions  in  solution,  but  on  the  other  hand,  many  characteristics  of  the 
dynamics  in  liquids  are  not  important  for  the  gas  phase  and  therefore  the  liquid  does  not 
behave  like  a  gas  of  high  density  [Schroeder,  1987].  For  instance,  solvent  molecules  can 
affect  the  motion  of  reactants  in  a  minor  to  a  profound  manner  [Elies,  2006],  The  effects  of  a 
solvent  on  a  reaction  rate  may  be  either  chemical  or  physical.  If  the  solvent  molecules 
participate  directly  in  the  reaction  mechanism,  the  effect  is  chemical.  If  the  role  of  the  solvent 
is  merely  to  modify  the  interaction  between  reacting  species,  its  effect  may  be  termed 
physical.  A  reaction  between  two  molecules  in  a  solution/liquid  can  be  thought  as  occurring 
in  three  well  defined  steps: 

(1)  diffusion  of  the  molecules  to  each  other 

(2)  the  actual  chemical  transformation 

(3)  diffusion  of  the  products  away  from  each  other 

The  first  way  in  which  a  solvent  medium  can  affect  a  reaction  rate  is  via  the 
modification  of  the  activation  energy  by  stabilization  of  the  transition  state  as  illustrated  in 
Figure  8. 
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Transition  state 


Figure  8.  Energy  diagram  for  a  bimolecular  reaction  in  the  gas  phase  and  in  solution 
[Elies,  2006]. 


If  activation  energies  are  impacted,  it  means  that  both  thermochemistry  and  kinetics 
are  impacted  as  shown  in  Figure  9.  Solvents  can  also  play  a  role  by  initiating  reactions, 
stabilizing  intermediates  (both  ionic  and  radical),  and  moderating  the  collisions  between 
molecules.  In  fact,  the  collision  frequency  is  lower  in  solution  than  in  the  gas  phase  because 
solvent  molecules  get  in  the  way,  but  the  net  effect  is  that  the  presence  of  a  solvent  often 
increases  the  rate  of  a  reaction,  especially  the  reaction  rates  of  ions.  Generally,  solvents  have 
only  a  small  effect  on  reactions  of  neutral  species. 
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Figure  9.  Energy  changes  for  the  reaction  Cf  +  CH3Br  — >CH3C1  +  Br'  in  gas  phase 
and  in  liquid  phase  (DMF).  CP  and  TS  stand  for  complex  and  transition  state,  respectively. 
Data  are  from  Reichardt  and  Welton  [Reichardt,  2011]. 


The  impact  of  solvation  on  thermochemistry  and  kinetics  has  been  discussed  by  Jalan 
et  al.  [Jalan,  2010]  among  others.  To  treat  these  issues,  two  approaches  are  possible.  The 
first  one  deals  with  experimental  measurements,  which  continues  to  be  the  preferred 
approach.  The  second  one  consists  of  applying  solvent  corrections  to  the  gas  phase  data.  This 
treatment  is  presented  herein  in  sections  (A)  and  (B)  below. 

(A)  Solvent  effects:  impact  on  the  thermochemistry 

Corrections  are  illustrated  by  the  thermodynamic  cycle  given  in  Figure  10.  It  includes 
individual  corrections  applied  from  a  real  gas  mixture  to  a  concentrated  solute  in  solvent. 
Corrections  can  be  divided  into  four  main  steps  [Orozco,  2000]: 
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(1)  ideal  state  corrections  in  which  gases  and  liquids  are  “moved”  from  concentrated 


to  diluted  conditions 


(2)  cavity  fonnation  corrections 


(3)  Van  der  Waals  corrections  (also  called  dispersion-repulsion  forces)  and 


(4)  electrostatic  contributions,  which  refer  to  the  charge  distribution  in  the  solvent. 
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Figure  10.  Thermodynamic  cycle  representing  a  solvation  process.  Grey  arrows 
represent  what  is  commonly  computed;  the  free  energy  change  between  ideal  gas  to  diluted 
solute  in  solvent.  Diagram  from  Jalan  et  al.  [Jalan,  2010]. 
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The  accuracy  of  the  results  obtained  is  indeed  a  function  of  each  correction,  which  can 
be  evaluated  separately  by  using  molecular  modeling  or  corrected  collectively  by  using  semi- 
empirical  methods.  The  accuracy  of  the  methods  depends  on  the  description  of  the  solute- 
solvent  interactions  in  both  semi-empirical  methods  with  solute-solvent  descriptors  and  in 
computational  modeling  methods  using  different  solvent  description  theories.  The  level  of 
description  of  solvation  in  molecular  modeling  is  presented  in  Figure  1 1 .  These  solvation 
models  can  be  divided  into  two  approaches.  In  the  first  one,  solvent  molecules  surrounding 
solutes  are  represented  explicitly  (see  Figures  1 1  a,  b,  and  c).  In  the  second  one,  a  continuum 
model  is  used  to  model  all  (Figure  1  If)  or  part  of  the  solvent  (Figures  1  Id  and  e). 


Figure  11.  Solvent  model  illustration  for  discrete  (a)  to  continuum  models  (f)  [Jalan, 
2010].  Grey  represents  continuum  dielectric. 


In  these  two  approaches,  molecules  can  be  either  described  using  quantum  mechanics 
(QM)  or  molecular  mechanics  (MM).  Evaluations  of  these  methods  and  their  accuracies  are 
reported  in  Jalan  et  al.  [Jalan,  2010].  As  is  usual  in  modeling  studies,  accuracy  of  results  and 
computational  costs  need  to  be  balanced.  The  number  of  electrons  in  the  molecules  that  are 
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included  in  the  model  is  the  main  limiting  factor.  For  this  reason,  semi-empirical  methods, 
like  the  LSER  (linear  solvation  energy  relation)  method  described  hereafter,  are  still  widely 
used.  These  methods  require  fewer  parameters  that  are  based  on  the  physical  properties  of 
both  the  solute  and  solvent  to  obtain  solvation  corrections.  Many  solute/solvent  pairs  can  be 
investigated  and  the  solvation  correction  is  obtained  using  a  small  set  of  equations. 

However,  only  stable  species  parameters  can  be  obtained  experimentally  and 
calculations  are  still  required  for  radicals  and  other  short  lifetime  species.  The  liquid  phase 
thermochemistry  can  be  calculated  in  two  steps  using  the  LSER  method.  The  first  one  (Eq. 
5),  requires  calculating  Guq,  the  liquid  phase  free  energy  of  the  species  estimated  by  applying 
free  energy  of  solvation  correction  AGsoiv,  to  its  gas  phase  free  energy  Ggas.  AGsoiv  is  obtained 
by  using  Eq.  6  and  the  partition  coefficient  Kvapor-soivent  of  a  given  solute/solvent  pair.  The 
LSER  method  based  on  Abraham  et  al.’s  models  (Eq.  7)  [Abraham,  2010]  provides  an 
estimation  of  this  partition  coefficient  at  298  K.  In  Eq.  7,  the  first  set  is  referring  to  solute 
descriptors  (E,  S,  A,  B,  and  L)  and  the  second  to  solvent  descriptors  (c,  e,  s,  a,  b,  and  I).  The 
next  step  is  to  calculate  the  enthalpy  of  solvation  AHsoiv  at  298  K,  using  Mintz’s  model  also 
based  on  Abraham  et  al.’s  model  just  described  (Eq.  8).  With  Mintz’s  model,  the  solute 
descriptors  remain  the  same  and  the  new  solvent  descriptors  (c’,  a’,  b’,  e’,  s’,  and  I’)  are 
obtained  through  another  empirical  correlation.  It  is  then  possible  to  obtain  ASsoiv  according 
to  AG  =  AH  -  TAS.  AGsoiv  can  be  expressed  as  a  function  of  temperature  following  Eq.  9, 
wherein  AHsoiv  and  ASsoiv  are  assumed  to  be  temperature  independent.  These  assumptions 
provide  good  agreement  (with  root  mean  square  error  of  0.5  kcal/mol)  for  many  systems  as 
presented  by  Jalan  et  al.  [Jalan,  2013]  with  over  900  data  points. 
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Eq.  5 
Eq.  6 
Eq.  7 
Eq.  8 
Eq.  9 


In  Eqs.  7  and  8,  the  sS  and  s’S  terms  correspond  to  the  electrostatic  interactions,  the 
eE  and  e’E  account  for  residual  contributions  from  dipolarity/polarizability,  the  IL  and  I’L 
terms  are  related  to  the  cavity  formation  and  the  dispersion  interactions,  aA,  a’A  and  bB,  b’B 
terms  stand  for  the  ability  to  donate  protons  and  accept  protons,  respectively. 

(B)  Solvent  effects:  impact  on  the  kinetics 

Solvent  effects  may  impact  transition  state  theory  (TST)  calculations  according  to 
three  different  ways: 

(1)  by  increasing  or  decreasing  the  activation  energy  Ea  as  presented  in  Figure  8. 
Many  observations  of  kinetics  solvent  effects  are  presented  in  the  literature  where  the  rate 
constants  may  be  increased  or  decreased  compared  to  their  gas  phase  values.  Kinetics  solvent 
effect  is  important  for  Menschutkin  reactions  as  shown  in  Figure  12.  It  is  relevant  to  reactions 
where  ions  are  formed  from  neutral  species.  When  the  transition  state  (TS)  is  more  polar  than 
the  reactants,  the  use  of  polar  solvents  stabilizes  the  TS  and  decreases  the  value  of  the 
activation  energy  Ea. 
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Figure  12.  Solvent  effect  for  Menschutkin  reactions.  A  polar  solvent  stabilizes  the 


transition  state. 


(2)  by  solvent  friction. 

As  presented  by  Hynes  [Hynes,  1985],  solvents  may  impact  the  TS  due  to  solvent 
friction.  In  case  of  competing  reactions,  any  one  reaction  that  is  less  impacted  by  solvent 
friction  may  become  predominant  compared  to  the  situation  observed  in  the  gas  phase. 


(3)  by  local  (at  cage  scale)  diffusion-reaction  competition. 

Finally,  solvents  may  greatly  affect  the  mass  transfer  rates  during  reactions.  The  key 
feature  is  that  the  solvent  slows  down  the  diffusion  rate  of  the  reactants,  which  affects  the 
collision  rate.  In  addition,  the  solvent  tends  to  trap  the  reactants  in  a  solvent  cage,  which  has  a 
dramatic  effect  on  the  rate.  Qualitative  predictions  of  the  effects  of  solvents  on  the  rate  are 
possible  through  Hughes-Ingold  rules  for  ionic  reactions  [Masel,  2001]: 
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•  if  the  TS  for  a  reaction  has  a  larger  charge  than  the  reactants  do,  then  the  rate  of 


reaction  will  increase  as  the  polarity  of  the  solvent  increases 

•  if  the  TS  for  a  reaction  has  a  smaller  charge  than  the  reactants,  then  the  rate  of  reaction 
will  decrease  as  the  polarity  of  the  solvent  increases 

•  if  the  net  charge  remains  the  same,  but  the  charge  is  dispersed,  then  there  will  be  a 
small  decrease  in  rate  as  the  polarity  of  the  solvent  increases 

•  if  the  net  charge  remains  the  same,  but  the  charge  is  localized,  then  there  will  be  a 
small  increase  in  rate  as  the  polarity  of  the  solvent  increases 

Semi-quantitative  predictions  are  also  possible  through  four  key  semi-quantitative 
models: 

•  the  double-sphere  model 

•  the  single-sphere  model 

•  Debye-Huckel  theory 

•  Regular  solution  theory 

Essential  connections  between  chemical  dynamics  in  gases  and  liquids  are  described 
by  Elies  and  Crim  [Elies,  2006].  Existing  solution  reaction  types  are:  (a)  radical 
recombinations,  (b)  precipitation  reactions,  (c)  acid-base  reactions,  (d)  electron  transfer 
reactions,  (e)  oxidation-reduction  or  redox  reactions,  (f)  inorganic  exchange  and  substitution 
reactions,  (g)  organic  exchange  and  substitution  reactions,  (h)  elimination  and  decomposition 
reactions,  (i)  atom  transfer  reactions,  (j)  isomerization  reactions,  (k)  charge  transfers  and 
isomerizations,  and  (1)  tunneling  reactions,  which  are  further  described  below. 
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(a)  Radical  recombination  reactions: 


Radical  recombination  reactions  are  simple  because  they  are  generally  barrierless.  One 
example  of  such  reactions  is  the  iodine  atom  recombination: 

I  +  I  — >  I2. 


(b)  Precipitation  reactions: 

Precipitation  reactions  are  reactions  in  which  a  precipitate  is  formed.  All  potential 
precipitates  are  considered  by  writing  all  possible  combinations  between  cations  and  anions 
present  in  the  liquids.  Then  solubility  rules  are  used  to  detennine  which  combinations  are 
insoluble  and  able  to  precipitate. 


(c)  Acid-base  reactions: 

H+  +  NH3  -»•  NH4+ 

(d)  Electron  transfer  reactions: 

Fe2+  +  Fe  3+  — >  Fe3+  +  Fe*2+  (isotopic  exchange  electron  transfer  reactions)  where  the 
asterisk  denotes  a  radioactive  isotope  [Marcus,  1993]. 


(e)  Oxidation-reduction  or  redox  reactions: 

Fe2+  +  Ce4+  — >  Fe3+  +  Ce3+  (electron  transfer  reactions  between  two  different  redox 
systems  [Marcus,  1993]) 
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(f)  Inorganic  exchange  and  substitution  reactions: 


[Co(H20)6]2+  +  4  Cl  ->  [CoCl4]2"  +  6  H20 

(g)  Organic  exchange  and  substitution  reactions: 

These  are  the  so  called  SnI  and  Sn2  reactions  or  nucleophilic  substitution  reations 

such  as 

RBr  +  OH'  ->  ROH  +  Br¬ 
ill)  Elimination  and  decomposition  reactions: 

(i)  Atom  transfer  reactions: 

Atom  transfer  reactions  such  as  AB  +  C  — >  A  +  BC.  It  is  the  transfer  of  an  atom  or  a 
group  B  from  AB  to  form  BC  [Marcus,  1993]. 

HI  +  Br  — >  I  +  HBr 

(j)  Isomerization  reactions: 

(k)  Charge  transfers  and  isomerizations: 

(l)  Tunneling  reactions: 

7.  Reactions  in  solution  (in  liquids)  of  excited  states 
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Thomas  provides  a  review  on  the  reactions  of  excited  states  in  liquids  [Thomas,  1970]. 
The  occurrence  of  excited  states  is  not  assessed  in  the  liquid  phase  but  excited  OH  radicals  are 
certainly  fonned  in  the  gas  phase  when  H2O2  is  used  as  an  oxidant  as  well  as  excited  NO2 
radicals  are  probably  fonned  when  the  oxidant  is  NTO  or  WFNA.  However,  these  reactions 
are  probably  marginal  because  usually  species  exist  in  their  ground  state.  Excited  state  ion 
reactions  are  also  observed  such  as: 

0+(2D)  +  n2  -+  n2+  +  o 


8.  Ion-molecule  reactions 


Finally,  other  ions  (positive  or  negative)  can  also  be  fonned  upon  heating  due  to 
exothermic  reactions,  independently  from  the  ions  initially  present  in  the  ILs.  Various  types  of 
reactions  are  possible  in  the  family  of  ion-molecule  reactions  in  the  gas  phase  [Ferguson, 
1975].  Either  chemically  stable  neutral  gases  or  “unstable”  neutrals  such  as  atoms  or  free 
radicals  in  their  ground  state  or  excited  species  can  participate.  For  instance,  charge  transfer 
reactions  such  as: 

N2"1"  +  O  — »  N2  +  0+ 

generally  compete  with  ion-atom  interchange  reaction  such  as: 

N2+  +  O  ->  NO+  +  N 

or  in  associative-detachment  reactions,  where  the  energy  to  detach  an  electron  is 
provided  by  the  formation  of  a  new  bond  between  the  reactants,  may  be  fast  such  as  in: 

H'  +  H  — >■  H2  +  e' 

OH'  +  H  — >■  H20  +  e' 
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and  charge  transfer  reactions  in  doubly  charged  ions  may  be  important  such  as  in: 

Ca2+  +  NO  — >  Ca+  +  NO+ 

However,  most  ionic  reactions  occur  in  solutions  or  on  solid  catalysts.  The  major 
difference  between  a  reaction  in  the  gas  phase  and  in  solution  is  that  the  ionic  species  are 
stabilized  by  the  presence  of  the  solution  resulting  in  ions  that  are  much  more  stable  than 
those  in  the  gas  phase.  The  ion  is  generally  surrounded  by  a  solvent  cage.  In  order  for  a 
reaction  to  occur,  one  needs  to  first  distort  the  ion  cage  so  the  reactants  can  get  close  enough 
to  react.  Because  the  exact  nature  of  the  distortion  is  unclear,  it  is  at  present,  very  difficult  to 
accurately  predict  how  much  the  value  of  the  intrinsic  barrier  changes  as  a  result  of  the 
solvent.  Fortunately,  the  change  in  the  intrinsic  barrier  does  not  really  change  the  mechanism 
of  a  reaction.  All  of  the  same  elementary  reactions  observed  in  the  gas  phase  are  also  seen  in 
solution  phase.  Olah  et  al.  [Olah,  1997]  showed  that  one  can  often  predict  mechanisms  of 
reaction  without  considering  the  changes  in  the  intrinsic  barriers  because  of  the  solvent.  It  is 
just  that  the  relative  rates  of  reactions  change  in  solution.  Often,  the  intrinsic  barriers  to  all  of 
the  reactions  of  a  species  are  changed  by  a  similar  amount. 

9.  Chemical  reactions  at  gas-liquid  interfaces 


As  underlined  by  Benjamin,  [Benjamin,  2015]  chemical  reactions  that  take  place  at  the 
interface  between  a  liquid  and  a  second  phase  (vapor,  another  immiscible  liquid  or  a  solid)  are 
of  major  interest.  Ordinarily,  the  gas-liquid  interface  is  generally  considered  to  play  the  role 
of  a  “gateway”  into  the  liquid,  and  is  believed  not  to  play  the  role  of  a  reaction  medium. 
However,  it  is  shown  in  this  report  that  the  gas-liquid  interface  is  a  distinct  region  where  often 
unique  reactions  can  occur.  Since  Daimon  et  al.  and  Tanaka  et  al.  have  previously  presented 
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evidence  that  a  gas-liquid  interface  exists  just  before  ignition,  it  may  therefore  be  possible  that 
reactions  taking  place  in  this  region  are  the  cause  of  ignition.  It  is  clear  that  during  the 
ignition  delay,  reactions  occur  in  the  liquid  phase  and  that  reactions  occur  in  the  gas  phase, 
but  the  gas-liquid  interface  aspects  have  not  been  discussed  at  all  in  any  of  the  previous 
papers  on  hypergolic  ignition.  It  is  therefore  of  interest  to  examine  the  processes  at  the  gas- 
liquid  interface,  especially  as  to  how  gas  phase  species  enter  into  the  liquid  phase  and  most 
importantly  how  they  react.  The  contrary,  the  entering  into  the  gas  phase  of  a  species  formed 
in  the  liquid  phase,  is  simple  to  understand;  if  the  species  formed  is  gaseous  or  has  a  high 
vapor  pressure,  it  will  be  liberated  as  is  into  the  gas  phase. 


a.  Collision  of  a  gas  phase  molecule  with  a  liquid  surface 


Nathanson  [Nathanson,  2004]  has  presented  a  review  of  experiments  perfonned  to 
clarify  the  role  of  the  surface  region  during  interfacial  and  bulk  phase  reactions.  When  a  gas 
phase  molecule  collides  with  a  liquid  surface  it  can  be  sent  back  to  the  gas  phase  through 
inelastic  scattering  or  it  can  be  trapped,  i.e.,  its  kinetic  energy  is  absorbed  and  it  is  able  to  bind 
to  the  surface.  The  probability  of  the  trapping  event  depends  on  many  parameters:  the  kinetic 
energy  of  the  impinging  molecule,  the  mass  of  this  molecule  and  the  ones  at  the  surface,  the 
liquid  temperature,  its  surface  tension,  etc.  Experiments  have  demonstrated  that  gas-liquid 
reactions  can  occur  as  soon  as  the  gas  phase  molecule  enters  the  interfacial  region  and  this 
region  is  therefore  a  reaction  medium. 


b.  Collision  of  a  gas  phase  strong  acid  with  a  ionic  liquid  surface 
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Figure  13  illustrates  the  pathways  observed  after  a  gaseous  acid,  HX  collides  with  the 
surface  of  a  solution  of  deuterated  sulfuric  acid,  D2SO4  in  deuterated  water,  D20. 


HX  HX  HX  DX 


Exchange 


Further  Reaction 


Figure  13.  Observed  pathways  for  an  acidic  molecule  HX  colliding  into  D2O/D2SO4 
[Nathanson,  2004]. 


The  HX  molecules  are  sent  back  to  the  gas  phase,  or  some  fraction  are  trapped  and 
then  desorb  unchanged.  The  remaining  molecules  react  in  the  acid  interfacial  region  and/or  in 
the  acid  bulk  as  intact  HX  or  as  H+  and  X'  after  dissociation.  These  processes  finally  lead  to 
the  formation  of  DX  molecules  in  the  gas  phase  once  desorbed  from  the  liquid  surface.  If  X  is 
Cl,  then  the  gas  phase  acid  is  HC1  (see  Figure  14).  Further  reactions  cited  in  Figure  13  are 
unspecified. 

120 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


Figure  14.  Observed  pathways  for  HC1  colliding  into  D2O/D2SO4  [Nathanson,  2004], 


HC1  enters  deuterated  sulfuric  acid  solution  in  two  steps:  a  fraction  of  the  impinging 
molecules  first  thermalize  upon  collision  with  the  surface  and  then  a  fraction  of  these  trapped 
molecules  enter  the  acid  as  intact  HC1  or  as  H+  and  Cl  .  Experiments  show  that  at  low 
incident  energy,  nearly  all  HC1  molecules  are  trapped,  but  the  amount  of  DC1  formed 
indicates  that  only  1  in  9  HC1  molecules  dissociates  in  the  interfacial  or  bulk  regions  of  the 
acid  and  that  the  remaining  8  out  of  9  trapped  HC1  molecules  desorb  from  the  acid  before 
reacting.  These  results  are  particularly  interesting  for  ILs/acid  hypergolic  systems  if  one 
considers  that  the  impinging  molecule  is  WFNA  (a  strong  acid  as  is  HC1)  at  the  surface  of  an 
ionic  liquid  (as  is  the  D2O/D2SO4  solution).  Most  likely,  HNO3  will  dissociate  into  H+  and 
NO3'  ions,  which  then  react  with  the  ions  from  the  ionic  liquid. 
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c.  Collision  of  a  gas  phase  strong  acid  with  a  neutral  liquid  surface 


The  interfacial  and  bulk  reactions  of  a  gas  phase  acid  with  a  neutral  liquid  is  also  of 
interest  for  hypergolic  ignition  as  is  the  case  with  the  non- ionic  liquid,  MMH  and  HNO3 
system.  Gas  phase  hydrogen  chloride  HC1  and  deuterated  hydrogen  chloride  DC1  have  been 
studied  with  glycerol,  a  neutral  liquid.  Experiments  show  that  most  of  the  DC1  molecules 
dissolve  into  glycerol.  Of  the  DC1  molecules  that  adsorb  (or  thermalize),  20%  are  found  to 
desorb  before  dissolution  and  7%  of  DC1  are  converted  to  HC1  in  the  interfacial  region  and 
desorb  before  they  can  be  dissolved  into  the  bulk.  It  is  presumed  that  the  HC1  molecules 
observed  in  the  gas  phase  during  DCl/liquid  phase  glycerol  reaction  are  formed  according  to 
an  ionic  pathway  described  in  Figure  15. 
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Figure  15.  Scattering  pathways  for  DCl  collisions  with  glycerol.  Only  the  partitioning 
between  direct  scattering  and  trapping  depends  on  Einc  and  tinc.  7%  of  the  trapped  DCl 
undergoes  rapid,  near-interfacial  exchange  followed  by  immediate  HC1  desorption 
[Nathanson,  2004]. 


Another  interesting  feature  is  that  HC1  molecules  dissolve  exothermically  in  glycerol. 
Although  this  is  neither  surprising  nor  new  (the  same  holds  for  the  dissolution  of  HC1  in 
water),  it  is  however  worthy  to  note  because  this  will  lead  to  an  increase  in  the  temperature  of 
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the  interfacial  region,  which  can  potentially  play  a  beneficial  role  in  the  energy  transfer 
process  during  ignition  delay. 

Molecular  dynamics  studies  can  provide  infonnation  on  the  fate  of  the  solvated  solute 
at  the  water  liquid-vapor  interface.  Figure  16  shows  that  the  small  Li+  ion  is  able  to  keep  its 
hydration  shell  intact  when  it  is  moved  to  the  liquid-vapor  interface.  The  ability  of  the  larger 
I'  ion  to  do  this  is  significantly  diminished.  In  contrast  with  the  case  of  charged  particles, 
there  is  depletion  in  the  number  of  water  molecules  in  the  first  hydration  shell  for  neutral 
atoms.  The  second  hydration  shell  is  significantly  diminished  as  the  solute  is  transferred  to 
the  interface  [Benjamin,  2015].  Although  not  thoroughly  discussed  the  transfer  of  ions  from 
the  liquid  phase  to  the  gas  phase  is  probably  possible. 


124 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


Z(A) 


Figure  16.  Relative  (to  the  bulk)  number  of  water  molecules  in  the  first  solvation  shell 
of  different  ions  and  neutral  atoms  as  a  function  of  the  distance  from  the  liquid-vapor 
interface.  Z=0  is  the  location  of  the  Gibbs  dividing  surface. 

10.  Reactions  dynamics  at  liquid  interfaces 


The  interface  between  a  liquid  and  a  second  phase  (vapor,  another  immiscible  liquid, 

and  solid)  is  a  narrow  and  a  few  nanometer  thick  highly  anisotropic  region,  where  chemical 

reactions  may  take  place.  Benjamin  presented  a  recent  review  on  this  topic  [Benjamin,  2015], 

which  was  mostly  devoted  to  the  reactions  in  the  interface  of  two  immiscible  electrolyte 
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solutions  (IESs).  For  hypergolic  ignition,  which  results  from  liquid-liquid  contact,  this 
interface  exists  even  if  the  liquids  are  mixed  in  a  turbulent  manner.  The  characteristic 
molecular  structure  and  dynamics  of  the  liquid  interfaces  are  expected  to  influence  the  rate  of 
interfacial  chemical  reactions.  This  assumption  is  valid  both  for  non-ionic/oxidant  and 
ionic/oxidant  hypergolic  systems.  In  fact,  the  miscibility  or  immiscibility  of  the  liquid 
reactants  has  not  been  thoroughly  discussed  in  the  frame  of  hypergolic  systems.  If  both 
liquids  are  miscible,  then  the  reactions  will  take  place  in  the  bulk.  If  they  are  immiscible,  then 
reactions  at  the  interface  will  be  important.  On  the  other  hand,  it  could  well  be  that  these 
liquids  are  partially  miscible  depending  on  the  nature  of  the  reactants  further  increasing  the 
complexity  of  the  problem  as  discussed  below. 


a.  Ion  transfer  (IT)  reactions  across  a  liquid-liquid  interface  between  immiscible 
electrolyte  solutions  (IESs) 


These  transfers  may  be  important  because  they  could  occur  within  ILs/oxidant  systems 
when  the  oxidant  is  a  strong  acid  such  as  nitric  acid.  The  transfer  of  ions  shares  important 
characteristics  with  typical  charge  transfer  reactions. 


b.  Electron  transfer  (ET)  reactions  at  liquid-liquid  interface 

c.  Interfacial  Sn2  reactions 


1 1 .  Interpretation  of  non-ionic/oxidant  hypergolic  ignition 
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This  part  is  certainly  of  interest,  but  here  the  focused  is  on  ILs/oxidant  systems,  and  it 


will  therefore  be  considered  as  future  work. 


12.  Interpretation  of  IL/oxidant  hypergolic  ignition 


In  the  case  of  nitric  acid  as  the  oxidant,  the  species  that  reacts  either  in  the  bulk  IL  or 
at  the  liquid-liquid  interface  or  at  the  gas-liquid  interface  is  either  HNO3  (HONCb)  or  H+  and 
NO3".  The  IL  reactant  are  presumably  ions,  which  can  be  represented  as  A+  +  B‘.  One  task 
would  be  to  ascertain  what  happens  in  the  liquid  phase  and  simultaneously  in  the  gas  phase  by 
developing  a  detailed  chemical  kinetics. 

13.  Conclusions  and  future  works 


In  part  two,  our  literature  investigation  tends  to  show  that  reactions  at  gas-liquid  and 
liquid-liquid  interfaces  may  play  an  important  role  during  hypergolic  ignition  (see  Figure  17). 


127 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


a) 

b) 

c) 

d) 

Ionic  Liquid  \ 

\ 

V 

\ 

Oxidant  - »bd  bd  ■ 

t  mmamm  m.* 

Schneidei  et  al..  2008,  Energy  and  Fuels,  Vol.  22.  p.  2871-2872 


Ignition 


Figure  17.  Visualization  of  IL/oxidant  reactivity  showing  that  the  interface  plays  a  major  role 
in  hypergolic  ignition. 


The  effect  of  solvation  on  thermochemistry  and  on  kinetics  was  examined.  A  detailed 
chemical  kinetics  model  both  in  the  liquid  phase  and  in  the  gas  phase  to  explain  hypergolic 
ignition  remains  to  be  defined. 
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Part  three: 


DETAILED  CHEMICAL  KINETIC  MODELS  FOR 
HYPERGOLIC  IGNITION  OF  ILs  BASED  SYSTEMS 
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1.  Introduction 


In  the  frame  for  searching  new  hypergolic  (auto-igniting)  propellant  systems,  it  is 
necessary  to  discuss  their  reactions  in  both  the  liquid  and  gas  phases.  In  practical  hypergolic 
applications,  propellant  fuels  and  oxidants  are  liquids  and  upon  contact,  reactions  occur  in  the 
liquid  phase  and  also  in  the  gas  phase  above  the  liquids.  One  point  which  seems  clear  is  that 
ignition  occurs  in  the  gas  phase  even  if  the  reactants  come  into  contact  as  liquids.  The  issue 
of  interest  is  related  to  the  understanding  of  the  qualitative  and  quantitative  compositions  of 
both  phases  when  ignition  occurs.  In  the  frame  of  part  two,  it  appears  that  the  chemical 
reaction  dynamics  for  liquid  hydrazines/  liquid  NTO  (or  liquid  HN03)  hypergolic  ignition  is 
somehow  elucidated,  even  if  some  points  remain  to  be  worked  out  or  discussed. 

Daimon  et  al.  [Daimon,  1984],  Tanaka  et  al.  [Tanaka,  1985],  and  Daimon  et  al. 
[Daimon,  1991]  observed  a  N2H4  droplet  falling  into  a  NTO  pool  (see  Figure  1)  and  a  NTO 
droplet  falling  into  a  N2H4  pool  (see  Figure  2).  In  both  cases,  they  found  that  a  layer  of  NTO 
was  vaporized  upon  contact.  Their  interpretations  of  the  visualizations  given  in  Figures  1  and 
2  are  reported  in  Figures  3  and  4,  respectively.  In  the  first  case,  ignition  is  supposed  to  be  due 
to  gas-gas  reactions  at  the  interface  of  the  two  liquids  whereas  in  the  second  case,  NTO  vapor 
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penetrates  into  the  liquid  hydrazine  and,  although  not  evoked  explicitly  by  the  authors, 
however,  a  gas-liquid  interface  is  fonned. 


(TIME) 


ENLARGED  PHOTO 


I _ i _ I 

0  5  10mm 


i _ 1 

0  10  mm 


C.  :  CONTACT  OF  DROPLET  AND  POOL  SURFACE 

D. L.;  APPEARANCE  OF  DISPERSED  LAYER 
EX.  :  OCCURRENCE  OF  EXPLOSION 


135 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


Figure  1.  High  speed  motion  pictures  (92900  frames  per  second)  of  the  contact  of  a 


N2H4  droplet  falling  into  a  NTO  pool. 
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Figure  2.  High  speed  motion  pictures  (48500  frames  per  second)  of  the  contact  of  a 


NTO  droplet  falling  into  a  hydrazine  pool. 
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Figure  3.  Interpretation  of  Tanaka  et  al.  [Tanaka,  1985]  to  explain  hypergolic  ignition 
when  a  droplet  of  hydrazine  is  falling  into  a  pool  of  liquid  NTO. 
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POOL 


Figure  4.  Interpretation  of  Tanaka  et  al.  [Tanaka,  1985]  to  explain  hypergolic  ignition 
when  a  droplet  of  NTO  is  falling  into  a  pool  of  liquid  hydrazine. 


For  N2H4  and  its  methylated  derivatives  MMH  and  UDMH  (unsymmetrical  dimethyl 
hydrazine)  and  for  two  oxidizers  NTO  and  FNA  (fuming  nitric  acid)  it  is  reported  that  the 
explosion  phenomena  observed  can  be  classified  into  three  categories  [Daimon,  1991]: 

In  the  case  of  hydrazine/NTO,  sudden  gasification  of  a  superheated  liquid  layer 
formed  at  the  boundary  of  the  two  liquids  occurs  spontaneously  and  a  detonation-like 
reaction  proceeds  in  the  reactive  mixture  produced. 

In  the  cases  of  MMH/NTO  and  UDMH/NTO,  sudden  gasification  is  caused  by  the 
shock  from  a  local  ignition  source  and  a  turbulent-combustion  reaction  proceeds  in  the 
reactive  mixture  produced. 
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In  the  cases  of  hydrazine  type  fuels/FNA,  sudden  gasification  occurs  spontaneously  as 
in  the  case  of  hydrazine/NTO,  but  it  is  not  augmented  by  a  chemical  reaction  and  in 
these  cases  the  observed  explosion  is  weak. 

In  all  these  cases,  the  vapor  layer,  which  is  formed  between  the  reactive  fuel  droplet 
and  the  liquid  pool,  plays  an  important  role  for  the  occurrence  of  the  explosion.  The 
composition  of  this  layer  (pure  oxidant,  pure  fuel  or  mixture  of  both)  is  not  clearly  discussed 
and  all  these  studies  do  not  identify  unambiguously  the  chemical  cause  of  ignition  or 
explosion.  In  the  previous  report  item  2,  the  chemical  dynamics  and  kinetics  in  the  gas  phase, 
in  the  liquid  phase,  at  the  gas-liquid  interface,  and  at  the  liquid-liquid  interface  have  been 
examined  to  get  a  deeper  understanding  of  the  hypergolic  ignition  phenomenon.  The  pure  gas 
phase  chemical  dynamics  and  kinetics  are  not  further  examined  here  because  this  reactivity 
(or  these  reactivities)  can  be  considered  as  sufficiently  well  known  even  if  this  statement  is 
somehow  contestable.  More  interesting,  if  one  considers  the  above  experimental  evidences 
and  interpretations,  the  chemical  dynamics  and  kinetics  at  gas-liquid  interface  are  certainly 
important.  As  underlined  by  Benjamin,  [Benjamin,  2015]  chemical  reactions  that  take  place  at 
the  interface  between  a  liquid  and  a  second  phase  (vapor,  another  immiscible  liquid  or  a  solid) 
are  of  major  interest.  Ordinarily,  the  gas-liquid  interface  is  generally  considered  to  play  the 
role  of  a  “gateway”  into  the  liquid,  and  is  believed  not  to  play  the  role  of  a  reaction  medium. 
However,  part  two  showed  that  the  gas-liquid  interface  is  a  distinct  region  where  often  unique 
reactions  can  occur.  Since  Daimon  et  al.  and  Tanaka  et  al.  have  previously  presented 
evidence  that  a  gas-liquid  interface  exists  just  before  ignition,  it  may  therefore  be  possible  that 
reactions  taking  place  in  this  region  are  the  cause  of  ignition. 

Nathanson  [Nathanson,  2004]  has  presented  a  review  of  experiments  perfonned  to 
clarify  the  role  of  the  surface  region  during  interfacial  and  bulk  phase  reactions.  When  a  gas 
phase  molecule  collides  with  a  liquid  surface  it  can  be  sent  back  to  the  gas  phase  through 
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inelastic  scattering  or  it  can  be  trapped,  i.e.,  its  kinetic  energy  is  absorbed  and  it  is  able  to  bind 
to  the  surface.  The  probability  of  the  trapping  event  depends  on  many  parameters:  the  kinetic 
energy  of  the  impinging  molecule,  the  mass  of  this  molecule  and  the  ones  at  the  surface,  the 
liquid  temperature,  its  surface  tension,  etc.  Experiments  have  demonstrated  that  gas-liquid 
reactions  can  occur  as  soon  as  the  gas  phase  molecule  enters  the  interfacial  region  and  this 
region  is  therefore  a  reaction  medium.  The  present  work  is  devoted  to  the  replacement  of 
hydrazines  by  ionic  liquids  (ILs)  and  therefore  the  interaction  of  gas-phase  molecule  with  a 
ionic  liquid  is  of  great  interest.  Part  two  emphasizes  experimental  results  devoted  to  the 
interaction  of  a  strong  acidic  molecule  represented  as  HX  with  a  ionic  liquid  in  fact  deuterated 
sulfuric  acid  D2S04  diluted  in  deuterated  water  D20.  This  scenario  can  be  representative  of 
what  happens  when  a  liquid  oxidant  (such  as  NTO  or  nitric  acid)  droplet  collides  with  a  liquid 
pool  or  droplet  of  ILs.  It  is  tenable  that  gaseous  NTO  or  nitric  acid  is  formed  upon  collision 
and  enter  the  liquid  according  to  reaction  scheme  described  in  Fig.  5  which  illustrates  the 
pathways  observed  after  a  gaseous  acid,  HX  collides  with  the  surface  of  a  solution  of 
deuterated  sulfuric  acid,  D2S04  in  deuterated  water,  D20. 
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HX  HX  HX  DX 


\  H  / 


Exchange 


Further  Reaction 


Figure  5.  Observed  pathways  for  an  acidic  molecule  HX  colliding  into  D2O/D2SO4 
[Nathanson,  2004]. 


The  HX  molecules  are  sent  back  to  the  gas  phase,  or  some  fraction  are  trapped  and 
then  desorb  unchanged.  The  remaining  molecules  react  in  the  acid  interfacial  region  and/or  in 
the  acid  bulk  as  intact  HX  or  as  H+  and  X'  after  dissociation.  These  processes  finally  lead  to 
the  formation  of  DX  molecules  in  the  gas  phase  once  desorbed  from  the  liquid  surface.  If  X  is 
Cl,  then  the  gas  phase  acid  is  HC1  (see  Figure  6).  Further  reactions  cited  in  Figure  5  are 
unspecified. 
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HC1  (Einc  =  100  kJ/mol) 


i  HC1  HC1  DC1 
45°  1 45  V  ^  v 


d2o/d2so4/^ 

HC1  +  D+-*-  DC1  +  H+ 


Figure  6.  Observed  pathways  for  HC1  colliding  into  D20/D2S04  [Nathanson,  2004]. 


HC1  enters  deuterated  sulfuric  acid  solution  in  two  steps:  a  fraction  of  the  impinging 
molecules  first  thermalize  upon  collision  with  the  surface  and  then  a  fraction  of  these  trapped 
molecules  enter  the  acid  as  intact  HC1  or  as  H+  and  Cl".  Experiments  show  that  at  low 
incident  energy,  nearly  all  HC1  molecules  are  trapped,  but  the  amount  of  DC1  formed 
indicates  that  only  1  in  9  HC1  molecules  dissociates  in  the  interfacial  or  bulk  regions  of  the 
acid  and  that  the  remaining  8  out  of  9  trapped  HC1  molecules  desorb  from  the  acid  without 
reacting.  These  results  are  particularly  interesting  for  ILs/acid  hypergolic  systems  if  one 
considers  that  the  impinging  molecule  is  WFNA  (a  strong  acid  as  is  HC1)  at  the  surface  of  an 
ionic  liquid  (as  is  the  D20/D2S04  solution).  Most  likely,  HNO3  will  dissociate  into  H+  and 
NO3'  ions,  which  then  react  with  the  ions  from  the  ionic  liquid  but  HNO3  can  also  react  as 
H0N02  molecule  in  its  neutral  form.  In  that  case  nitric  acid  will  be  present,  in  the  liquid  layer 
at  the  surface,  as  H+  and  NO3"  ions  and  as  H0N02  neutral  molecule  whereas  the  ionic  liquid 
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itself  is  present  in  the  layer  as  A+  and  B'  ions  (without  considering  further  separately  protic 
and  aprotic  ILs).  Therefore  it  is  possible  to  propose  a  general  chemical  detailed  kinetic 
scheme  in  this  layer.  For  this  purpose,  it  remains  to  examine  existing  reactions  between  ions 
in  the  liquid  state  and  reactions  between  ions  and  neutrals  in  the  liquid  state. 

2.  Interpretation  of  IL/oxidant  hypergolic  ignition:  detailed  chemical  kinetic 
model  in  liquid  phase 


Based  on  all  the  above  features  and  interpretations  or  conclusions,  in  the  case  of  nitric  acid  as 
the  oxidant,  the  species  that  reacts  either  in  the  bulk  IL  or  at  the  liquid-liquid  Il/nitric  acid 
interface  or  at  the  gas(nitric  acid)/IL  interface  is  either  neutral  HNO3  (HONO2)  or  H+  and 
NO3".  The  IL  reactants  are  ions,  which  are  abbreviated  hereafter  as  A+  and  B‘.  Therefore  the 
reactants  for  all  these  systems,  whatever  their  physical  states,  are  HONO2  (hereafter 
abbreviated  as  HNO3),  H+,  NO3",  A^  and  B\  Reactions  between  HONO2  and  H+  and  reactions 
between  HONO2  and  NO3'  can  be  considered  as  marginal,  even  if  they  exist.  For 
completeness,  these  reactions  are  considered  here.  This  leads  to  the  proposal  of  the  following 
detailed  chemical  kinetic  scheme. 


Reactions  between  nitric  acid  and  H+ 


HONO2  +  H+  -»•  H2  +  N03+ 
HONO2  +  H+  -»•  H2+  +  N03* 
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H0N02  +  H+  -»•  H20  +  N02+ 


H0N02  +  H+  -»•  H20+  +  N02* 


However,  none  of  these  reactions  seem  to  be  realistic.  Species  such  as  NC>3+,  for  instance, 
have  never  been  observed  in  the  liquid  phase  as  well  as  H2+  for  instance. 

Reactions  between  nitric  acid  and  N03‘ 

hono2  +  N03'  -»•  H0N02  +  N03' 
hono2  +  N03'  -»•  hoono2  +  NOT 
hono2  +  N03'  -*■  02N0N02  (N205?)  +  OH' 


In  this  list,  only  the  probable  reactions  are  reported,  i.e.  reactions  forming  stable  ions  such  as 
NO3",  OH'  and  molecules.  The  first  reaction  will  not  be  reported  in  the  detailed  chemical 
kinetic  model  because  reactants  and  products  are  the  same.  Concerning  the  two  remaining 
reactions,  existence  of  N2Os  is  assessed.  Interestingly,  the  second  reaction  forms  H-0-0-N02 
a  species  belonging  to  the  chemical  family  of  peroxides  containing  0-0  bond  which  is  known 
to  be  very  sensitive  to  heat  and  able  to  lead  to  the  0-0  bond  breaking  fonning  two  radicals 
according  to,  here, 

H00N02  ->  OH*  +  N03* 

and  then  to  ignition,  even  if  present  at  low  level.  It  is  also  possible  that  this  species  formed  in 
the  liquid  state  is  able  to  vaporize  and  decompose  in  the  gas  phase  leading  to  ignition 
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Reactions  between  nitric  acid  and  A+ 


H0N02  +  A+— >■  AH  +  N03+ 
HONO2  +  A+  -»■  AH+  +  N03* 
HONO2  +  A+  -»•  AOH  +  N02+ 
HONO2  +  A+  -»•  AOH+  +  N02* 


These  reactions,  as  well  as  the  ones  with  H+  above,  do  not  seem  very  significant. 


Reactions  between  nitric  acid  and  B 


HONO2  +  B'  — >  BH  +  NOT 
HONO2  +  B'  — >  BOH  +  N02' 
HONO2  +  B'  — >  BN02  +  OH' 
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All  these  reactions  are  particularly  interesting  and  depending  on  the  nature  of  B'  (to  be 
defined)  significant  species  can  be  formed  to  explain  ignition. 


Reactions  between  H+  and  B',  N03'  and  A+  and  between  species  formed  in  the  previous 
reactions 

H+  +  B'  — *  BH  +  other  reactions  depending  on  B’ 

A+  +  NO3'  — >  ANO3  +  other  reactions  depending  on  A+ 

N03*  +  N03*  ->  N206  (O2NOONO2) 

N03*  +  N02#  ->  N2O5  (O2NONO2) 

N03*  +  OH*  -»■  H02*  +  N02* 

OH*  +  OH*  ->  H202  (HOOH) 

N02*+N02*  ->  N2O4  (O2NNO2) 

N02‘  +  no3*  -»■  N03'  +  N02* 

OH'  +  N03*  -»•  OH*  +  N03' 

H02*  +  N02*  -»•  HOONO2 

Interestingly,  the  last  reaction  forms  H00N02  again,  the  peroxide  compound  previously 
described. 

OH*  +  N02*  ->  HONO2 
OH*  +  N02*  -»•  H02*  +  NO* 
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etc. 


An  infinity  of  reactions  may  be  written.  Unfortunately,  almost  none  of  these  reactions  have 
been  studied  in  the  liquid  phase  and  few  data  are  reported  in  compilations 


3.  Conclusions  and  future  works 


In  this  part  three,  our  literature  investigation  tends  to  show  that  reactions  at  gas-liquid  and 
liquid-liquid  interfaces  may  play  an  important  role  during  hypergolic  ignition  (see  Figure  17). 
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Schneider  et  a).,  2008.  Energy  and  Fuels.  Vol.  22,  p.  2871-2872 

Figure  17.  Visualization  of  IL/oxidant  reactivity  showing  that  the  interface  plays  a  major  role 
in  hypergolic  ignition. 


A  detailed  chemical  kinetics  model  both  in  the  liquid  phase  and  in  the  gas  phase  to  explain 
hypergolic  ignition  remains  to  be  defined  explicitly.  Even  if  it  is  possible  to  write  such  a 
detailed  chemical  kinetic  model  and  therefore  to  know  what  are  the  species  formed  in  the 
liquid  phase,  it  is  certainly  difficult  to  assess  quantitatively  their  presence  in  the  gas  phase. 
According  to  this  study  it  seems  that  some  peroxy  compounds  may  be  formed  in  the  liquid 
phase  and  that  these  species  may  be  present  in  the  gas  phase  where  their  thermal 
decompositions  are  able  to  induce  ignition. 
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